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E INTRODUÇÃO 
n CONTROLE ANALÓGICO 


Geralmente é uma malha de controle com realimentacáo negativa. 


e Busca atender especificações para o período transitório e para o 


regime permanente. 


e Necessidade de sensor para medir a variável controlada e comparar 


com um sinal de referéncia. 


e Acáo de controle depende de um sinal de erro entre a variável de 


interesse e o sinal de referéncia. 
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E: INTRODUÇÃO 
= CONTROLE ANALÓGICO 


e Manipula o sinal de erro. 


e Lida com variáveis e sistemas: 
e Elétricos 
e Fluidos 
e Pneumáticos 


e Mecânicos 


e Sistemas com entradas e saídas analógicas: 


e Valores definidos em um intervalo contínuo e com uma faixa de 
amplitude continua. 
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E: INTRODUÇÃO 
e CONTROLE DIGITAL 


e Nas ultimas décadas os controladores digitais estão substituindo 


gradativamente os controladores analógicos. 


Entradas e saídas de um controlador digital são definidas em 
intervalos de tempo discretos. 


e Exemplos de controladores digitais: 
e Circuitos Digitais. 
e Computadores Digitais. 


e Microprocessadores e Microcontroladores. 
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INTRODUÇÃO 
= CONTROLE ANALÓGICO X CONTROLE DIGITAL 


e Qual o melhor tipo de controle”? 


e a) Controladores que monitoram continuamente uma variável. 
e b) Controladores que se baseiam em sinais amostrados de uma variável. 


e Variáveis de controle que mudam continuamente sáo melhores controladas do 


que as que mudam periodicamente. 


e Supondo condições idênticas para o controle analógico e o controle digital: 


e Controle Analógico é superior ao Controle Digital 


e Então, qual o motivo da popularidade do Controle Digital nas últimas 


décadas? 
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E: INTRODUCAO 
ш” POR QUE CONTROLE DIGITAL? 


e Principais vantagens que explicam a popularidade do Controle Digital: 


e Precisáo 

e Erros de Implementacáo 
e Flexibilidade 

e Velocidade 


e Custo 
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E: INTRODUÇÃO 
e PRECISÃO 


e Sinais Digitais são representados usando “0” e “1º com 16 bits ou mais 


para representar um único número. 


e Erro muito menor do que os que aparecem em sinais analógicos onde o 


ruído está sempre presente. 
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E INTRODUÇÃO 
т ERROS DE IMPLEMENTAÇÃO 


e Processamento digital de sinais de controle envolve adição e 
multiplicação de valores numéricos armazenados: 
e Os erros que resultam da representação digital e manipulação 


aritmética são desprezíveis. 
e Processamento de sinais analógicos é realizado utilizando componentes 


como resistores e capacitores, com valores reais que variam 


significativamente dos valores nominais de projeto. 
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E: INTRODUÇÃO 
үө FLEXIBILIDADE 


e Um controlador analógico é difícil de se modificar ou reformular uma vez 
implementado em hardware: 


e Circuitos eletrónicos com amplificadores operacionais. 


e Geralmente sáo estruturas de controle simples. 


e Um controlador digital é implementado no firmware ou software: 
e Sua alteração é possível sem uma substituição completa do controlador 
original. 
e Estruturas náo sáo limitadas a formas simples, que sáo normalmente 
utilizadas no controle analógico: 
Estruturas de controle complexas envolvendo algumas operações 


aritméticas extras podem ser “facilmente” implementadas. 
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E: INTRODUÇÃO 
ee VELOCIDADE E CUSTO 


e Hardware mais “sofisticados” permitem curtos períodos de 


amostragem (taxas de amostragem rápidas). 


e Com a curtos períodos de amostragem, controladores digitais podem 


monitorar as variáveis controladas quase continuamente. 


e Avanços na tecnologia VLSI (Very Large Scale Integration) fornecem 
circuitos integrados melhores, mais rápidos e confiáveis, a preços 
mais baixos. 

e Relação custo x benefício permite a sua utilização mesmo em 


aplicações de pequeno porte. 
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E INTRODUÇÃO 
= ESTRUTURA DE UM SISTEMA DE CONTROLE DIGITAL 


e Para controlar um sistema físico ou processo utilizando um controlador digital, é 


necessário que o controlador: 
e 1) Receba as medições do sistema. 


e 2) Processe essas informações. 


e 3) Envie o sinal de controle para o atuador, para aplicar a ação de controle. 
Na grande maioria das aplicações, planta e atuador são sistemas analógicos. 


e Controlador e sistema a ser controlado não falam a mesma língua: 
e Necessário “traduzir” essa informação. 


e Conversores A/D e D/A. 
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E: INTRODUÇÃO 
= ESTRUTURA DE UM SISTEMA DE CONTROLE DIGITAL 


e Conversor Digital-Analógico (Digital-to-analog converter - DAC): 
e Traduz os sinais do controlador digital para a linguagem do processo 


analógico. 


e Conversor Analógico-Digital (Analog-to-digital converter - ADC): 
e Faz o processo inverso, traduzindo os sinais do processo analógico 


para a linguagem do controlador digital. 


e Sensor: 


e Necessário para monitorar a variável controlada que é utilizada na 


realimentação. 
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INTRODUÇÃO 
= ESTRUTURA DE UM SISTEMA DE CONTROLE DIGITAL 


Entrada (Referência) 


Computador ou 
microprocessador 


Variável Controlada 


Atuador e 
Processo 
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E INTRODUÇÃO 
= ESTRUTURA DE UM SISTEMA DE CONTROLE DIGITAL 


e Outras configurações de controle: 


e Várias entradas de referência, variáveis controladas, e malhas 


fechadas. 


e Existência de malha interna com controle digital ou analógico. 
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E INTRODUÇÃO 
= EXEMPLOS DE SISTEMAS DE CONTROLE DIGITAL 


e Sistema de Medicação Automática em Malha Fechada 
e Carson, E. R.; T. Deutsch. A spectrum of approaches for controlling 


diabetes in Control System Magazine, 12(6):25-31, 1992. 
e Controle do Motor à Jato de uma Aeronave 


e Para atingir o alto desempenho necessário para as aeronaves, 


empregam-se estratégias de controle sofisticados via computador. 
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E INTRODUÇÃO 
= EXEMPLOS DE SISTEMAS DE CONTROLE DIGITAL 


e Controle de um manipulador Robótico 
e Os manipuladores robóticos são capazes de executar tarefas repetitivas em 


velocidades e exatidões melhores do que os operadores humanos. 


e Para executar suas tarefas com precisão e confiabilidade, posição e 


velocidade da mão do manipulador são controladas digitalmente. 


e Cada movimento ou grau de liberdade (DOF - degree of freedom) do 
manipulador é posicionado utilizando um sistema de controle independente de 


posição. 
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Ë INTRODUCAO 
7"* SISTEMA DE MEDICAÇÃO AUTOMÁTICA EM MALHA FECHADA 


Computer 
| Regulated 
Drug 

or Nutrient 


Reference Regulated 


Blood Drug 


Level or Nutrient 
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INTRODUCAO 
CONTROLE DE UM MOTOR A JATO DE UMA AERONAVE 
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E: INTRODUÇÃO 
нїттлонр С 


ONTROLE DE UM MANIPULADOR ROBÓTICO 


NET TUTDFED 


E INTRODUÇÃO 
е DIAGRAMA DE BLOCOS: CONTROLADORES DIGITAL X ANALÓGICO 


Controle contínuo 
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E: INTRODUÇÃO 
"= DIAGRAMA DE BLOCOS: CONTROLADORES DIGITAL X ANALÓGICO 


Di gi tal Controller 
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E MODELO DE TEMPO DISCRETO PARA 
CONTROLE DIGITAL 


e Principal objetivo de um sistema de controle: 
e Controlar um sistema físico (planta) para que ele se comporte de acordo 
com determinadas especificações de desempenho. 
e Geralmente: 
e Planta > contínua (analógica). 


e Controle desejado (referência) > constante por partes. 


e Ação de controle > atualizada periodicamente. 


e Uma boa aproximação permite que um modelo simples possa ser obtido. 


e Neste caso, é possível representar esse tipo de sistema através de um 
modelo de tempo discreto. 
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E: EXEMPLO: SISTEMA DE CONTROLE DE 
т= NÍVEL EM UM TANQUE 


e Exemplo 1: Deseja-se manter o nível de fluido no tanque, ajustando-se o 
fluxo de entrada de líquido 


e H > nível do fluido em regime permanente. 

e h > perturbação no nível do seu valor nominal. 

e Q > fluxo em regime permanente através do tanque. 

e qi > perturbação no fluxo de entrada do valor nominal. 
e qo > perturbação no fluxo de saída do valor nominal. 
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E SISTEMA DE CONTROLE DE NÍVEL EM UM TANQUE 
e SOLUÇÃO 


e Obter um modelo matemático contínuo (analógico) do tanque e, a partir 


desse modelo, obter um modelo de tempo discreto para o sistema. 


e Supor que: 
e Entrada é constante por partes (fluxo de entrada q;). 


e Saída é a perturbação em torno do nível nominal (h). 


e Observação: Embora o sistema seja não-linear, um modelo linear pode 
descrever satisfatoriamente o sistema supondo que o nível do fluido é 


regulado em torno de um valor constante. 
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E SISTEMA DE CONTROLE DE NÍVEL EM UM TANQUE 
aee SOLUÇÃO 


e O modelo linearizado para a válvula de saída é análogo a uma 
resistência elétrica: 


e onde: 


e h > perturbação no nível nominal do tanque. 


e qo > perturbação na saída do tanque de um nível nominal Q. 


e R > resistência do fluido da válvula. 
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E SISTEMA DE CONTROLE DE NÍVEL EM UM TANQUE 
e SOLUÇÃO 


e Supondo que o fluido é incompressível, o princípio da conservação de 
massa reduz-se ao equilíbrio volumétrico: 
e taxa de aumento do volume de fluido = taxa de volume de fluido entrando 
- taxa de volume de fluido saindo. 


dC(h + H) 


22 = (qi + Q) — (qo + Q) 


e onde: 
e C > área do tanque (capacitáncia do fluido). 


e Lembrando que H é constante (derivada nula) e que o termo Q se 
cancela, sobram na equacáo apenas termos relacionados com a 


perturbacáo. 
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E SISTEMA DE CONTROLE DE NÍVEL EM UM TANQUE 
aee SOLUÇÃO 


e Substituindo o modelo linearizado da válvula de saída na equação de 


balanço volumétrico do fluido, tem-se: 


e onde т = RC é a constante de tempo do fluido para o tanque. 


e A solução da equação diferencial é dada por: 


1 t 
h(t) = e 0/Th(to) + = | et qd 


to 
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E SISTEMA DE CONTROLE DE NÍVEL EM UM TANQUE 
e SOLUÇÃO 


e Considerando q; constante em cada período de amostragem T: 


qi(t) = qi(k) = constante para t no intervalo [kT, (k + 1)T]. 


e Então, pode-se resolver a equação continua para qualquer periodo de 


amostragem e obter o modelo de tempo discreto para o sistema: 


h(k +1) =e T/*h(k) + R|1—e77/7]q,(k) 


e onde as variáveis nos instantes kT são denotadas por simplicidade 


apenas por k. 
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E SISTEMA DE CONTROLE DE NÍVEL EM UM TANQUE 
e SOLUÇÃO 


e Observações importantes: 


e O modelo de tempo discreto obtido neste exemplo é conhecido 


como equação a diferenças. 


e Para sistemas contínuos invariantes no tempo, a equação a 


diferenças também será linear e invariante no tempo. 
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- EQUACOES A DIFERENCAS 


e Equações Diferenciais: 


e Representam sistemas contínuos no tempo. 


e Equações a Diferenças: 
e De forma análoga, representam sistemas discretos no tempo. 
e O valor da saída no instante de tempo y(kT) é função de todas as 


entradas e saidas anteriores à kT. 


e Uma equação a diferenças aparece em problemas onde a variável 
independente, geralmente o tempo, tem a possibilidade de assumir um 


conjunto discreto de valores possíveis. 
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- EQUACOES A DIFERENCAS 


e Tipos: 
e 1) Equacáo a diferencas náo-linear de ordem n: 


y(k +n) = Рук -n— 1), y(k - n — 2), ..., y(k + 1), y(k),u(k + n), 
u(k + тп — 1), ..,u(k + 1), u(k)] 


e 2) Equacáo a diferencas linear de ordem n: 
y(k+n)+ a, 4y(k -n — 1) t: a4y(k +1) + agy(k) 
= bíu(k +n) + b, 4u(k - n — 1) t + b4u(k + 1) + bgu(k) 
e Observação: 
e Quando os coeficientes a;,b;, i = 0,1,2,... são constantes, diz-se que a 
equação a diferenças é linear e invariante no tempo (LTI). 
e Para os casos onde a função de excitação u(k) é nula, a equação a 


diferenças é dita ser homogênea. 
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E: EQUACÓES A DIFERENCAS 
т EXEMPLOS 


e Para as equações a diferenças abaixo, determine: 
e sua ordem. 
e seu tipo (linear ou não-linear; variante ou invariante no tempo; 


homogénea ou náo-homogénea). 
1. y(k + 2) + 0.8y(k + 1) + 0.07y(k) = u(k) 
2. y(k + 4) + sin(0.4k)y(k + 1) + 0.3y(k) = 0 


3. y(k + 1) = —0.1 y? (k) 
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EQUAÇÕES A DIFERENÇAS 
E SOLUÇÃO: EXEMPLO(1) 


1. y(k + 2) + 0.8y(k + 1) + 0.07y(k) = u(k) 


e Ordem: 


e Equacáo de 22 ordem. 


e Coeficientes constantes: 


e Equacáo linear e invariante no tempo. 


e Existe uma entrada forçada u(k): 


e Equação não-homogênea 
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EQUAÇÕES A DIFERENÇAS 
Ё. SOLUÇÃO: EXEMPLO(2) 


2. y(k + 4) + sin(0.4k)y(k + 1) + 0.3y(k) = 0 


e Ordem: 


e Equacáo de 42 ordem. 


e O 22 termo é dependente do tempo, mas todos os termos sáo 
lineares: 


e Equacáo linear e variante no tempo. 


e Náo existe uma entrada forcada: 


e Equacáo homogénea. 
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EQUAÇÕES A DIFERENÇAS 
Ё. SOLUÇÃO: EXEMPLO(3) 


3. y(k + 1) = —0.1 y?(k) 


e Ordem: 


e Equacáo de 12 ordem. 


e O lado direito da equacáo é claramente náo-linear, mas náo apresenta 
termos que dependam do tempo explicitamente: 


e Equação não-linear e invariante no tempo. 


e Não existe uma entrada forçada: 


e Equação homogênea. 
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E: EQUACOES A DIFERENCAS 
a SOLUÇÃO: EXEMPLOS 


1. y(k + 2) + 0.8y(k + 1) + 0.07y(k) = u(k) 
2. y(k + 4) + sin(0.4k)y(k + 1) + 0.3y(k) = 0 
3. y(k + 1) = —0.1 y? (k) 


e Equacáo 1: 


e 22 ordem, LTI, não-homogênea. 


e Equação 2: 


e 42 ordem, linear variante no tempo, homogênea. 


e Equação 3: 


O e 12 ordem, não-linear, invariante no tempo, homogênea. 
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Е. _ TRANSFORMADA 2 


e A transformada Z е uma ferramenta importante рага апаііѕе e projeto 


de sistemas discretos. 


e Com a transformada Z, é possível converter equações a diferenças 
lineares e invariantes no tempo para equações algébricas, além de 


transformar operações de convolução em multiplicação. 


e Pode-se dizer que a transformada Z tem a mesma importância para 
sistemas discretos que a transformada de Laplace tem para sistemas 


contínuos. 
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E: TRANSFORMADA Z 
= DEFINIÇÃO 1 


e Dada uma sequência causal (nula para instantes de tempo negativos) 


do tipo: 
(uo, Us, Us, ..., Up, e] 


e Asua transformada Z é definida como: 


U(Z) = uo + uzt + uz 2 +... +lukZ É 


00 


U(z) = > uyz * 


k=0 


e Observação: A variável z”! pode ser entendida como um operador de 


atraso de um intervalo de amostragem. 
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E: TRANSFORMADA Z 
= DEFINIÇÃO 2 


e Dada a seguinte representação de um trem de impulsos de um sinal 


discreto: 
ut) = ugó(t) + u4ó(t — T) + uó6(t — 2T) +--+ и, ó (t — КТ) 
u'(t) = >. uxÓ(t — КТ) 
k=0 


e A sua transformada de Laplace pode ser escrita como: 


U*(s) = ug +u,e ST + изе 7T +... + ирет T +... 


005) = У unes) 


k=0 
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E: TRANSFORMADA Z 
= DEFINIÇÃO 2 


e Aplicando-se a seguinte igualdade: 


Ae 


e Chega-se a expressão apresentada anteriormente: 


U(zZ) = uo + uz !+u,z2+-+upz* 


00 


U(z) = > uz ^ 


k=0 
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E TRANSFORMADA Z 
" EXEMPLO 


e Para a sequência abaixo, obtenha sua transformada Z: 


(ux)? = (1,3,2,0,4,0,0,0, ... ) 


e Aplicando a Definição 1: 


U(z) =1+3z 1+2z 2 + 4z 4 
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E: TRANSFORMADA Z 
: CONCLUSÓES 


e Embora a mesma transformada Z possa ser obtida pelas duas 


maneiras mostradas anteriormente, cada uma tem suas vantagens e 
desvantagens. 


e Definicáo 1: 


e Evita o uso de impulsos e da transformada de Laplace. 


e Definicáo 2: 
e Permite tratar Z como uma variável complexa. 


e Possibilidade de usar as diversas propriedades da transformada de 


Laplace (linearidade, deslocamento no tempo, etc.) 


Eng. de Controle e Automação – Rejane de Barros Araújo/André Ferreira 


E: TRANSFORMADA Z 
ss CONCLUSOES 


e Fica claro que é possível utilizar a transformada de Laplace para 


estudar sistemas discretos, sistemas contínuos e sistemas mistos. 
e Entretanto, a transformada Z oferece uma grande simplificacáo na 


notacáo para sistemas discretos, simplificando bastante o esforco 


para análise e projeto desses sistemas. 
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E _ TRANSFORMADA Z 


IDENTIDADES USADAS PARA DERIVAR TRANSFORMADAS Z DE FUNÇÕES BÁSICAS 


1— 
B ad^ = а=1 


2.138 ак = — lal < 1 
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E TRANSFORMADA Z 
= IMPULSO UNITÁRIO 


e Considerando o impulso discreto: 


E k = 0 


u(k) = у= Ë k = 0 


б) Eng. de Controle e Automação – Rejane de Barros Araújo/André Ferreira 


E TRANSFORMADA Z 
= IMPULSO UNITÁRIO 


e Aplicando a Definicáo 1, a transformada Z do impulso unitário é dada 
por: 

U(z) = 1 

e A partir da Definição 2, pode-se fazer u“(t) = ó(t), o que leva à 


seguinte transformada de Laplace: 
U*(s) = 1 


e A substituição por z=e* não causa nenhum efeito. Assim, a 
transformada Z obtida usando a Definição 2 apresenta o mesmo 
resultado da Definição 1. 
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E: TRANSFORMADA Z 
e DEGRAU UNITÁRIO 


e Considere a seguinte sequência: 


(Jo = 11,1,1,1,1,1,1,1,...} 
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E TRANSFORMADA Z 
= DEGRAU UNITÁRIO 


e Usando a Definição 1: 


U(z) = 1 „+71 +77? + 3773 + ... + Zk +... 
U(z) = qa 
k=0 


e Aplicando-se a identidade (11), tem-se a seguinte forma fechada para a 


transformada Z do degrau unitário discreto: 
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E TRANSFORMADA Z 
т” DEGRAU UNITÁRIO 


e É importante notar que a identidade (II) é válida apenas para || < 1. Isto 
significa que a transformada Z obtida apresenta uma região de 
convergência fora dessa região que não é válida. 


e А região de convergência deve ser dada claramente tomando-se a 
transformada nos dois lados de funções que apresentam valores não- 
nulos para tempos negativos. 


e Porém, para funções temporais que são nulas para o tempo negativo, 
pode-se aplicar a transformada Z em apenas um dos lados e estender a 
região de convergência, com a transformada Z obtida sendo válida 
em todo o plano-Z. 
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E TRANSFORMADA Z 
nd EXPONENCIAL DISCRETA 


e Seja: 
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E: TRANSFORMADA Z 
=" EXPONENCIAL DISCRETA 


e Aplicando-se a Definicáo 1: 
О(2) = 1+ а2 l+a2z +... + akz Kk +... = > (9/2)* 
k=0 


e Usando-se a identidade (Il), pode-se escrever: 


Z 


1 
MARE z-a 


e Da mesma forma que para o degrau, pode-se estender a validade da 
transformada Z para todo o plano-Z, já que a função é nula para valores 


negativos. 
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E TRANSFORMADA Z 
i PROPRIEDADES 


e Atransformada Z pode ser obtida a partir da transformada de Laplace, 


como mostrado na Definicáo 2. 


e Desta forma, a transformada Z compartilha de uma série de 


propriedades que a transformada de Laplace apresenta. 


Ө Eng. de Controle e Automação – Rejane de Barros Araújo/André Ferreira 


E TRANSFORMADA Z 
" PROPRIEDADES 


e Linearidade: 


Z(afi(k) + P f, (k); = аР, (2) + BF (2) 


e Exemplo 2: 


e achar a transformada Z da sequéncia causal: 

FU) = 2xA(K) + 46(K), k = 0,1,2, ... 
e Solucáo: 

F(z) = Z(2x1(k) + 48(k)} = 2Z(1(k)) + 4Z1[6(k)) 


2Z 6z — 4 
F(z) = ——+4= 


z—1 z = 1 
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E TRANSFORMADA Z 
" PROPRIEDADES 


e Atraso no tempo: 


Zif (k —n)) = z "F(z) 


e Exemplo 3: 


e achar a transformada Z da sequéncia causal: 


CE 2 9... 
caso contrário 


“Já, 
к= Ë: 
e Solucáo: 


F(z) = Z(4x1(k — 2)) = 47 “Z(1(k)) 


4z 4 
F(z) = z 2 


z—1 z(z-1) 
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E TRANSFORMADA Z 
i PROPRIEDADES 


e Avanco no tempo: 


ZU (k + 1); = zF(z) — zf (0) 


Z(f(k + 1)) = z"F(z) - z"f(0) - z" fA) – р (п - 1) 


e Prova: 
e Aplicando a Definicáo 1 para uma funcáo discreta no tempo 


avancada em um intervalo de amostragem, tem-se: 


Z(f(k--1) = > f(k+ 1)z =z > f(k+ 1)z 0+D 
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E TRANSFORMADA Z 
" PROPRIEDADES 


e Adicionando-se e subtraindo-se a condição inicial: 


Z(f(k + 1)) = ro + > f(k + pe - го) 
k=0 


e Alterando-se o índice do somatório para m=Hk+1, pode-se 


reescrever a transformada Z como: 


ZU(k+1)) = 2 |> fen = rol = zF (z) — zf(0) 
m=0 
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E TRANSFORMADA Z 
" PROPRIEDADES 


e Exemplo 4: Utilizando a propriedade de avanco no tempo, ache a 


transformada-Z da sequéncia causal: 


(f(k)) = (48,16, ...) 


e Solucáo: 


e a sequência pode ser escrita na forma: 
für e2* = g(k +2), k= 0,1,2,... 
e onde g(k) é a função exponencial g(k) = 2*, k = 0,1,2, ... 


e utilizando a propriedade de avanco no tempo: 
47 


F(z) = zºG(z) —z^g(0) — zg(1) = z = : 


A = 
= Z“ == = 
== 2 Z—2 
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E TRANSFORMADA Z 
" PROPRIEDADES 


e Exemplo 4: solução mais simples: 


e Reescrever a sequéncia na forma: 


(f(k)) = 4(1,24, ...) = {4.25} 


e Aplicar a propriedade de linearidade da transformada Z: 


47 
z — 2 


F(z) = 
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E: TRANSFORMADA Z 
ges CONVOLUCAO (LEMBRETE!) 


e Convolucáo é um operador que, a partir de duas funcóes, produz uma 
terceira. O conceito de convolucáo está ligado ao de média móvel, 


sendo crucial no estudo de sistemas lineares invariantes no tempo. 


e A notação para a convolução de fe g é f » g. Еа é definida como a 
integral do produto de uma das funcóes com uma cópia invertida, com 
relacáo a um determinado plano, da outra. A funcáo resultante depende 


do valor deste deslocamento. 


e A definição de convolucáo pode ser aplicada tanto para funções de 


domínio contínuo quanto para funções de domínio discreto. 
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E TRANSFORMADA Z 
" PROPRIEDADES 


e Convolucáo Discreta: 


ZAR) * fo(k)) = zo. -o fii) f2 (k — i)! = Е, (2)F, (z) 


e Prova: 
e Suponha que o resultado da convolução seja dado por fy(k)), de tal forma 


“E YC) = у(0) +y()z 1 + у(2)272 + + y(k)z + 


Y) = rr = Y VORE- | 


k=0 | í=0 


e com f(k —i) 20, i > К (sequência causal) 
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E TRANSFORMADA Z 
ni PROPRIEDADES 


e Assim, pode-se colocar Y (z) na forma: 
YG) => AO > folk — pe 
1=0 k-i 


e Alterando o índice do somatório de k para j — (k — i), tem-se: 


Y(z) = Y AG У ечи 


i=0 j=0 


Y(z) = Dio fili)z” [izo fz |-F (DF, (2) 
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E TRANSFORMADA Z 
" PROPRIEDADES 


e Exemplo 5: 
e Ache a transformada Z da convolução entre dois degraus 


amostrados. 


e Solucáo: 
e Utilizando o teorema da convolução, basta multiplicar as 


transformadas Z de cada uma das funções: 


F(z) = S) 1) E G : ) 
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E TRANSFORMADA Z 
ni PROPRIEDADES 


e Multiplicação por um exponencial: 
Z(a *f(k)) = F (az) 
e Prova: 


Lado esquerdo = > q fuz" 


k=0 


= > fŒ (a2)-* 


k=0 


= F(az) 
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E TRANSFORMADA Z 
е PROPRIEDADES 


e Exemplo 6: 


e Obtenha a transformada Z da seguinte sequéncia exponencial: 
f(k)=e7®T, k=0,12,.. 


e Solucáo: 


e Lembrando que a transformada Z do degrau amostrado vale: 
F(z2)=(1—z 1) 1 
e Assim, f (k) pode ser reescrita como: 


f(k) = (exi, k = 0,1,2,... 
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E TRANSFORMADA Z 
" PROPRIEDADES 


e Aplicando-se a propriedade da multiplicacáo por um exponencial: 


Z 


Z((e?T)-*f(k)) = [1 — (e?7z) 1] = 


z—e- at 


e Observacáo: Lembre-se que conforme visto anteriormente: 


_ јак, k20 
109 = |9 ‚ k«0 


2 


1 
И =I 2-а 
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E TRANSFORMADA Z 
nd PROPRIEDADES 


e Diferenciação Complexa: 
d m 
Zik" f (k); = (25) F(z) 


e Prova: 
e Fazendo inicialmente m = 1: 


00 00 q 
107009) = > kfü0z-* = > fQO Е =) z- 
k=0 k=0 


die d 
Z(kf (k)) = d go |) Ро 


E TRANSFORMADA Z 
" PROPRIEDADES 


e Supondo agora que m pode assumir qualquer valor e definindo a 


seguinte sequéncia: 


Ё. (К) = К"Р(К), k= 0,1,2,... 


e Pode-se obter a seguinte transformada: 


00 00 q 
Zik fm (k)) = > k fm(k)z ~" = > fa (k ) (-z =) z 
k=0 k=0 


а\ҳ I d 
Zik fon (10) m [z NL к= (25) Fa (2) 
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E TRANSFORMADA Z 
" PROPRIEDADES 


e Substituindo E, (z), tem-se: 


d т+1 
ZÁk fm (k); = (-z =) F(z) 
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E TRANSFORMADA Z 
" PROPRIEDADES 


e Exemplo 7: 
e Achar a transformada Z da rampa amostrada, representada pela 
sequéncia: 
f(k)=k, k=0,1,2,.. 
e Solucáo: 
e Lembrando que a transformada Z do degrau amostrado vale: 


Z 


POS d 


e Pode-se reescrever f (k) como: f (k) = kx1, k = 0,1,2,... 
e Aplicando-se a propriedade da diferenciacáo complexa: 


= d 7 (z—1)—-z 
ZAkx1) = 07 Zíkx1] = (-22) Е = r) = (—z) (z E 1)? 
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m _ TRANSFORMADA Z INVERSA 


e Principal propósito da transformada Z: 
e simplificar a solução de problemas de tempo discreto. 
e torna-se necessário calcular a transformada Z inversa para obter a 


solucáo no domínio do tempo. 


e Transformada Z é análoga a transformada de Laplace: 
e uma integral complexa pode ser utilizada para se obter a transformada 
inversa. 
e integral geralmente trabalhosa para ser resolvida e raramente 


necessária em aplicações de engenharia. 
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m _ TRANSFORMADA Z INVERSA 


e Duas maneiras mais simples de se obter a transformada Z: 
e 1) Divisão longa (contínua): 
e pode fornecer qualquer número de termos da série que forem 


desejados. 


e 2) Expansão em frações parciais: 


e similar à transformada de Laplace inversa. 
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E TRANSFORMADA Z INVERSA 
- DIVISAO LONGA 


e Baseia-se na Definição 1, que relaciona uma série temporal com sua 


transformada Z diretamente. 


e Utiliza-se a divisáo longa para se obter tantos termos quanto forem 


desejados a partir da expansáo da transformada Z. 


e Depois, utilizando os coeficientes da expansão, pode-se escrever a 


sequência no tempo. 
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E TRANSFORMADA Z INVERSA 
- DIVISAO LONGA 


e Procedimento: 
e 1) Utilizando a divisáo longa, expanda F(z) como uma série para 


obter: 


FE(z)-—f,t-fz54«e—fzi-ofz* 
otf 2. : 


e 2) Atransformada Z inversa pode ser representada pela sequéncia: 


Uo fv a s — 
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E TRANSFORMADA Z INVERSA 
edi DIVISÃO LONGA 


e Exemplo 8: 
e Obtenha a transformada Z inversa da função abaixo: 


z+1 


F(z) = == 
(2) z? + 0.27 + 0.1 
e Solução: 
e 1) Divisão longa: 2+1 | 22+0,22+0,1 
-z-0,2-0,1271 z1+ 0,8 z — 0,26 z` +... 
0,8- 0,1 z! 


— 0,26 z1— 0,08 22 
0,26 2' + 0,52 22+ 0,026 z 


0,44 22+ 0,026 z? 
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E TRANSFORMADA Z INVERSA 
„й DIVISAO LONGA 


e Exemplo 8: continuação 
e Considerando apenas os 4 primeiros termos da sequência, tem- 


Se: 


F(z) = 0 + z 1 + 0.827? — 0.26z ? 


e 2) Aplicando-se a transformada inversa: 


(f) = (0,1,0.8, —0.26, ...) 
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E TRANSFORMADA Z INVERSA 
=  EXPANSAO EM FRAÇÕES PARCIAIS 


e Como a maioria das funções que representam sistemas discretos 


possuem o termo z no seu numerador, geralmente é mais conveniente 


expandir ^ ?/;, ao invés de somente F(z). 


e À expansão em frações parciais permite que se escreva a função 
original como uma soma de funções mais simples, cujas transformadas 


Z sejam funções no tempo discreto conhecidas. 


e As funções no tempo mais utilizadas, assim como suas transformadas Z 


são mostradas a seguir: 
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Discrete Timer — Z-Transform 


6e) 


z(z+1) 
(z—1y 


2(22+42+1) 


(z -1) 
5+0 z—a 
7 l-e” a 1—а* (I-a)z 
s(s+0t) (2-1)(2-а) 


(аЬ) = 
(z—a)(z—b) 


10 sin( a£) On sin(@, Ë) sin (0), ) Z 
s +0; 2-2cos(0,)7+1 
11 costo, t) s cos( a, ) z|z—cos(0, )] 


52+ 0} 2? 2с05(0,)2 +1 
Eng. de Controle e Automação – Rejane de Barros Araújo/André Ferreira 


TRANSFORMADA Z 


Continuous Laplace 
Time Transform 


e“ sin (ei) ED, 


(s-- Сод, Y +0 


e 9^! cos(at) s +C, 


(s-- a, Y +0 


cosh( Bir) 


Discrete Time* 


e ^" sin (0), i) 


er" cos(m,k) 


cosh( Bl) 


Z-Transtorm 


e C0 sin(m,)z 


e” cos(m,)z-e +” 
zlz-e*” cos(m,)] 
z^—2e*"" cos(m,)z+e En 


sinh(B)z 


z^—2cosh(f)z 1 


z|z-cosh(8)] 
z*-2cosh(Bjz+1 


*The discrete time functions are similar to but not always a sampled form of the continuous time 


functions. 


“Sampling г gives RT with a transform that is obtained by multiplying the transform of k by T. 
**"The function e "^" is obtained by setting a = e". 
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E TRANSFORMADA Z INVERSA 
ий EXPANSÃO EM FRAÇÕES PARCIAIS 


e Procedimento para obter a transformada Z inversa: 
e 1) Achar a expansão em frações parciais de ED, ou F(z). 


e 2) Obter a transformada inversa f(k) utilizando a tabela de 


transformadas Z. 


e Deve-se considerar três tipos de funções F(z) no domínio Z: 
e a) funções com polos reais simples e desiguais. 


e b) funções com polos reais simples e polos complexos conjugados. 


e C) funções com polos iguais. 
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É. TRANSFORMADA Z INVERSA 


EXPANSÃO EM FRAÇÕES PARCIAIS — POLOS REAIS SIMPLES 


e Utiliza-se o método dos resíduos. 


e O resíduo de uma função F(z) para um polo simples z; é dado por: 


A; = (z — z;)F(z) 
292} 
n 


А 


e onde A; é o i-ésimo termo da expansão F(z) = у "m 
x i 
1=1 


e Observação: Сото a maioria das transformações па tabela de 


transformadas Z inclui um termo z no numerador, é conveniente em 


muitos casos fazer a expansão de * (22/2, como será mostrado а seguir. 
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Ë TRANSFORMADA Z INVERSA 


EXPANSÃO EM FRAÇÕES PARCIAIS — POLOS REAIS SIMPLES 


z+1 


e Exemplo 9: Obter a transformada Z inversa de F(z) = 7240374002 


e Solução: 
e a) Expansão em frações parciais de ua 


F(z) |. z+1 . B E C 
z | z(z24-0.32-0.02 z z+01 z-02 


e Os resíduos (A, B e C) são dados por: 


Fo ÃO Е Е(2) | 1-02 ` 
кес O =н С = @+02)— ШИ 0202010) ^ 
F(z) 1- 0.1 

B = (z + 0.1) —90 


o E0DOL — 
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Ë TRANSFORMADA Z INVERSA 


EXPANSÃO EM FRAÇÕES PARCIAIS — POLOS REAIS SIMPLES 
Exemplo 9: continuaçao 


Desta forma, multiplicando-se na, por z, tem-se: 


И у= 907 407 
DC x z401 z+02 


Utilizando-se a tabela de transformadas Z, tem-se: 


f(k) = asa — 90(—0.1)* + 40(—0.2)*, 


k > 0 
0, 


k «0 


Mas como f(0) = 0, f (k) se reduz a: 


_ |—90(—0.1)¥ + 40(-0.2)*, k>0 
8 | 0, k «0 
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É. TRANSFORMADA Z INVERSA 


EXPANSÃO EM FRAÇÕES PARCIAIS — POLOS REAIS SIMPLES 


e Exemplo 9: continuação 
e Solução: 
e b) Expansão em frações parciais de F(z): 


z+1 A B 


F — D— T ——— 
G)= 07730374002 2+01 2+02 


e Os resíduos (A e B) podem ser calculados usando as expressões: 


1 — 0.1 
А = (z + oiFG| _ = = 9 


В = (z + DE] = == -8 


Ө Eng. de Controle e Automação – Rejane de Barros Araújo/André Ferreira 


Ë TRANSFORMADA Z INVERSA 


EXPANSÃO EM FRAÇÕES PARCIAIS — POLOS REAIS SIMPLES 


e Exemplo 9: continuação 
e Solução: 
e Desta forma, a expansão de F(z) vale: 


9 8 
z+01 2+0.2 


F(z) = 


e Utilizando-se a tabela de transformadas Z, tem-se: 


F(z) = 92 E 8Z ЕЕ 
^7 ROL 24027 


e Pode-se obter a transformada inversa utilizando a propriedade de 


atraso no tempo, com f(k) sendo dada por: 


 .|9(-0.1)*1—8(-02)* 1 k 20 
fuss | 0, k «0 
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É. TRANSFORMADA Z INVERSA 


EXPANSÃO EM FRAÇÕES PARCIAIS — POLOS REAIS SIMPLES 


e Recapitulando: 


- = k = k 
° Utilizando ^/;: f(k) zi 90(—0.1)*  40(-0.2*, k>0 


0, k «0 


9(—0.1)*1 — 8(-0.2)* "1, k>0 


e Utilizando F(z): f(k) = | 0 k«0 


e Variando-se k nas duas expressões, verifica-se que elas são 


completamente equivalentes. 


e Observacáo: Embora seja mais simples e direto obter a expansáo em 
fracóes parciais sem dividir a funcáo por z, existem situacóes onde esta 


divisáo irá simplificar bastante os cálculos. 
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É. TRANSFORMADA Z INVERSA 


EXPANSÃO EM FRAÇÕES PARCIAIS — POLOS REAIS SIMPLES 
e Exercício: Encontre a transformada Z inversa de 


Z 


Р) = T04DG £026 + 03) 


e Solucáo: 


01K. E PN. ak 
ao = (50 0.1) кы + 50(—0.3)*, e 


88 
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RANSFORMADA 4 INVERSA 
= u EXPANSÃO EM FRAÇÕES PARCIAIS — POLOS REAIS E 
COMPLEXOS CONJUGADOS 


e Para os polos complexos conjugados, utilizam-se as seguintes 


relações: 


e “sin(wa)z 
Z(e “*sin(k = — — OOO 
{е ѕіт(коа)) z^ — 2е-&с05(04)2 + e * 


z|z — e “cos(wa)] 
{е “*cos(k = ZII 
ja z^ — 2e *cos(og)z + ee 
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RANSFORMADA NVERSA 
É. EXPANSÃO EM FRAÇÕES PARCIAIS — POLOS REAIS E 
E COMPLEXOS CONJUGADOS 


e Exemplo 10: Obter a transformada Z inversa de 


z?-2z41 


dc (z — 0.1)(z2 + z + 0.5) 


e Solucáo: 
e Fazendo-se a expansão em frações parciais de F(z)/z: 


F(z) _ 23 + 22+ 1 A, А, " Az +В 
z z(z—041)(z22+z+05) z z-01 z2+z+0.5 
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RANSFORMADA INVERSA 
Ë EXPANSÃO EM FRAÇÕES PARCIAIS — POLOS REAIS E 
COMPLEXOS CONJUGADOS 


e Exemplo 10: continuação 
e Os resíduos dos polos reais podem ser calculados como mostrado 


para o caso anterior: 


F 


e Os coeficientes 4 e B podem ser encontrados igualando-se os 
coeficientes da função F(z)/z original com os da função na forma 


expandida. 
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RANSFORMADA NVERSA 
É. EXPANSÃO EM FRAÇÕES PARCIAIS — POLOS REAIS E 
E COMPLEXOS CONJUGADOS 


* Exemplo 10: continuação 
* Desta forma, após algumas manipulações algébricas, tem-se: 
aa + А, +A=1 
z2:0.94, + 4, — 01A+ B = 0 
2 0.44, + 0.54, — 0.1B = 2 
29: —0.054, = 1 
e De onde se obtém os valores de A e B: 


A = 1.311 = —1.557 
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RANSFORMADA INVERSA 
É. EXPANSÃO EM FRAÇÕES PARCIAIS — POLOS REAIS E 
COMPLEXOS CONJUGADOS 


e Exemplo 10: continuação 


e Multiplicando-se a expansão obtida por z, tem-se: 


19.6897  1.311z^ — 1.557z 


F(z) = —20-4 ———— 4 A A 
(2) z- 0.1 22 +2 + 0.5 


e Os dois primeiros termos podem ser localizados facilmente na 


tabela de transformadas. 


e O terceiro termo deve ser “arrumado” para que a tabela possa ser 


utilizada. 
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RANSFORMADA INVERSA 
Ë EXPANSÃO EM FRAÇÕES PARCIAIS — POLOS REAIS E 
COMPLEXOS CONJUGADOS 


e Exemplo 10: continuação 
e Pode-se reescrever a fração envolvendo os pólos complexos 


conjugados como: 


1.3112? — 1.5577 _ 1311[z — e “cos(w4)] — Cze *sin(oq) 
z2+2(-0.5)z+ 05 . z^ — 2e “cos(wg)z+te ** 


e Os termos do denominador podem ser facilmente obtidos: 


e * = v0.5 = 0.707 cos(wa) = о ui = — 0.707 


e Com: wq = ?*/, e sin(wa) = 0.707 
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RANSFORMADA INVERSA 
É. EXPANSÃO EM FRAÇÕES PARCIAIS — POLOS REAIS E 
COMPLEXOS CONJUGADOS 


e Exemplo 10: continuação 


e Igualando-se os coeficientes de z no numerador, obtém-se: 
—1.311e “cos(wa) — Ce “sin(wa) = —0.5(C — 1.311) = — 1.557 


e Resultando em: С = 4.426. 


e Utilizando-se a tabela de transformadas Z, pode-se obter f(k), 
para k > 0: 


f(k) = —208(k) + 19.689(0.1)* + (0.707) [1.311 cos(37k/,) z 4.426sin(37k/,)| 
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É. TRANSFORMADA Z INVERSA 


EXPANSÃO EM FRAÇÕES PARCIAIS — POLOS REPETIDOS 


e Para o caso de polos de multiplicidade r: 


T n 
N(z Ај; А; 
(Z — Z1) ID=r+17 — Z; 2. (Z — Z1) к 


e Os coeficientes para os polos repetidos podem ser obtidos utilizando-se 


a expressão: 


1 di | 
Аі; — (iD dzi 1 (z = 21) F(z) : l = 1,2, тее g 


2—21 
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É. TRANSFORMADA Z INVERSA 
sil EXPANSÃO EM FRAÇÕES PARCIAIS — POLOS REPETIDOS 
1 


e Exemplo 11: Obter a transformada Z inversa de F(z) = —Ñ—  ——o—a 
z^(z — 0.5) 


e Solucáo: 
e Fazendo-se a expansáo em fracóes parciais de F(z)/z: 


AO CN VECES A 
Z Z3(z—0.) 23 z? z z—05 


e O resíduo do polo em 0,5 vale: 


z=0.5 
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É. TRANSFORMADA Z INVERSA 


EXPANSÃO EM FRAÇÕES PARCIAIS — POLOS REPETIDOS 


e Exemplo 11: continuação 
e Já o resíduo dos polos repetidos podem ser obtidos como 


mostrado a seguir: 


F(z) 1 
a ES - "E 
Mies 2 — 0.5 É 
=0 
1d „F(z d 1 —1 
A. = ——z? ) = —— = — — —4 
1!dz Z dz z — 0.5 = (z — 0.5)? M 
2 = _1\(— 
„сл PO) dd 31 | ICDD) g 
2! dz 2 loe 2 dz (z — 0.5) ЛИР 2 (z — 0.5) 6 
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É. TRANSFORMADA Z INVERSA 


EXPANSÃO EM FRAÇÕES PARCIAIS — POLOS REPETIDOS 


e Exemplo 11: continuação 


e Multiplicando-se F(z)/z por z, pode-se escrever que: 


1 OZ ЕЕ e 
F(z) = 27 ^—42 – 8 


(2-0.5) z-05 


e Utilizando-se a tabela de transformadas, obtém-se: 


8(0.5)F — 28(k — 2) — 46(k — 1) — 86(k), k 20 


pao =| 0 k «0 
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É. TRANSFORMADA Z INVERSA 


EXPANSÃO EM FRAÇÕES PARCIAIS — POLOS REPETIDOS 


e Exemplo 11: continuação f(0)=8-8=0 


e Calculando-se f(k) рага К = 0,1e 2, tem-se: f(1) = 8(0.5) – 4 = 0 
f(2)=8(0.5)2-2=0 


e Assim, a transformada inversa também pode ser expressa por: 


“(05 3, k23 
pao =| 0 ‚ k«3 


e Esta solucáo poderia ter sido obtida diretamente, utilizando-se a 
propriedade de atraso no tempo, se ao invés de expandir F(z) em 


frações parciais, a função fosse colocada na forma: 


F(z) = z 
— ре 
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ЭЖ TEOREMA DO VALOR FINAL 


e O valor inicial de um sinal causal (sinal nulo para instantes de tempo 


negativos) é dado por: 


u(0) = lim U(z) 


e Da definicáo de sinais causais: 


e Prova: 


00 


U(z) = > u(k)z k = u(0) +u(1)z l+ u(2)z ? + -- 


k=—co 


e Aplicando-se o limite: 
lim U(z) = u(0) 
Z — 00 
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m _ TEOREMA DO VALOR FINAL 


e Se uma sequência possui um limite constante quando k tende a infinito, 
entáo: 


f (со) = tim f(k) = ul re 


e Limitações (casos onde limite não existe): 
e Sequência ilimitada (unbounded). 


e Sequência oscilatória 


я Deve ser utilizada apenas рага sistemas estáveis! 
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E: TEOREMA DO VALOR FINAL 
Е EXEMPLO 1 


e Verifique a validade do teorema do valor final para ита sequéncia 


exponencial decrescente: 
Z 7 
aan ES 
z— ear 
e Solucáo: 


e Supondo a > 0 e utilizando a função no domínio do tempo, tem-se: 
fs = lim poU =) 


e Utilizando-se o teorema do valor final: 


f (e) = lim (=) (+) = 
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E: TEOREMA DO VALOR FINAL 
Е EXEMPLO 2 


e Obter o valor final da sequéncia cuja transformada Z é dada por: 


z^(z —a) 


FU --0G-BG-o) 


e Solucáo: 


e Aplicando-se o teorema do valor final, tem-se: 
m z—1 z? (z — a) Е 1-а 
iE: Cz) fe -I(E — b)G — 5 “a-Da-o 


o Eng. de Controle e Automação – Rejane de Barros Araüjo/André Ferreira 


E TEOREMA DO VALOR FINAL 
Е EXEMPLO 2 


e Solucáo: 
e Observação: Expandindo-se F(z) em frações parciais, tem-se três 
termos: 
e Degrau unitário. 
e Exponencial b“ 


e Exponencial c“ 


e A validade do resultado depende da sequéncia convergir para uma 


sequência constante quando k tende a infinito: 


Ib| < 1 Ic| < 1 
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- COMANDOS ÚTEIS DO MATLAB 


e Suponha a seguinte funcáo de transferéncia no domínio Z: 


5(z + 3) 


EQ. 
2 23 + 0,172 + 0,47 


>> пит = 5*[1 3] 
>> den = [1 0.1 0.4 0] 
>> roots(den) 
апѕ = 
0 
-0.0500 + 0.6305 
- 0.0500 - 0.63051 
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- COMANDOS ÜTEIS DO MATLAB 


>> [r, p, k] = residue(num,den) 


r= 

- 18.7500 - 2.47831 
- 18.7500 + 2.4783i 
37.5000 


- 0.0500 + 0.6305i 
- 0.0500 - 0.6305i 


[] 
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- COMANDOS ÜTEIS DO MATLAB 


>> g=tf(4,[1 2 4]) 
Transfer function: 


$^2+25+4 


>> gd=c2d(g,0.1,'zoh') 
G(s) = 


Transfer function: s2+25+4 


0.01867 z + 0.01746 


2^2 - 1.783 z + 0.8187 
Sampling time: 0.1 


>> 
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- COMANDOS ÜTEIS DO MATLAB 
e Outras funções: 


e Transformação do domínio Z para o domínio s: 
e d2c 


e Transformada Z e Transformada Z inversa (utilizando a toolbox de 
manipulação simbólica): 
e Ztrans 


e |ztrans 


e syms kn wz 
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É. SOLUÇÃO DE EQUAÇÕES A DIFERENÇAS USANDO 
di TRANSFORMADA Z 


e Resolver a equacáo a diferencas linear abaixo: 
x(k + 2) — (9/;)x(k + 1) + (1/5)x(K) = 10) 


com as seguintes condições iniciais x(0) = 1, x(1) = 5/2 


e Solucáo: 
e Aplica-se a transformada Z em ambos os lados da equacáo: 
[z2X (2) — z?x(0) — zx(1)] — (?/5)IzX C2) — zx(0)] + (1/5)x 2 = 7/, 4 
e Isolando-se a variável de interesse X(z): 


z|1 + (z + 1)(z — 1)] Ze 


XD GE TDE=DE=05 021020205) 
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SOLUÇÃO DE EQUAÇÕES A DIFERENÇAS USANDO 
TRANSFORMADA Z 


HITTUTOFID 
mu 


e Aplica-se a expansão em frações parciais de X(z)/z: 


X(2) _ z? A11 A12 + 42 
Z (2 1)2(2— 0.5) - (z - 1)? (z—1) (2-0.5) 


e Calcula-se os resíduos: 


EM ze 1 
Au SUPE CREU 01-087. 
"n d 20 22 Ze _ 2 1 _ 
= ile- ) zu - As nu ,, 05 (052 
o Е ne Е (0.5)? Е 
“6057 | ED. 05-125 
zs 2=0.5 


e Eng. de Controle e Automação – Rejane de Barros Araújo/André Ferreira 


SOLUÇÃO DE EQUAÇÕES A DIFERENÇAS USANDO 
ni TRANSFORMADA Z 


e O que resulta em: 


MP Dt eda 


e Aplicando a transformada Z inversa, obtém-se a solucáo da 


equacáo de diferenca como: 


x(k) = 2k + (0.5)* 
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7 ^N 
s. RESPOSTA TEMPORAL DE SISTEMAS DISCRETOS 


e А resposta temporal de um sistema linear discreto é a solução da 


equação a diferenças que “governa” o sistema. 


e Para sistemas LTI, a resposta devido às condições iniciais e devido à 
entrada podem ser obtidas separadamente. Estas respostas somadas 


fornecem a resposta total do sistema. 


e A resposta devido a entrada (resposta forcada) é o somatório da 
convolucáo (convolucáo discreta) para esta entrada com a resposta ao 


impulso unitário. 
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7 ^N 
s. RESPOSTA TEMPORAL DE SISTEMAS DISCRETOS 


e A resposta de um sistema discreto ao impulso unitário é chamada de 


sequéncia de resposta ao impulso. 


e Esta sequéncia pode ser utilizada para representar a resposta do 


sistema LTI para uma sequéncia de entrada arbitrária: 


(u(k)) = (u(0),u(1), ...,u(k), ...) 


Representando a sequéncia de entrada como um conjunto de 
impulsos discretos 


00 


u(k) = > TO — D 


1=0 
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^ ^ 
s „RESPOSTA TEMPORAL DE SISTEMAS DISCRETOS 


e Para um sistema linear, o princípio da superposição é válido e a saída do 
sistema para uma entrada pode ser representada como um somatório de 


sequéncias de resposta ao impulso: 
UC) = th(t)ju(0) + {hlt — 1)}и(1) + (h(t— 2)juQ) + eee A — Dyu (i) ++ 
e Assim, a saída em um instante k vale: 
y(k) = h(k) *u(k) = У h(k — i)u(i) 
En 


convolução 
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(a N 
uu. RESPOSTA TEMPORAL DE SISTEMAS DISCRETOS 


e Para um sistema causal, a resposta ao impulso no tempo i é uma 
sequência que começa no instante i, e cuja resposta atrasada h(k — i) 


satisfaz a relação: 
h(k — i) = 0, {>К 


e Resposta de um sistema LTI discreto a um impulso em iT: 


e(k — i) | fh(k— i) | 
iT iT 
LTI System 
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Eu ^N 
s RESPOSTA TEMPORAL DE SISTEMAS DISCRETOS 


e Assim, um sistema causal é aquele cuja resposta ao impulso é uma 


sequéncia causal. 


e Desta forma, a equação anterior y(k) = h(k) * u(k) = >. h(k — iju(i) 
reduz-se a: n 


y(k) = u(0)A(k) + u(1)A(k — 1) + u(2)h(k — 2) +++ u(k)h(0) 


k 


y() = > иһ — D 


1=0 
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7 ^N 
s. RESPOSTA TEMPORAL DE SISTEMAS DISCRETOS 


e Fazendo uma mudanca na variável do somatório para j = k — i, tem-se: 


y(k) = u(0)A(K) + u(1)A(k — 1) + u(2)h(k — 2) + -- + u(k)A(0) 
k 


y() = > u()h(k — D 


1=0 


| 


y(k) = u(k)h(0) + u(k — 1)A(1) + u(k — 2)h(2) + + u(0)h(k) 
k 
y() = > u(k Da) 


j=0 
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RESPOSTA TEMPORAL DE SISTEMAS DISCRETOS. 
ой ТЕОВЕМА 


e À resposta de um sistema LTI discreto para uma sequência de entrada 
arbitrária é dada pelo somatório da convolução da sequência de entrada 


e da sequência de resposta ao impulso do sistema. 


e Exemplo: k = 2 Ж Иен = > мк TO 
i=0 


| 


2 
yQ) = » uD - D 
і=0 


y(2) = и(0)А(2) + u(1)h(1) + u(2)A(0) 
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e ^N 
nam | RESPOSTA TEMPORAL DE SISTEMAS DISCRETOS 
e Exemplo: continuação 
| Entrada | Saa | cor _ 
Y u(0) h(2) 
DES 
o Eng. de Controle e mon и Rejane de Barros Araújo/And?s Ferreira A 


І RESPOSTA TEMPORAL DE SISTEMAS DISCRETOS ` 
sili TEOREMA DA CONVOLUCAO 


e A transformada Z da convolucáo de duas sequéncias no tempo é igual 
ao produto de suas transformadas Z. 


e Prova: 
e aplicando a transformada Z em: 


y(k) = МЕ) «u(k) = > hlk дид) 
ї=0 


үх) = y(k)z = Ү(х) = | u(i)h(k — 5 q 
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(a N 
« RESPOSTA TEMPORAL DE SISTEMAS DISCRETOS 


e Mudando a ordem dos somatórios e fazendo j = k — i: 
Y(z) = > > u(Dh(gj)z 8+ 
і=0 j--i 


e Usando a propriedade de causalidade, tem-se: 


Y(z) = n «| У e 
i=0 j=0 
y 


Y(z) = H(z)U(z) 
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С ^N 
- FUNCAO DE TRANSFERÉNCIA PULSADA 


e Afuncáo H(z) é chamada de função de transferência pulsada. 


e Quando se utiliza o teorema da convolucáo na resposta de um sistema 
LTI, é possível usar a transformada Z para obter a saída do sistema: 


e 1) Aplicando-se a transformada Z na sua entrada; 


e 2) Multiplicando-se a transformada Z da entrada pela função de 


transferência; 


e 3) Calculando-se a transformada Z inversa, para obter a 
sequência de saída temporal. 
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С ^N 
E _ FUNÇÃO DE TRANSFERÉNCIA PULSADA 


e Exemplo: Achar a resposta ao impulso h(k) do seguinte sistema 


discreto, supondo y(0) = 0: 


y(k + 1) — ay(k) = u(k) 
e Solucáo: 
e a) A partir da equacáo a diferencas 
e Fazendo u(k) = ó(k), tem-se: 
y(1)21 
у(2) = ау() = а 

y(3) = ау(2) = а? 

acolat . 151283, 
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С ^N 
- FUNCAO DE TRANSFERÉNCIA PULSADA 


e Exemplo: continuacáo 
e Solução: 
e a) Usando a transformada Z: 


Y(z 1 
U(z) z-a 


y(k+1)—ay(k)=u(k) —>  H(z)- 


e Aplicando-se o teorema do atraso: 


dat , К=1,23,... 
һо) = | 0 , k«1 
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(< К 
- EXERCÍCIOS 


e Dada a equação a diferenças abaixo, determine sua função de 


transferência pulsada e sua resposta a um degrau unitário discreto: 


y(k +1) — y(k) = u(k + 1), y(0) = 
e Solução: 
YZ) 7 
Hi U(z) z-1 


ro = (EEA =: 


Fi t1, к= 0,1,2,3,... 


y(k) = ; k < 0 
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RESPOSTA TEMPORAL DE SISTEMAS DISCRETOS ` 
sil CONCLUSÃO 


e À resposta de um sistema discreto pode ser obtida com a 
utilização da transformada Z. Esta abordagem é mais rápida e 
simples do que a utilização direta da equação a diferenças, 


especialmente para sistemas de ordem elevada. 
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f E 
- TRANSFORMADA Z MODIFICADA 


e Aamostragem e a transformada Z capturam somente os valores de uma 


funcáo contínua em determinados pontos de amostragem. 


e Para obter os valores da funcáo contínua entre intervalos de 
amostragem, é necessário atrasar a forma de onda a ser amostrada de 
uma fracáo do intervalo de amostragem antes que a amostragem 


ocorra. 


e Mudando-se esse atraso de tempo, podem-se variar os pontos onde o 


sinal é amostrado. 


e Este procedimento recebe o nome de transformada Z modificada. 
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f ^ 
- TRANSFORMADA Z MODIFICADA 


e Considere um sinal causal contínuo y(t), que é amostrado a cada T 


segundos. 


e Inserindo-se um atraso Т; < T antes do amostrador, tem-se: 


| = T SET 
жд AU X (KT) 
» Delay T; — p 


e A saída do elemento que causa o atraso é o sinal contínuo: 


_fyt-—Ta) , t20 
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м 


- EFEITO DO ATRASO 


e Sequência causal: 


e Amplitude é sempre nula para tempo negativo. 


e O atraso sempre resultará em um valor nulo para t=0, 


independente da amplitude do sinal. 


e Deve-se dividir o atraso em um avanço de uma determinada 
fração do período de amostragem, seguido do atraso de um 
período de amostragem completo, para se evitar valores 


inapropriados para t = 0. 
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e Exemplo: 
ө T, = 0,3 S 
eT = 1,0 S 
em = 0,7 S 


DECOMPOSIÇÃO DO ATRASO 


Та =Т= тт О<т<1 


| 


Т, 
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Е. _ DECOMPOSICAO DO ATRASO 
° y_ (t) = y(t) atrasado de um período de amostragem, de tal forma que: 


Yat) = y(t - T +mT) = y_-1(t + mT) 


e Amostrando a forma de onda atrasada com um período de amostragem 


T, tem-se: 


ya (KT) = y_¿(kT + mT), k = 0,12, ... 
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ы 
- TRANSFORMADA Z DO ATRASO 


ya) =2"Y(2) 


| 


Y(z, m) = Z my UT); 


Y (z,m) = z 1Z(y(kT + тТ)} Eq. (1) 
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7 E TRANSFORMADA Z DO ATRASO ^N 
"а= EXEMPLO 1 


e Degrau unitário discreto: 
e Como essa função apresenta um valor constante para todos os 
intervalos de tempo, deslocando ou atrasando o instante de 
amostragem não afeta os valores obtidos. 


e Neste caso, a transformada Z modificada será sempre igual à 
transformada Z, sendo dada por: 


1 


Ym) = z 
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7 E: TRANSFORMADA Z DO ATRASO ^N 
e EXEMPLO 2 


e Exponencial: y(t) = Ae ?* 
e Um avanco no instante de amostragem mT multiplica o sinal 


amostrado por e PT. 


y(kT + mT) = Ае "Т = Ag-pmTo-pkT k = 0,52, ... 


= —pmT 
Z{y(kT + тТ)} = Ae Урт 


Ae ?mT 


Y (z, m) = 7 — o pT 


o- pT Utilizando a Eq. (1) 
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F: TRANSFORMADA Z DO ATRASO 
L EXEMPLO 2 


e Exponencial: continuacáo 
e Para p = 4,Т = 0,2 s, A = 1 e T = 0,7, tem-se: 
e m = 0,3 
«kg PMI = 0787 


e g PT = g-08 = 0,449 
e Atransformada Z modificada: 


Ae РТТ 


Y(z,m) = -oT ——>  Y(zm)- 


0,787 
2 — 0,449 
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С E: TRANSFORMADA Z DO ATRASO 
mo EXEMPLO 2 


e" (k +m—1) Т 


0 


e"? (k+m)T 


0 T 2T 3T КТ 
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¡GA 


UN 
RESPOSTA EM FREQUÉNCIA DE SISTEMAS DISCRETOS 


мттитонр 
mu 


Resposta de um sistema discreto a uma senóide amostrada. 
e Senóide de saída com a mesma frequência da senóide de 
entrada, com  possibilidade de deslocamento de fase e 


modificacáo da magnitude. 


O deslocamento de fase e modificacáo (escalonamento) da magnitude 
definem uma funcáo complexa da frequéncia: 


e Resposta em frequência do sistema. 
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(a N 
a... RESPOSTA EM FREQUÊNCIA DE SISTEMAS DISCRETOS 


мттитонр 
mu 


e Lembrando que, para sistemas contínuos: 
Н(јо) = H(s)| _ 
S= Jo) 


e Considere o seguinte sinal discreto de entrada obtido como uma 
representação amostrada por impulso: 


u'(t) = У u(kT)ó(t — КТ) 


e onde u(kT) é o valor de u no instante kT e ó(t — kT) representa um 
delta de Dirac no instante KT. 
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(a N 
a... RESPOSTA EM FREQUÊNCIA DE SISTEMAS DISCRETOS 


мттитонр 
mu 


e Aplicando-se a transformada de Laplace no sinal discreto, tem-se: 


00 


U*(s) = 2 u(kT)e ETIS 


k=0 
e Pode-se definir uma funcáo de transferéncia para entradas amostradas 
como: 


Ү (5) 
U*(s) 


H*(s) = 


condições iniciais nulas 
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e 7x 
- RESPOSTA EM FREQUÉNCIA DE SISTEMAS DISCRETOS 


e Desta forma: 
H*(jo) = Н"(5)| _ 
S= Jo) 


e Lembrando que z = es”: 


H(z) = H* (s) Aia 


e Assim, a resposta em frequência é dada por: 
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(a N 
a... RESPOSTA EM FREQUÊNCIA DE SISTEMAS DISCRETOS 


мттитонр 
тми 


e Uma outra forma de se conseguir a resposta em frequéncia, sem 
recorrer a amostragem por impulso, pode ser obtida considerando a 


exponencial complexa amostrada: 
= ikwoT 
u(kT) = uge7*o 


u(kT) = ug[cos(k«gT ) + jsen(kogT )], k == 0:12... 


e Atransformada Z desta sequéncia de entrada vale: 


UG) = uo Jer 
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(a N 
a... RESPOSTA EM FREQUÊNCIA DE SISTEMAS DISCRETOS 


мттитонр 
mu 


e Supondo-se que a funcáo de transferéncia possa ser colocada na 


forma: 
N(z) 


H(z) = === 
e onde N(z) é um polinômio de no máximo ordem n e que todos os polos 


p; estão dentro do círculo unitário. 


e Desta forma, a saída do sistema sujeito à entrada exponencial vale: 
Y(z) №2) 7 

z) = || Un —— 

¡=1(Z — Pi) 02 — ejooT 
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(a N 
a... RESPOSTA EM FREQUÊNCIA DE SISTEMAS DISCRETOS 


мттитонр 
mu 


e Expandindo-se em frações parciais: 


Az - Biz 
Y(z) = A —— 


— ejwoT — 1). 
z egg z — p; 


e Aplicando-se a transformada Z inversa: 
n 
y(kT) = AeJ*oT + 2 B;p;", k = 0,1,2,... 
і=1 
e Lembrando аа suposicáo que os polos estáo dentro do círculo unitário: 


Yss (КТ) = AelkwoT 


@) * рага k suficientemente grande. 
SS 
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(a N 
a... RESPOSTA EM FREQUÊNCIA DE SISTEMAS DISCRETOS 


NiTTUTOFID 
mu 


e O termo A da expansão em frações parciais vale: 


А = 2 — еЇ®Т) 


z=ejwoT 


A = Н(еі®Туи, 


e Assim, pode-se escrever a saída em regime permanente na forma: 


Yes (KT) = |Н(е1®еТу|ц„еЛк®оТ+ ue] para k grande. 


| 


yss (KT) = |H(e/“) luycos[kwyT + ZH(ej%oT)]| 
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(a N 
a... RESPOSTA EM FREQUÊNCIA DE SISTEMAS DISCRETOS 


мттитонр 
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e Pode-se substituir o cosseno na expressáo anterior por um seno, 
podendo-se entáo concluir que a resposta a uma senóide amostrada é 
uma senóide de mesma frequéncia, escalonada em magnitude e com 


um deslocamento de fase devido à: 


|H (eJ 9 oT) | ZH (eJ oT) 


e Este é o mesmo resultado que é obtido utilizando-se o sinal de entrada 


amostrado por impulso. 
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(a N 
a... RESPOSTA EM FREQUÊNCIA DE SISTEMAS DISCRETOS 


мттитонр 
mu 


e Exemplo: Obtenha a resposta em regime permanente do sistema 
discreto: 
1 
H(z) =  — 
(2) = (бту — 05) 
e devido à uma entrada senoidal u(kT) = 3 * cBls(0,2k). 


e Solucáo: 
ys (KT) = |H(e1%2)| + 3 + cos(0,2k + ZH (eJ9?)) 


1 


1 
Yss (KT) = (2702 0,1) (e92 — 0,5) * З ж cos c + LIDOS) 


(ej0,2 — 0,1)(eJ02 — 0,5) 
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e Ys (KT) = 6,4 * cBs(0,2k — 0,614) 


RESPOSTA EM FREQUÉNCIA DE SISTEMAS DISCRETOS ` 
silii PROPRIEDADES 


e 1) Ganho DC: 
e O ganho DC é igual a H(1). 


e Prova: 


H* (jo) = H'(Z)l, jer = Н(е?°Т, 


| 


H(ejo7, M - H(z) |. - H(1) 
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RESPOSTA EM FREQUÉNCIA DE SISTEMAS DISCRETOS ` 
silii PROPRIEDADES 


e 2) Natureza Periódica: 
e À resposta em frequência de um sistema discreto é uma função 


periódica da frequência, com período igual à w, = 27/т rad/s. 
e Prova: 


eJoT = cEs(wT) + jsen(wT) 


e É uma função periódica com período о, = 27/;. rad/s. Assim H(e/%7) 


também será uma função periódica com esse mesmo período. 
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RESPOSTA EM FREQUÉNCIA DE SISTEMAS DISCRETOS. 
m PROPRIEDADES 
° 3) Simetria: 
e Para funções de transferência com coeficientes reais, a magnitude da 


função de transferência é uma função par da frequência, enquanto 


sua fase é uma função ímpar desta mesma frequência. 


e Prova: 


e Para uma frequência negativa: H(e7107) = H(eloT) 


e Para coeficientes reais: H(e/47) = H(e/07) 
A fien) = CeT) 


ZH(e-J9T) =-/H(e1%7) 


e Desta forma: H(e 7*7) = H(ejeT) 
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RESPOSTA EM FREQUÉNCIA DE SISTEMAS DISCRETOS ` 


i PROPRIEDADES 


e Resposta em frequência de um sistema discreto de 2º ordem: 


3 


2.5 


0 200 400 600 800 1000 1200 
w rad/s 
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RESPOSTA EM FREQUÉNCIA DE SISTEMAS DISCRETOS ` 
silii PROPRIEDADES 


e Observações: 


e 1) Apenas é necessário obter o valor de H(e/%7) entre o ganho DC e 


25/5 
e 2) As frequências negativas podem ser obtidas por simetria. 
e 3) A resposta em frequência se repete рага о > “s/, 
e 4) Para os casos onde a resposta em frequência é desprezível para 


w >ºs/,, não há superposição (overlap) de ciclos de resposta em 
frequência repetidos. 
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RESPOSTA EM FREQUÊNCIA DE SISTEMAS DISCRETOS ` 
PROPRIEDADES 


e Observações: 
e 5) Amostragem sem overlap: repetição periódica da resposta em 
frequência de um sistema contínuo. 


6) As respostas em frequência de sistemas físicos não possuem 
limitação de banda: Superposição (overlapping) de ciclos de resposta em 
frequência repetidos > folding. 


e 7) “s/, = frequência de folding. 


8) O folding resulta em uma distorção na resposta em frequência que 
deve ser evitada ou minimizada pela escolha adequada da frequência de 
amostragem ou através de filtragem. 
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qo \ 
E _ TEOREMA DA AMOSTRAGEM DE SHANNON 


e A amostragem deve refletir a natureza do sinal analógico. 

e Um sinal contínuo que foi amostrado pode ser recuperado a partir de 
uma sequência de amostras, desde que tenha sido convenientemente 
amostrado. 

e Quanto mais rápido o sinal for amostrado, garante-se a reconstrucáo 
adequada desse sinal. 

e O teorema da amostragem fornece um limite inferior na taxa de 
amostragem necessária para um sinal com banda limitada (sinal com 


uma largura de banda conhecida). 


IT 


"IT 
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q UN 
- TEOREMA DA AMOSTRAGEM 


e Um sinal com banda limitada com: 


fŒ 9 F(jo), F(jo) + 0, —oq € о € Om 


F(jo) = 0, para outras frequéncias 


onde F representa a transformada de Fourier do sinal, pode ser 
reconstruído a partir da seguinte sequéncia discreta: 


г = > 7080-7) 


k=— co 
se, e somente se, a frequência de amostragem w, = 2лт/Т satisfaz a 
condição о, > 20. 
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ы 
- TEOREMA DA AMOSTRAGEM 


e O espectro da forma de onda contínua pode ser recuperado usando um 


FPB ideal de largura de banda dentro da faixa wm < wp < wç/2. 


©; (Us 
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q UN 
- TEOREMA DA AMOSTRAGEM 


e Para poder recuperar o sinal original, sem distorcáo, a frequéncia de 


amostragem 


fr = = (Sendo Т =período de amostragem) 


tem que ser pelo menos duas vezes maior que todas as frequéncias 
presente no espectro do sinal original. Do contrário, ocorrerá o fenómeno 


de aliasing. 
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ГАВ 
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ES 


"a TEOREMA DA AMOSTRAGEM 


e Em engenharia de controle, recomenda-se um amplo coeficiente de 
segurança ao escolher fr. 
e Geralmente de 5 a 10 vezes maior que a maior frequência presente 


no sinal original. 


e Exemplo: Para um sinal com largura de banda de 10 rad/s, escolha 
uma frequência de amostragem adequada e seu correspondente 


período de amostragem. 


wp = 10 rad/s 


ws=60rad/s — T= be = 0.1s 


= ESCOLHA DA TAXA DE AMOSTRAGEM 


мттитонр 
ми 


e Limite inferior especificado pelo Teorema da Amostragem. 
e Regra prática: о, = kwm, 5 <k <10 


e Limitações: 
e À frequência de amostragem possui um limite superior que depende 
do atraso no sensor. 
e Tempo computacional deve ser levado em conta para atualizar o sinal 
de controle. 


e Lembrar sempre: a taxa de amostragem deve ser suficientemente 
rápida para fornecer uma boa representação das variáveis físicas 
(analógicas). 
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E MODELAGEM DE SISTEMAS DIGITAIS 


e Configuração Típica de um Sistema Digital 


Computador ou 
microprocessador 


Entrada externa Saída analógica 


Subsistema 
analógico 
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E... MODELO DE UM CONVERSOR A/D IDEAL 


мттитонр 


e Suposições: 
e As saídas do conversor A/D são exatamente iguais em magnitude as 
suas entradas (erros de quantização desprezíveis); 
e A saída digital do conversor A/D é obtida instantaneamente; 


e À amostragem é perfeitamente uniforme. 


e Conversor A/D ideal = amostrador ideal 
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la 
f MODELO DE UM CONVERSOR A/D IDEAL 


e As suposições anteriores são idealizações que são possíveis porque: 
e Erros de quantização são quase sempre pequenos, mas não são 


nulos; 


e Variações na taxa de amostragem ocorrem, mas também quase 


sempre são desprezíveis; 


e Conversores A/D apresentam um tempo de conversão finito. 


e De uma maneira geral, o modelo do amostrador ideal é aceitável para a 


maioria das aplicações em engenharia. 
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- MODELO DE UM CONVERSOR D/A 


e Suposições: 


e As saídas do conversor D/A são exatamente iguais em magnitude às suas 
entradas. 


e Asaída digital do conversor D/A é obtida instantaneamente. 


e As saídas do conversor D/A são constantes em cada período de amostragem. 
e Segurador de ordem zero (Zero-order-hold — ZOH). 
e Outros tipos: 
e Segurador de primeira ordem (First-order-hold — FOH): reconstrói o sinal 
usando retas. 


e Segurador de segunda ordem (Second-order-hold — SOH): reconstrói o 
sinal usando parábolas. 
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- MODELO DE UM CONVERSOR D/A 


e A relacáo entrada-saída de um conversor D/A ZOH pode ser expressa 


por: 


(u(k))2u(),  KkT<t<(k+DT, к= 0,1,2... 


onde u(k) é a sequéncia de entrada. 


kT Zero-Order KT (k+1)7 
Hold 


Positive Step 


Negative Step 
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E... FUNÇÃO DE TRANSFERÊNCIA DO ZOH 


e Usando a transformada de Laplace de um degrau unitário e o teorema 


do atraso para a transformada de Laplace: 


1 
L{1(t)} = m 
per 
L{1(1 -T)} = E KT (EDT 


e Entáo a funcáo de transferéncia do ZOH pode ser calculada como: 


1—e-^3T 
Gzou(s) = 
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a N 
a ANÁLISE DO ZOH NO DOMÍNIO DA FREQUÊNCIA 


мттитонр 


e Da expressão anterior: 


1—gJoT 


G l) = ——o D. 
zou (jo) jo 


e Reescrevendo a resposta em frequência na forma: 


jo 


‚ T "EU . T 
EA [e — duc) 
W J 
Hm. Т 
е ?°2 Т =, TSsen[o > 
= " [25en (03) = Te jur sen [o z) 
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e 


Gzou jo) = 


a UN 
a ANÁLISE DO ZOH NO DOMÍNIO DA FREQUÊNCIA 


мттитонр 


e Lembrando que sinc(x) = UR x 


wT T 
IGzog (jo )|ZGzogu (jo) = T |sinc (5) gene 


2 


e Conclusão: 
e Ângulo: 
e Diminui linearmente com a frequência. 
e Magnitude: 
e Proporcional à função sinc. 
e Oscilatória com valor de pico igual ao período de amostragem, 


ocorrendo na frequência zero. 
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f ^ 
- ANÁLISE DO ZOH NO DOMÍNIO DA FREQUÉNCIA 


e Magnitude da resposta em frequência do ZOH para T = 1 s. 


0 | 
-30 —20 -10 0 10 20 30 
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f E: EFEITO DO AMOSTRADOR EM UMA FUNÇÃO DE | 
"= TRANSFERÉNCIA EM CASCATA 


e Sistemas contínuos em cascata: 


U(s) X(s) Y (5) 


Y(s) = H2(s)X (s) = Ho (s)Hi(s)U(s) 


e Cada sub-sistema tem um amostrador: 
e Cada bloco tem saída e entrada amostrada, assim como uma função 


de transferéncia no domínio Z. 


Y (2) = H(z)U(z) = Ну(2)Н (2) ...Н.(2)0(2) 
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f E: EFEITO DO AMOSTRADOR EM UMA FUNÇÃO DE ^ 
e TRANSFERÊNCIA EM CASCATA 


e Os sub-sistemas são primeiramente agrupados, e depois utiliza-se um 
único amostrador: 


e Componentes não podem ser separados depois da amostragem. 


Y (2) = H(z)U(z) = (11H) ...H,)(Z)U (z) 


UG) <ç 0) |a Y(s) 
— HS 
T 
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f E: EFEITO DO AMOSTRADOR EM UMA FUNÇÃO DE ^ 
= TRANSFERÊNCIA EM CASCATA 


e Exemplo: Ache a sequência impulsiva amostrada e a função de 
transferência pulsada para um sistema em cascata formado pelas 


funções 


1 
Н; (5) = — H> (s) = 


S+2 s+4 


e Para: 
e 1) Os sistemas estáo conectados diretamente. 


e 2) Os sistemas estão separados por um amostrador. 
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q E: EFEITO DO AMOSTRADOR ЕМ UMA FUNÇÃO DE ^ 
"= TRANSFERÉNCIA EM CASCATA 


e Exemplo: solucáo 1. 
e Na auséncia de amostradores entre os sistemas, a funcáo de 
transferéncia do sistema total vale: 


2 


Эе EEG 


e Aresposta impulsiva pode ser obtida como: 
h(t) =e c — £e 
e Com a resposta impulsiva amostrada sendo: 
h(kT) = e-2kT — g-4kT c0 12, 


e E a função de transferência pulsada valendo: 
7 (e — e Dz 


2 
© A gre Eee) 
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f E: EFEITO DO AMOSTRADOR EM UMA FUNÇÃO DE ^ 
= TRANSFERÊNCIA EM CASCATA 


e Exemplo: solução 2. 

e Se os sistemas analógicos são separados por um amostrador, então 
cada um deles tem uma função de transferência pulsada: 

H(z) = —— 

2 1(2) z—e- ?T 


1 
@+2)(s+4) —— 27 


Haz) = 7 


H(s) = Н; (s)H2(s) = 


e Afuncáo de transferéncia pulsada do conjunto vale: 


222 


HG) = ¿NE 
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f E: EFEITO DO AMOSTRADOR EM UMA FUNÇÃO DE ^ 
e TRANSFERÊNCIA EM CASCATA 


e Exemplo: solução 2. 
e Aplicando a transformada Z inversa: 


2 
h(kT) = AT [а а = кы ad] 


2 
h(kT) = 2-2T = g-4T FB 3 = proe k = 0,1,2, — 
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É E: REPRESENTAÇÃO DISCRETA DO SUB-SISTEMA: ^ 
"а= D/A + PROCESSO + A/D 


e Representação genérica de um controlador digital: 


Computador 


| | 
І І Ргосеѕѕо 


r(t) + e(t) e(kT) u(kT) u(t) y(t) 
° 
= І I 


Algoritmo 
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f E: REPRESENTAÇÃO DISCRETA DO SUB-SISTEMA: ^ 
"а= D/A + PROCESSO + A/D 


e Supondo que o relógio (que náo é mostrado por simplicidade) trabalha 
com um período de amostragem constante e que a referéncia esteja 
sendo gerada internamente pelo computador: 

Computador 


Processo 


e Mistura de variáveis nos domínios s e Z: 
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f E: REPRESENTAÇÃO DISCRETA DO SUB-SISTEMA: ^ 
es D/A + PROCESSO + A/D 


e Deve-se obter a função de transferência discreta do sub-sistema: 
Processo 


e O sistema (b) é o modelo discreto do sistema (a) que possui a entrada 
discreta u(kT) e a saída discreta y(kT). 
e No modelo discreto (b), não está acessível o y(t) presente no modelo 


Ө ° só estando disponível y(kT). 
sa 
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"ы ^N 
E CONVERSOR D/A DE ORDEM ZERO 


e O conversor D/A de ordem zero aproxima os valores amostrados por um 


polinómio de ordem zero: 


Polinômio de 
ordem zero 


D/A 


Ordem zero 


u(t) =u(kT), para kT <t < (k + 1)T 
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Q E: APLICAÇÃO DA ENTRADA IMPULSIVA EM: ^ 
т= D/A — PROCESSO — A/D 


e Fazendo u(kT) ser uma entrada impulsiva, a resposta u(t) do conversor 
D/A será: 


u(kT) 
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С E APLICAÇÃO DA ENTRADA IMPULSIVA EM: ^ 
ge D/A — PROCESSO - A/D 


e Como visto anteriormente, a resposta do conversor D/A é um pulso, que 
corresponde a uma combinacáo de degraus. 


e Chamando o degrau de d(t): 


O pulso é dado por: 


u(t) = d(t) – d(t — T) 
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E: APLICAÇÃO DA ENTRADA IMPULSIVA EM: B 
- D/A — PROCESSO - A/D 


e Para determinar a resposta de G(s) (processo) a essa entrada, é 


necessário aplicar a transformada de Laplace em u(t). 


e Lembrando que: 
LY (t а)) = е ®Е(5), sendo F(s) = £U (t)) 
e Pode-se escrever: 


Ltu(t)) = даб) — d(t — T)j 
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f E: APLICAÇÃO DA ENTRADA IMPULSIVA EM: 
z D/A — PROCESSO - A/D 


U(s) = £id(t)) — Lid(t — T)j 


1—e-^5T 
U(s) — 


o Eng. de Controle e Automação – Rejane de Barros Araújo/André Ferreira 


f E: APLICACÁO DA ENTRADA IMPULSIVA EM: 
um D/A — PROCESSO — A/D 


e Lembrando que a resposta Y (s) do processo devido à entrada dada por: 


Y(s) = G(s)U(s) 


1—e 53 
Y (s) = G(s) — — 


e Desta forma: 


— a T 
y(t) = 17 [== 19) 
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ES 


° 


E APLICAÇÃO DA ENTRADA IMPULSIVA EM: 
n D/A — PROCESSO - A/D 


e Esta resposta y(t) passa pelo conversor A/D e gera y(kT). 


e Desta forma: 
1 — e sf 
y(kT) = (f E 29) 


e Atransformada Z de y(kT) será: 


(1— e^t 
Y(z) = zty(kT)) = 3 (< 2 | E 29) 
YO e 2 2 =. fe Du 
5 Mer P 
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t=kT 


a 
a 


E APLICAÇÃO DA ENTRADA IMPULSIVA EM: 
т D/A — PROCESSO — A/D 
o Fazendo: 
fo = c (E 
o Tem-se: 


YE) = a0 (0 - fE- T))| = aU (T) - f(kT - T) 
Y(z) = zu (k)j — sU (k — 1)) 
YG) = zif} - z ^ z(f (00) 


Y(z) = (1-2 *)2{7(КТ)) 


Q............ — Rejane de Barros Araújo/André Ferr 


o Entáo: 


ES 


f E: APLICAÇÃO DA ENTRADA IMPULSIVA EM: ^ 
e D/A — PROCESSO - A/D 


«(^£ 
t-kT 

Y(z) = (1 — z `1)Z (с E) | 
t=kT 


* Deve ser interpretada como a transformada Z da sequência obtida pela 


o Assim: 


amostragem do sinal 
f(t) = f =, com t = kT 
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f E: APLICAÇÃO DA ENTRADA IMPULSIVA EM: ^ 
e D/A — PROCESSO - A/D 


o Como a resposta ao impulso Y (z) é igual à funcáo de transferéncia do 


N 


sistema H(z), tem-se: 


H(z) = (1 — z Dz (e 22 
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q E: APLICAÇÃO DA ENTRADA IMPULSIVA EM: 
L D/A — PROCESSO — A/D 


o Assim, para calcular H(z) deve-se seguir os seguintes passos: 


e 1) Achar a transformada inversa de G(s)/s: 
-1 fG 
r= (88) 


• 2) Fazer f(kT) = f(t)|z==r e encontrar a transformada 4: 


F(z) = sU (kT)) 


e 3) A funcáo de transferéncia discreta do sistema D/A — G(s) — A/D 


será: 
H(z) = (1-z"Dr(z) 
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f E: APLICAÇÃO DA ENTRADA IMPULSIVA EM: ^ 
@ 


т D/A — PROCESSO — A/D 
o Exemplo: Para o sistema abaixo, calcule a função de transferência 
pulsada Н(2): 
u(kT) u(t) a y(t) y(kT) 
D/A A/D = O 
e Solução 1: 
а 
G(s) = 
s+a го 102091) — licen 
GG) _ a S s(s + a) 
(S (s + a)s 
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С E: APLICACÁO DA ENTRADA IMPULSIVA EM: 
ш” D/A — PROCESSO — A/D 


° Solução 2: f (kT) = 1 — e-**T 


o Aplicando-se a transformada Z: 


= —aT 
FG) = ZU (TD) = zü e) = се О” 


+ 3) Assim: 


(1-e Dz 


Hio cbe gu = ре уле ma 


2 — =) (1 — e “Dz 


z J(z—1)(z — e) 


H(z) = | 
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é E: FUNÇÃO DE TRANSFERÉNCIA DISCRETA EM 
B= MALHA FECHADA 


C(z)G(z) 
1 + C(z)G(z) 


Ga(z) = 


o Polinómio característico: 


1-4 C(z)G(z) = 0 
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Е. ERRO DE REGIME PERMANENTE 


мттитонр 


2) — ganho de malha aberta 


o Entradas padrão: 
e Degrau amostrado 
e Rampa amostrada 
e Parábola amostrada 


o Pode-se calcular o erro do sistema como: 
R(z) 

E — — 

RICO AÇO 


o Aplicando-se o teorema do valor final: 


— IDR 
(DOT 
z=1 
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(a N 
Е. ERRO DE REGIME PERMANENTE 


o O limite existe se todos os termos (z — 1) no denominador se cancelam, 
dependendo do: 
e a) Sinal de referência 


e b) Ganho de malha aberta 


o Pode-se reescrever o ganho de malha aberta como: 


N(z) 


Sana "79 


n 


e Com N(z) e D(z) não possuindo raízes unitárias. 
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Е. ERRO DE REGIME PERMANENTE 


мттитонр 


o Tipo do sistema: 
• É o número de polos unitários na função de transferência discreta do 


sistema. 


o Assim, partindo das equacóes anteriores, tem-se: 


_ (z - 1)R(Z) C 
(O) = 1025) - O 
| (z- 1" D(2RG) 
e) = xa x Na) + G — D"DG)|,_, 


(z — D""D(DR(Z) 
N(1) + (z - 1)*(1)|,_, 
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е(со) = 


Е. ERRO DE REGIME PERMANENTE 


мттитонр 


o 1) Degrau Amostrado: 


7 
R(z) = 


7—1 


e(00) — K,-L() 


— ar. 


A+] A 


o A constante de erro de posição (K,) é finita para sistemas do tipo O e 


infinita para sistemas do tipo 1 ou superior: 


1 


е(оо) 2414 L(1) 
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= ERRO DE REGIME PERMANENTE 


мттитонр 


o 2) Rampa Amostrada: 


TE T 
CER: 


T 1 


CE) MA 


1 
— K, = = (z — DL(z) 
T z=1 
o À constante de erro de velocidade (K,) é nula para sistemas do tipo 0, 
finita para sistemas do tipo 1 e infinita para sistemas do tipo 2 ou 


superior: 
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Е. ERRO DE REGIME PERMANENTE 


мттитонр 


o 3) Parábola Amostrada: 


e De maneira análoga, tem-se: 
K : ( 1L 
= — (Z — 7 
а T2 ( ) 221 


о А constante de erro de aceleração (K,) é nula рага sistemas do tipo O ou 
1, finita para sistemas do tipo 2 e infinita para sistemas do tipo 3 ou 


superior: 


С ^N 


7 
ERRO DE REGIME PERMANENTE (RESUMO) 


TABLE 4-5 STATIC ERROR CONSTANTS FOR TYPICAL CLOSED-LOOP 
CONFIGURATIONS OF DISCRETE-TIME CONTROL SYSTEMS 


Closed-loop configuration Values of Kp, K,, and Ka 


K, = lim GHG) 


K, = lim 
z—1 


— ;-1y2 
К, = tim @—#-У©Н@) 
2=1 T? 


(1 —z-1)GH(z) 
T 


Kp = lim G(2)H (2) 


K, = lim 


2-1 


E =m a — 


z=) 


(1 —z-')G(z)H(z) 
T 


ER lim С,(2)НС (т) 


= (1 —z "DG (2)HG 2) 
‚ = lim T 


z-1 


K. = lim (1 —z12G,(z)HG.;(z) 


z—1 T? 


Kp = lim Gi(z)Gdz)H(z) 


K, = lim 1 —z73)6Gi(z)G»(z)H (z 


2-1 T 


К, == lim Ë — Ga(z)H(z) 


2=1 
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= EXEMPLO 


мттитонр 


o Para o sistema com as funções de transferência pulsadas mostradas а 
seguir, calcule o erro estático para uma entrada: 
• а) degrau amostrado 


e b) rampa amostrada 


|. . K(z*a) ú 
66) = DG D С) = еа 


Kc(27b) 0 <a,b,c < 1 
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- EXEMPLO 
o Solucáo: 
e O ganho de malha aberta é dado por: 
L(z) = C(z)G(z) = ССО 
(z — 1) (х5 ^ с) 
e Qual o tipo do sistema? 
o Tipo 1 
e Para entrada degrau: 
o K, = L(1) = oo, е(со) = 0 
e Para entrada rampa: 
K, => DO e) = roses 7 xe ss) 
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мттитонр 


ES 


Z 


o Considere a seguinte exponencial amostrada e sua transformada Z 


associada, com p sendo real ou complexo: 


рк, k = 0,1,2,.. Z — 


2-р 


о O valor da sequência, рага k suficientemente grande vale: 


0 , |р|<1 
Iplh 241 , Ipl=1 
о, |р|>1 
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E 


o Considere agora uma sequéncia do tipo: 


Z 
Z — Pi 


n n 
fü)-5 Ар,  k-012. eF()-» A 
[21 i-1 


o Pode-se concluir que: 
e À sequência decai exponencialmente quando todos os pólos estão 


dentro do círculo unitário aberto (dentro do círculo). 


e Para o caso de múltiplos pólos em cima do círculo unitário, a 


sequéncia será ilimitada. 
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ы 
E, LOCALIZAÇÃO DE POLOS ESTÁVEIS NO PLANO- 


Circulo Unitario 


Dé 


Região ' 
Estavel Re[z] 


Região 
Instavel 
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- ESTABILIDADE DE SISTEMAS DIGITAIS 
o Critério BIBO (Bounded Input — Bounded Output) 
o Definição: 

• Um sistema possui a propriedade de estabilidade externa se toda 
sequência de entrada limitada produz uma sequência de saida 
limitada. 

o Lema: 
e Um sistema linear, discreto e invariante no tempo, com resposta 


impulsiva g(k) é “BIBO estável” se, e somente se: 


+00 
> 0001 < + 
k=0 
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ES 


e ^ 
Be _ EXEMPLOS 


o 1) Determine se o sistema abaixo é ou não estável. 
G (z) 
U(z) Z Y (z) 
z=l 
o Solução: Integrador 


e Da definição de FT discreta, sabe-se que G(z) é a resposta 


impulsiva do sistema. Assim, para se determinar g(k) pode-se usar a 


divisão longa: z =: 
szyi ltz «EE. 
0+1 
-ip 
047" 
sg? 


0+2? 
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Е. EXEMPLOS 


мттитонр 


o Solucáo: 
e А resposta impulsiva g(k) vale então: 


g(0) =g(1) = 9(2) = 9(3) = ~: = 1 


e Е tem-se: 
+00 +00 
> ladol= > 1= + 
k=0 k=0 


e CONCLUSÃO: 
o Sistema Instável! 
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Ë EXEMPLOS 


мттитонр 


o 2) Determine se o sistema abaixo é estável: 


z+1 


e. E аа Н 
(2) = A062 401 


o Solucáo: 
e Os pólos de G(z) são as raízes do denominador: 


—0,6 + /—0,04 
z2 + 0,6z + 0,1 = 0 212 = ——=———= —0,3 + j0,1 
— 2 
A= 0,6? — 4 * 0,1 = —0,04 |z, | = [21 = 40,32 + 0,12 = 0,3162 < 1 


e Conclusão: quando as raízes têm módulo menor que 1 (estão 


LÁ 


dentro do círculo unitário), o sistema é “BIBO estável”. 
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A 
E. CRITÉRIO DE ROUTH-HURWITZ 


o Náo pode ser usado diretamente para sistemas discretos. 


o Deve-se utilizar a Transformação Bilinear: 


1-40 _ 2—1 


Z = c w = 
1— w z+1 


e transforma a região interna do círculo unitário do plano-Z no semi- 
plano esquerdo do plano-s. 


o Para um polinômio de ordem n: 


_ а к=? ju itoy 
F(Z) = anz” + а„—12? 1 + ... + Go | „(25 + an1 up t üg 


e TRABALHOSO! 
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Н. CRITÉRIO DE JURY 
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о Estuda а estabilidade de sistemas discretos sem a necessidade de 
determinar os polos. 


o Considerando que o polinómio característico do sistema é dado por: 
F(z) = anz” + а„—17? 5 L... + aZ + ao = 0, аһ > 0 


o As raízes de F(z) estão dentro do círculo unitário se: 


(1) F(1)»0 

(2  (C1YrF(-1)»0 
(3) lag| < an 
(4) Ibo] > |bn-1| 


(5) [Col > Icn-2| 


O mn» mn 
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CRITÉRIO DE JURY 


o Onde as n + 1 condições são calculadas a partir da tabela: 


LINHA 20 P E " | жап > sã z 
1 As а, а, i. e. s Ant а, 
2 а, Ana "P 24. а, а, а 
3 bo b, b, b, + Des 
4 р, 1 b, 2 [P b, Do 
5 Co €i C2 Cn-2 
6 C2 Сыз 6:4 Со 
2п-– 5 So Sı $2 $5 
2n-4 $3 5 Sı So 
2n-3 Yo r, r; 
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2 
Be _ CRITÉRIO DE JURY 


o Com os valores da tabela sendo calculados por: 


Ao An—k 
by = | | k =0,1,..,п—1 
An ак 
b é _ 
ск = |, ° | k=0,1,..,n—2 
bn-1 bg 
" zi M x” = [0 ll to Щ 
0 sz Sol 3 [53 sl Z Z |53 S; 
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- CRITÉRIO DE JURY 


o Observacoes: 
e 1) A primeira linha da tabela de Jury é formada pelos coeficientes 


de F(z) em ordem crescente. 


e 2) O número de linhas da tabela vale (2n — 3), sendo sempre 
impar. Os coeficientes de cada linha par são os mesmos da linha 


impar imediatamente acima, com a ordem invertida. 
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2 ^N 
- CRITÉRIO DE JURY 


o Observacoes: 
e 3) As condições 3 a (2n — 3) são calculadas usando o coeficiente da 
1% coluna juntamente com o último coeficiente da linha anterior. O 
coeficiente do meio da última linha não é usado e não precisa ser 


calculado. 


e 4) Condições 1 e 2 são calculadas usando F(z) diretamente. Se uma 
dessas condições não for atendida, conclui-se que o sistema possui 
raízes fora do círculo unitário sem a necessidade de se construir a 


tabela de Jury ou de testar as outras condições. 
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E: _ 
ми 


o Tabela de Jury: 


CRITÉRIO DE JURY 
EXEMPLO 


o Verifique a estabilidade do polinómio F(z): 


F(z) = 2º + 2.67! — 0.56z? — 2.0522 + 0.0775z + 0.35 = 0 


LINHA 2º 2! 2? “= 2^ E 
1 0.35 0.0775 -2.05 -0.56 2.6 1 
2 1 2.6 -0.56 -2.05 0.0775 0.35 
3 -0.8775 -2.5729 -0.1575 1.854 0.8325 
4 0.8325 1.854 -0.1575 -2.5729 -0.8775 
z 0.0770 0.7143 0.2693 0.5151 
6 0.5151 0.2693 0.7143 0.0770 

-0.2593 -0.0837 -0.3472 


7 
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E CRITÉRIO DE JURY 
"ges EXEMPLO 


o 6 Testes (n + 1) — ver slide 98: 
e Condições 1 e 2 — diretamente usando F (z): 
1. F(1) 2 1 + 2.6 — 0.56 — 2.05 + 0.0775 + 0.35 = 1.4175 > 0 
2. (—1)°Е(—1) = (—1)(—1 + 2.6 + 0.56 — 2.05 — 0.0775 + 0.35) = —0.3825 < 0 


e Condições 3 a 6 — usando a tabela de Jury 
3./0.35| « 1 
4. |-0.8775| > |0.8325] 


5. [0.0770] < [0.5151 qm 
6.|-0.2593| < |-0.3472| «ааа 


o Condições 2, 5 e 6 violam o critério de Jury 


B ° Sistema Instável p = (z — 0.7) — 0.5)(z + 0.5)(z + 0.8)(z + 2.5) = 0 
is 
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la 
- ANÁLISE DE ESTABILIDADE PELO MATLAB 


o Declarando o polinómio e calculando as raízes: 
>> den = [1 2.6 -0.56 -2.05 0.0775 0.35] 
den - 
1.0000 2.6000 -0.5600 -2.0500 0.0775 0.3500 

>> roots(den) 
ans = 

-2.5000 

0.7000 

0.5000 


-0.8000 
-0.5000 
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n 
m ANÁLISE DE ESTABILIDADE PELO MATLAB 


o Outras funções: 


o »» Zpk(g) 
e g — declarado na forma de função de transferência 


o >> den = [1.0 0.2 0.0 0.4] 
o >> ddamp(den, 0.1) 


Eigenvalue Magnitude Eguiw. Damping Equiv. Freg. (rad/s) 
3.05e-001 + 6.33e-001i 7.O3e-001 3.ü00e-ü001 1.i18e4ü0ü01 
J.U55-DDl = 6.33e-001i т. 03e-001 з. 00e-001 1.18e+001 

-5.10e-001 o. 10e-001 6. “de-002 3.15e+001 
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= EXERCÍCIO 1 


мттитонр 


o Usando o critério de Jury, determine а faixa de valores para o ganho 
K, de tal forma que o sistema operando em malha fechada e com 


realimentacáo unitária permaneca estável. 


K 


6G) = ас 12) 


e Solucáo: 


0,12 < K < 0,52 
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o Determine as condições de a, e a, рага que o sistema com o 


EXERCÍCIO 2 


polinómio característico igual a D(z) seja estável: 


o Solucáo: 


D(z) — 22 + 412 + а» 


e Condições: 
а, T 
Regiáo que garante a 
E diti estabilidade de G(z), 
1 segundo os valores 
| de dj e 23. 
a, >-a,-1 
2 1 a, 


© Eng. de Controle e Automação – Rejane de Barros Araújo/André Ferreira 


O (| 
m 
Е 
EE 
un 


INSTITUTO FEDERAL DE 
EDUCACÁO, CIÉNCIA E TECNOLOGIA 
PARÁ 


Campus Belém 


Disciplina: Engenharia de Controle I! 
Unidade Il: Projetos de Sistemas de Controle 
Digital 


Prof.? Dr.? Rejane de Barros Araújo 


INSTITUTO FEDERAL DE EDUCAÇÃO, CIÊNCIA E TECNOLOGIA DO PARÁ 
É. CURSO DE ENGENHARIA DE CONTROLE E AUTOMAÇÃO 


METITUTO FED 


Referência Bibliográfica 


* FADALI, M.S. Digital Control Engineering: Analysis and 
Design. Academic Press, 2009. 


« FRANKLIN, G.F.; POWELL, J.D. & WORKMAN, M.L. Digital 
Control of Dynamic Systems. Addison Wesley, 1998. 


* ASTROM, K. J; WITTENMARK, B. Computer Controlled 
Systems: Theory and Design, Dover Publications, 2011. 


INSTITUTO FEDERAL DE EDUCAÇÃO, CIÊNCIA E TECNOLOGIA DO PARÁ 
É. CURSO DE ENGENHARIA DE CONTROLE E AUTOMAÇÃO 


Agenda 


* Análise dinámica de sistemas discretos 
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LOCALIZAÇ 


ES Eng. de Cor 


LOCALIZAÇÃO DOS POLOS X RESPOSTA TEMPORAL 
- SISTEMAS DISCRETOS 


Imizi 


E ESPECIFICAÇÕES DE SISTEMAS DE CONTROLE 


* Erro de regime permanente 
e É a precisão de rastreamento em regime permanente que o sistema 


deve ter. 


<+ b) Resposta dinâmica 
e É a precisão que o sistema deve ter durante o período transitório. 
Estabilidade; 


Tempo de subida; 


e 1 

ө 2 

e 3) Máximo overshoot ou sobre-sinal; 
4 


) 
) 
) 
) Tempo de acomodação. 
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E. ESPECIFICAÇÕES DE SISTEMAS DE CONTROLE 


cir) Toleráncia aceitável 
| EN | 
M | Y | 0,05 
р | 53 m ШЕ 
| x M zT RE A 
En ро Pt 002 
MEE | 
La | | 
і | і | 
| Ë | 
| | | 
05 ее 5 | | 
| ; | | 
| i | 
| | | | 
| | | | 
i 5 | 
i l | | 
| | | | | 
ü Í 
fp — 
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E ESPECIFICAÇÕES DE SISTEMAS DE CONTROLE 
uu TEMPO CONTÍNUO 


_ T$ 
PO=e 1—82 100) (sobre-sinal percentual - M,(%)) 


TU TU 
T. = m С — : 
р (tempo de pico) 
2 
2 Od ш„у1—ё 
L(s) = => 
52 + 28025 + 0%, 4 

T; = — (tempo de acomodação — 2%) 

Sw 


z l4 — ¿2 — 
п — areseny 1 — $^ — n—p (tempo de subida) 
Wn 1 — È? Wa 
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E ANÁLISE PARA SISTEMAS DE 2º ORDEM 


4 Casos Analisados: 


e 1) ¿w, constante (tempo de acomodação constante) 
e 2) é constante (PO constante) 


e 3) éw, constante (tempo de subida constante) 
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E. ANÁLISE PARA SISTEMAS DE 2º ORDEM 


* Pólos Dominantes (sistema sub-amortecido): 

5= 0 Ё јод = —$0 tjog 41—6^,  0<$<1 
4$ Lembrando que z = e'f, tem-se: 
_ -Eontion E) y = ¿Cor (tion 1-£2)T 


Z 
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E. ANÁLISE PARA SISTEMAS DE 2º ORDEM 


* Caso 1: £o, constante 


** Nesse caso, os pólos no domínio Z são dados por: 


y = gC Ei) (tion 1-£2)T 
L— 


modulo | 


e O módulo de z será: |z| = e *“n” = constante 


e A fase de z é dada рог: Zz = +o, T./1 — č? , que será qualquer valor. 
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E. ANÁLISE PARA SISTEMAS DE 2º ORDEM 


* Portanto, no plano Z tem-se os seguintes lugares geométricos que são 


circulos centrados na origem de raio 
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E. ANÁLISE PARA SISTEMAS DE 2º ORDEM 


* Caso 2: é constante 


z = eTSUnTetijwnTy 1-¿? 


% Se: 
<=] > LZz=0º e |z| > semi-eixo real 0 < z € 1 
¿=0 >  lz|=1 е Zz=V >° circulo unitário 


O<£<1 > | wn = 0 > |z|=1 e z=0 


— 


Om crescente АНТИ е Z Zcrescente 
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E ANÁLISE PARA SISTEMAS DE 2º ORDEM 


* Este lugar geométrico é uma espiral que inicia no ponto 1. Portanto, o 
lugar geométrico no plano Z que corresponde a um determinado 


coeficiente de amortecimento são espirais: 


Curvas 
espirais 


ЄЎ Observacáo: А figura acima está fora de escala. 


E ANÁLISE PARA SISTEMAS DE 2º ORDEM 


* Caso 3: о, constante 


7 = e $ ©пТ go tjonT, 1— ç 


¿=1 ә [| = e ?nl e Zz=0" > ѕеті-еіхогеа1 0 < z € 1 


¿=0 > |z|-1 е ZZ = wnT > circulo unitario 


@ Eng. de Controle e Automação – Rejane de Barros Araújo/André Ferreira 


E. ANÁLISE PARA SISTEMAS DE 2º ORDEM 


* Lugar geométrico para w, constante: 
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E. ANÁLISE PARA SISTEMAS DE 2º ORDEM 


“ Lugar geométrico para as especificações de desempenho no plano Z: 


IN 


cu 
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LOS 
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+ 
IE 
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d 
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ÉE _ RESPOSTA AO IMPULSO NO DOMÍNIO Z 


% O posicionamento do pólo no plano Z fornece as características de 


resposta temporal: 
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LUGARES GEOMÉTRICOS NO PLANO Z 


' = constante 
(POW) 


| Xe a 


АА A 


ES 
PAN 

tw 

Ü 02 


22 


. Cuwasde б =constante e Wy=constante Т = Período de amostragem бм é constante 


ж (5) 


_ ; Т 
Curvas de G= АА = constante: circulos de raio f = 5 n 


E Ñ = = Maior círculo — menor 7. 
Curvas de W¿= constante: raios com inclinação g=w, rad) E w,T p/ £20) ч 


Maior Wn = menor №. 
Maior © этели PO% 
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- LGR NO DOMÍNIO Z 


« A equação característica do sistema е dada por: 
1-4 C(z)G(z) = 0 


e onde C(z) é a função de transferência do controlador e G(z) é a função 


de transferéncia do conjunto “A/D + planta analógica + D/A”. 
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- LGR NO DOMÍNIO Z 


% Supondo que C(z) é formado por um ganho que multiplica uma função 
de transferência em z, a equação característica pode ser colocada na 


forma abaixo, onde L(z) é o ganho de malha aberta: 
1+KL(z) = 0 


4 É a mesma forma obtida рага o LGR no domínio s: 
e Regras para plotar o LGR sáo as mesmas. 
e Pode ser obtido utilizando o MATLAB (ou software com uma funcáo 


equivalente). 
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E LGR NO DOMÍNIO Z 
ge EXEMPLO 


* 1) Usando o LGR, obtenha o ganho crítico para o sistema abaixo: 


1 
L 7-7 


* Solução: usando o MATLAB 


e num = 1; 
e den = [1 -1]; 
e ts — 0.01; 


e G = tfí(num,den,ts); 


e rlocus(G) 
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LGR NO DOMÍNIO Z 
ges EXEMPLO 


E Figure 1 
File Edit View Insert Tools Desktop Window Help 


пыга aana 2 ПВ ац 


Root Locus 


„= 


- 


f 


System: g 
Gain: 2 
Pole: -1 
: Damping: -0.000568 
+ Overshoot (%): 100 
‚ Frequency (radisec): 314 


Imaginary Axis 
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E LGR NO DOMÍNIO Z 
ue EXEMPLO 


* 2) Usando o LGR, obtenha o ganho crítico para o sistema abaixo: 


1 


LG) = 01-05) 


% Solução: usando o Matlab 


e num = 1; 
e den = conv([1 -1],[1 -0.5]); 
e {5 = 0.01; 


e G=tf(num,den,ts); 


e rlocus(G) 
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LGR NO DOMÍNIO Z 
ges EXEMPLO 


> 


ElFieure1 (Ja) 


File Edit View Insert Tools Desktop window Help 


$a | 82479408 


Root Locus 


System: g 

Can: D.5 

Pole : 0.75 + 0.6621 
Damping: -0.000188 
Overshoot (%): 100 
Frequency (radisec): 72.3 


E 
2 
> 
m 
E 
Lj 
= 


E IMPLEMENTAÇÃO DIGITAL DE CONTROLADORES 
=> ANALÓGICOS (EQUIVALENTES DISCRETOS) 


« Efetua-se inicialmente o projeto de controlador no plano s, utilizando as técnicas ја 
conhecidas para sistemas contínuos. 


% Utiliza-se o mapeamento de polos e zeros do plano s para o plano z (z=e*), 
encontrando-se o controlador discreto equivalente ao contínuo (emulação). 


« Vantagem: 
e Possibilidade de se utilizar as técnicas de projeto já estudadas para sistemas 
contínuos. 


* Desvantagens: 
e Esse método ignora totalmente o fato de que os conversores A/D e D/A e o 
microcomputador serão utilizados. 
e Há a necessidade de que o período de amostragem 7 seja o menor possível para 
© evitar discrepáncia em ге!асао ao projeto contínuo. 
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E CONTROLADORES ANALÓGICOS 
aee EXEMPLO 


* Considere o sistema G(s) operando em malha fechada e com 


realimentação unitária. 


1 


G(s) = Os + D 


* Projete um controlador analógico e, utilizando um equivalente discreto, 
obtenha um controlador discreto que satisfaça as seguintes condições: 


o PO < 16% её, < 105 
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E CONTROLADORES ANALÓGICOS 
aee EXEMPLO 


« Solução: Projeto de controlador discreto usando emulação (equivalente 
discreto). 


*Mpdes = 0.16; 

* Csi = - log(Mpdes)/sqrt(pi^2--log(Mpdes)^2) 

“* tsdes = 10; 

“+ wn = 4/(csiíTsdes) 

** polodes = -csi*wn + j*wn*sqrt(1-csi^2) YPolo desejado 
* g1 = tf(1,[10 1 0]) % Sistema não-compensado 

4 rlocus(g1), axis([-0.12 0.02 -0.2 0.2]) 
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E CONTROLADORES ANALÓGICOS 
gre EXEMPLO 


*” LGR do sistema não-compensado 


m Figure 1 Qu 
м 


File Edit View Insert Tools Desktop Window Help Posicáo para o pólo 


Dia 1229/12 LE st desejado чаш 


RootLocus 


System: 91 
Gain: 0.0995 
Desa 
ng: 0.501 i 


16.2 ' 
Frequency (rad/sec): 0.0998; 


3 
> 
E 
? 
Е o 


*” Т ОК não passa pelo ponto desejado: 


O e o controlador não será apenas proporcional 
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E CONTROLADORES ANALÓGICOS 
aee EXEMPLO 


** Define-se o controlador C(s): 


5+7 
C(s) === 
S + pı 


* Fazendo z, = 0.1, pode-se escrever: 


s+01 1 _,s+01 01 
s+p11082+s ^ s+pņp s(s + 0,1) 


KC(s)G(s) =K 


0,1 
S(S + pı) 


O Eng. de Controle e Automação – Rejane de Barros Araújo/André Ferreira 


KC(s)G(s) = K 


E CONTROLADORES ANALÓGICOS 
m= EXEMPLO 


** Lembrando que, para determinado ponto pertencer ао LGR, é válida a 


> 6, - > 6, = —180? 


relacáo: 


* Para KC(s)G(s), tem-se: KC(s)G(s) = K A EE 
** Polo desejado: sg,; = —0,4 + j0,6857 
« Assim: 

УД; — “ру = —180? 


s=—0,4+j0,6857 
“ру = 180º — 120,2562° Zp, = 59,7438? 
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E CONTROLADORES ANALÓGICOS 
ges EXEMPLO 


** Desta forma, pode-se determinar o valor de р;: 


0,6857 
59,7438? 


Padi 


pr” 0,4 


0,6857 


tg(59,74389) = 
gl ) ARA 


=> p, = ӨН 
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E CONTROLADORES ANALÓGICOS 
e EXEMPLO 


** LGR do sistema compensado 


° K = 6,23 E Figure 1 Jae 


File Edit View Insert Tools Desktop Window Help 


paga 15220908 вп 


Root Locus 


System: g2 

Gain: 6.23 

Pole: -0.4 + 0.681 

Damping: 0.507 

Overshoot (965: 15.8 
Frequency (radisec): 0.789 


Imaginary Axis 
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E CONTROLADORES ANALÓGICOS 
2- EXEMPLO 


€ Sistema compensado (К = 6,23): 


Controlador Planta 


ES |ks+0oD| VC) 1 Ү (5) 
L — s(10s +1) Е 


% Reposta ао degrau сот К = 6,23 (Simular no matlab) 


* Reposta ао degrau сот К = 5,8 (Simular no matlab) 
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E CONTROLADORES ANALÓGICOS 
2- EXEMPLO 


* Então, a função de transferência do controlador que atende as 


especificações desejadas é: 


s+ 0,1 
s+ 0,8 


KC(s) = 5,8 


4 Fazendo o mapeamento de polos e zeros usando z = e?" e supondo 
T = 0,1 s: 


га =e 91912099 ра =e 080% = 0,9231 


« Assim: 


Z= Zi 2 — 0,99 
e KaCa(z) = Ka > с КаСа(2) = Ka — 9231 
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E CONTROLADORES ANALÓGICOS 
aee EXEMPLO 


* O ganho DC do controlador contínuo deve ser igual ao do controlador 


discreto: 


ganho DC = lim(KC(s)) = lim(KaCa(z)) 
S— Z— 
“ Desta forma: 


ú 5,8 * 0,1 * (1 — 0,9231) 


E = 5,5752 
d 0,8 * (1 — 0,99) 
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E CONTROLADORES ANALÓGICOS 
e EXEMPLO 


* Resposta ao degrau dos controladores: 


LIT JER Brew JOB 
Elle Edit View Insert Tools Desktop Window Help x Eile Edit View Insert Tools Desktop Window Help x 
Dek б &&9 9 ж 08 = El beds | Raro É DE stu 


Step Response Step Response 
T T T T T T 


Time (sec) Time (sec) 
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E CONTROLADORES ANALÓGICOS 
aee EXEMPLO 


% Sistema com controlador digital obtido através de um equivalente discreto: 


| | 
l Planta 


r(t) +m e(t) e(kT) и(КТ u(t) 1 y(t) 


* Obtendo o equivalente discreto usando o MATLAB: 


Resposta ao Degrau (sistema discretizado): (z — 0,99) 
G (z) = 5,5752 ——- 
(z — 0,9231) 
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m _ CONTORNOS NO PLANO Z 


4* Polos dominantes: 
F(s) — —SWn + jog 


F(z) —> eC $9ntjoq)T 


4* Para o caso discreto: 


212 = e $9nT 7 + Qa T 
(z = e CSuntiwa)T )(z = еС$®п-—/Ф®а)Т ) — 72 = 2соѕ(оаТ)е "Т7 + e ?$onT 
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m _ CONTORNOS NO PLANO Z 


Z1,2 — еп Z + Wal 


* оа constante no plano s e no plano Z: 


jæ Im(z) 
3 


3л/4 


ba 
т 


- m4 lr 
0 
—Л/4 4 
—2 
—3л/4 
3 G = 


(а) (b) 
e No plano Z: íz = wgT constante (linha radial) 
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Ë мттлонр 
ми 


CONTORNOS NO PLANO 2 


21,2 = esonT Z T Wal 


“* Parte real (-¿w,,) constante no plano s e no plano Z: 


(a) 


Im(z) 
3 


2 3 Relz) 


(b) 


e No plano Z: |z| = e *%nT constante (círculo) 
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ЭЖ CONTORNOS NO PLANO 2 


21,2 = esonT Z + Wal 


*» £ constante: 
e Espirais logarítmicas que ficam menores quanto maior for o valor 


do coeficiente de amortecimento. 


.. $0 -£(79"/180:) 
= 1-47 == 1—2 
|z] = е е V 


e onde |z| é a magnitude do polo e 0 sua fase. 


** wp constante: 


iz] - e (cT)? —0* 
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- CONTORNOS NO PLANO 2 


* wn e é constantes 
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m _ CONTORNOS NO PLANO Z 


* Observações: 
e 1) Para cada valor de é há 2 espirais que são imagens com relação ao 


eixo real, correspondendo aos ângulos positivos e negativos. 


e 2) Para uma dada espiral, a magnitude do polo cai logaritmicamente 


com seu ângulo. 


e 3) Para um mesmo ángulo 0, o aumento do coeficiente de 


amortecimento significa uma magnitude menor para o polo. 
e 4) Todas as espirais começam em 0 = 0 e |z| = 1. 
e 5) Para determinado é e |z|, o ángulo do polo é dado por: 


Je 
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- CONTORNOS NO PLANO 2 


** Usando o MATLAB: 


e zgrid(csi,wn) 


Paraé£203e207: 
e zgrid(0.3,[ ]) 
e zgrid([0.3 0.7],[ ]) 


- L | L 
-1 -0.8 -06 -. 
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- CONTORNOS NO PLANO Z 


* Para o, = 0,5 € o, = 0,7 : 
e zgrid([ ],0.5) 
e zgrid([ ],[0.5 0.7]) 
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- CONTORNOS NO PLANO 2 


* Para obter o contorno completo, use 'zgrid' sem parámetros: 


e ¿ (0,1 a 0,9) e 0, (7/57), сот N inteiro de 1 a 10 


Q.67/T | 0 
-1 | l l 0.5 l l l 
-1 0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 
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E: ESPECIFICACOES DE DESEMPENHO NO 
т DOMÍNIO DISCRETO 


* 1) Constante de Tempo 
e Constante de tempo para o decaimento exponencial da forma de 


onda amostrada. 


e O sinal amostrado náo é necessariamente igual a uma porcáo do 


valor final após uma constante de tempo. 


e É definida como: 


1 
T = —— 
Sw 
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E: ESPECIFICACOES DE DESEMPENHO NO 
т DOMÍNIO DISCRETO 


* 2) Tempo de acomodação: Para o critério de 2%, é dado por: 
4 


| = 
И? 


* 3) Frequência de oscilação o 
e Esta frequéncia é dada pelo ángulo do pólo complexo dominante 
dividido pelo período de amostragem. 


e Lembrar que: 


21,2 — e $enT7 T WaT 


* 4) PO (Mp (%)), о, e ¿.Análogos ao sistema contínuo. 
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B TRANSFORMAÇÃO BILINEAR 


* A relação abaixo é conhecida como transformação bilinear: 


7—1 
z+1 


S = C 


+ A partir da identidade z = e** e fazendo uma aproximação de primeira 


ordem (Regra de Tustin), chega-se em: 


Usando o matlab: help c2d 
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= _ METODOS DE DISCRETIZAÇÃO 


Backward Forward 


kT-T kT 


Trapezoidal 


kT-T kT 
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E: METODOS DE DISCRETIZAÇÃO 
ш TRAPEZOIDAL 


Plano s “Plano z 


E 
= т 2+1 
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E: METODOS DE DISCRETIZAÇÃO 


BACKWARD 


Plano Z 
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METODOS DE DISCRETIZAÇÃO 
FORWARD 


_ 2—1 
= т 
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PROJETO DE CONTROLADOR PROPORCIONAL NO 
are DOMÍNIO Z 


** O controle proporcional envolve a escolha de um ganho DC que atenda 


as especificações de desempenho desejadas. 


> Os polos dominantes em malha fechada devem estar localizados em 


uma posição adequada no plano complexo. 


* O projeto analítico de sistemas de baixa ordem é possível de ser 
realizado, embora seja mais trabalhoso do que o equivalente no 


domínio continuo. 
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PROJETO DE CONTROLADOR PROPORCIONAL NO DOMÍNIO Z 
m= EXEMPLO 


* Para o sistema G(s) abaixo, operando em malha fechada com período 
de amostragem T = 0,04 s e com realimentação unitária, projete um 
controlador proporcional discreto para que o sistema compensado 


possua as seguintes características: 


1 


GG) - S5) 


e a) е(о) = 10% para uma entrada rampa 
e b) č > 0,25 
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PROJETO DE CONTROLADOR PROPORCIONAL NO DOMÍNIO Z 
m= EXEMPLO 


« % Sistema analógico 


e num = 1; 

e den = [1 5 0]; 

e G=tf(num,den) 

e Ts = 0.04; 

e Gd = c2d(G,Ts) 


: _ 7.4923x10 ^ (z + 0.9355) 
e zpk(Gd) — a(z) = (z — 1)(z — 0.8187) 


Ө Eng. de Controle e Automação – Rejane de Barros Araújo/André Ferreira 


PROJETO DE CONTROLADOR PROPORCIONAL NO DOMÍNIO Z 
m EXEMPLO 


** Lembrando que: 


(z == е(75%п+1%аТ) (y == eC $9n-J0q)T) — 72 — 2cos(waT)e ° 9nTz + e 25 ФһТ 


* A equação característica do “sistema discreto + controlador proporcional" 


é dada por: 


1 + KGzas(z) = z? — (1.8187 — 7.4923x10"*K)z + 0.8187 — 7.009x10 ^K 


= z? —2cos(waT)e z+e 


* Igualando-se os coeficientes, existem 3 incógnitas (£, w, e К) e apenas 2 
equações (z! e z?) 


o Utiliza-se as especificações de projeto! 
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PROJETO DE CONTROLADOR PROPORCIONAL NO DOMÍNIO Z 


aee EXEMPLO 
7.4923x10 *(z + 0.9355) 
Ca (Z) = ——— N 
(z — 1)(z — 0.8187) 
% Qual o tipo do sistema? 
e Tipo 1 


* Lembrando que: 


K, = L(1) K, = =G — 1)2L(z) 


z=1 
a 
К» = TU = DL(z) 


z=1 
e Com L(z) — ganho de malha aberta 
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PROJETO DE CONTROLADOR PROPORCIONAL NO DOMÍNIO Z 
m= EXEMPLO 


7.4923x107*(z + 0.9355) 
(z — 1)(z — 0.8187) 


4* Constante de erro de velocidade: 


Ga(z) — 


1z-—1 
K, = KG(z) 
T 
z=1 
—4 
К, = К 74923410 (1+0.9355) _ К 


(0.04)(1—0.8187) 5 


* Рага um erro de 10%, tem-se: 


K 100 100 
== Ky === =10 — a K=50 
5 " e(0)% 10 
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PROJETO DE CONTROLADOR PROPORCIONAL NO DOMÍNIO Z 
m= EXEMPLO 


% % LGR 

$* rlocus(gd) 

$* zgrid(0.25,[]) 

« axis([-1.1 1.1 -1.1 1.1]) Roos 


0.8 
0.6 
0.4 - 

E 02 System: gd 

> Gain: 59.6 

© 0 

c H š 

5 | Damping: 0.25 

Ë 0.2) Overshoot (96): 44.4 

= ` \ Frequency (rad/sec): 7.67 | 
0.4 - / 
-0.6 
-0.8 


Real Axis 
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PROJETO DE CONTROLADOR PROPORCIONAL NO DOMÍNIO Z 
m= EXEMPLO 


* % Para erro estático de velocidade de 10% — K >= 50 - A partir do LGR 
+ K = 59.6; 

* % Malha Fechada 

* mfd = feedback(K"gd,1); 

*$ [nd,dd] = tfdata(mfd,'v'); 

« ddamp(dd, Ts) 


* % Como a especif. de amort. não é atendida, deve-se alterar o ganho 
4 K = 58; 

« mfd = feedback(K"gd,1); 

« [nd,dd] = tfdata(mfd,'v'); 

“• ddamp(dd,Ts) 
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PROJETO DE CONTROLADOR PROPORCIONAL NO DOMÍNIO Z 
m EXEMPLO 


Resposta ao degrau (č = 0,25 — PO = 44,4%): 


step Response 


| System: mid 
Peak amplitude: 1.44 
| Cwershoeot (6): 44.2 


‚ XAt time (sec): 0.44 


Amplitude 


Ü 0.5 1 1.5 2 2.5 d 


Time (sec) 
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PROJETO DE CONTROLADOR PROPORCIONAL NO DOMÍNIO Z 
тее EXEMPLO 


* Resposta a Rampa: 


Step Response 
T 


100 | | I I F 
Rampa de referência 


Sistema compensado 7 


Amplitude 


| 
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 


Time (sec) 
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E PROJETO DE CONTROLADOR DIRETAMENTE NO 
= DOMÍNIO Z 


« Vantagem sobre os equivalentes discretos: 
e Os polos do controlador náo ficam restritos ao mapeamento do 


plano s para o plano z. 
* Desvantagem: 


e Um controlador digital com um zero e nenhum polo não é 


realizável. 
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PROJETO DE CONTROLADOR DIRETAMENTE NO DOMÍNIO Z 
m EXEMPLO 


% 1) Projete um controlador digital com período de amostragem de 0,15 
para um motor DC com funcáo de transferéncia G(s). O sistema em 


malha fechada deve atender às seguintes especificações: 


1 


G(s) = TDS + то) 


e a) erro nulo em regime permanente рага uma entrada degrau. 


a 
e b) coeficiente de amortecimento de 0,7. 


e c) tempo de acomodação de 1s. 
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PROJETO DE CONTROLADOR DIRETAMENTE NO DOMÍNIO Z 
m EXEMPLO 


« % Especificações 

«csi = 0.7; 

* Мр = exp(-csi*pi/sqrt(1-csi^2)); 
« Wn = 4/(csi*1); 

# Ts = 0.1; 


** % Polo discreto desejado 
** pdes = exp(-csi*WnºTs + j"Wn*sqrt(1-csi^2)* Ts) 


* % Sistema analógico 
«• g-tf(1,conv([1 1],[1 10])) 
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PROJETO DE CONTROLADOR DIRETAMENTE NO DOMÍNIO Z 
m EXEMPLO 


% % Sistema discreto 
# Ts = 0.1; 

‹ gd = c2d(g, Ts); 

** zpk(gd) 


4$ C(Z): 
e Polo em 7 = 1 e um zero a ser definido. 


e Tentativa: cancelar o polo em z = -0,9048 


< % Controlador 


e cont = tf([1 -0.9048],[1 -1], Ts) 
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PROJETO DE CONTROLADOR DIRETAMENTE NO DOMÍNIO Z 


qe EXEMPLO 
% LGR: 
Root Locus 
1.53 | | 


System: gc 

. | Gain: 35.5 

\ Pole: 0.62 + 0.264i 
Damping: 0.7 
Overshoot (96): 4.58 
Frequency (rad/sec): 5.64 


E \ 


\ 

\ 

| 

| | 
SO» 


| 
/ 


0.5 - 


Imaginary Axis 
o 


-0.5 - 


ERN | | | | | | | | | 
-3.5 -3 -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 


Real Axis 
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PROJETO DE CONTROLADOR DIRETAMENTE NO DOMÍNIO Z 
qu EXEMPLO 


* Resposta ao Degrau 


step Response 


system: mf 
Peak amplitude: 1.05 
Overshoot (96у: 4.52 
: system: mf 
Ads Setting Time (sec): 1.06 


Amplitude 


ü 0.5 1 1.5 2 2.5 


Time (sec) 
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PROJETO DE CONTROLADOR DIRETAMENTE NO DOMÍNIO Z 
qu EXEMPLO 


* Lembrando que: >; 0, — 3,0, = —180º 


C(z) = = 215 = e 599 ¿ + waT 


* Zdes = 0,6153 + j 0,2660 


> Então: 


Oz. + Z(z + 0,6945)|,-0,6153 + j 0,2660 — (Z + 0,9048)|,-0,6153 + j 0,2660 
— Z(z + 0,3679) |,_о,15з + 0,2660 — Z(Z — Dl,=0,6153 + jo,2660 = — 180º 


07. + 11,4806? — 137,4231º — 47,0795? — 145,3382º = —180? 


04. = 138,3602º 
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PROJETO DE CONTROLADOR DIRETAMENTE NO DOMÍNIO Z 
qu EXEMPLO 


Im(z) 


Re(z) 


0.2660 0,2660 
T" , ——— Z + 0,6153 
E cades "c tan(41,64º) 
_ z — 0,914 
z, — 0,914 —— pio = 
Ec 
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PROJETO DE CONTROLADOR DIRETAMENTE NO DOMÍNIO Z 


— 
é LGR: 
1.5 
0.5 


Imaginary Axis 
o 


-0.5 


-1.5 


EXEMPLO 


Root Locus 


System: gc2 

. Gain: 35 

` Pole: 0.616 + 0.265i 
'"Damping: 0.702 
Overshoot (%): 4.5 
Frequency (rad/sec): 5.7 


| | 
| | B — 


| 
° 
| 


| 


/ 


-5 


-4 


Real Axis 
Eng. de Controle e Automação – Rejane de Barros Araüjo/André Ferreira 


PROJETO DE CONTROLADOR DIRETAMENTE NO DOMÍNIO Z 


NETTUTOFED 
" EXEMPLO 
* Resposta ao degrau: 
Step Response 
14 — ii 
System: mf2 
es Peak amplitude: 1.04 
12 -  Ovwershoot (%): 4.87 
At time (sec): 0.8 System: mf2 


Setting Time (sec): 1.07 


4 Exercício: Time (sec) 
e Refazer o projeto, usando um período de amostragem menor e/ou modificando o pólo do 


controlador PD. 
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PROJETO DE CONTROLADOR DIRETAMENTE NO DOMÍNIO Z 
m EXEMPLO 


* 2) Para o sistema G(s) abaixo, projete um controlador digital que atenda 


as seguintes especificações em malha fechada: 


1 


G(s) = s(s + 1)(s + 10) 


e a) coeficiente de amortecimento de 0,5; 
e b) tempo de acomodação de 5 s. 
e Usar o programa projeto PD.m e verificar o efeito do período de 


amostragem no projeto do controlador digital. 
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É. PROJETO DE CONTROLE DIRETO 


* Suponha que a função de transferência em malha fechada de um 


sistema Ga (2) seja conhecida. 


“ Projetar um controlador C(z) para que o sistema em malha fechada 


tenha a função de transferência desejada. 
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É. PROJETO DE CONTROLE DIRETO 


+» Desta forma: 
C(z)G (z) 
1+ C(z)G (>) 


cl (Z) 


G (2) сон 8 — Ga (Z) 


Ga (Z) = 


C(z) = 


* C(z) deve ser realizável e fazer com que o sistema seja 


assintoticamente estável. 
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É. PROJETO DE CONTROLE DIRETO 


4* Causalidade: 


| N(gz| _ Na(gz ! 
G(z) — DG) Ga (Z) = “Did 
E a 1 Ga(z) | р() Na(z)z ! 
©) = E -G4G)]  N G7 Da) - Na 


* Características de um controlador causal: 
e Nümero de polos 2 Nümero de zeros; 
e Náo ter avanco de tempo. 


% Sistema em malha fechada e planta devem ter : 
p 
e Mesma diferenca entre nümero de polos e zeros; 
o e Mesmo atraso de transporte. 
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É. PROJETO DE CONTROLE DIRETO 


* Para uma planta de fase náo-mínima: 
e Não deve ser feita o cancelamento de polos e zeros instáveis para que 


se garanta a estabilidade assintótica. 


<% Para um sistema com zero de fase náo-mínima: 


G(z) = Gy(z)(z — Z), |2| > 1 


goes Co CC 
dC 1-C()G(z 1-+C(@2G,(z)G — 2) 


© e Zero de malha fechada em z 
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É. PROJETO DE CONTROLE DIRETO 


* Para uma planta instável: 
e Não deve ser feita o cancelamento de polos e zeros instáveis para 


que se garanta a estabilidade assintótica. 


* Para um polo instável: 


G(z) — ZA |2] > 1 — > 1-6 (0=1- C(z)G (z) 


(2—2) 1 + C(z)G (>) 
Е 1 Е 7—7 
1+C(2) A z—z-4C(z)G(z) 


e e Zero de malha fechada em z 
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Ë PROJETO DE CONTROLE DIRETO 


HITTUTOFID 
mu 


*» Erro em regime permanente: 
e Para erro nulo em regime permanente para uma entrada degrau, tem- 


Se: 


y (co) = lim(z — DGa(2)— = Ga(1) 


Z 
7—1 
Go (1)=1 
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E PROJETO DE CONTROLE DIRETO 
ge RESUMO 


% Restrições para G./(z): 
e Mesma diferenca entre nümero de polos e zeros para sistema em 
malha fechada e planta; 
e Polos instáveis da planta devem ser zeros de malha fechada; 


e Zeros instáveis da planta devem ser zeros de 1 — Ga (Z); 


e Para erro de regime nulo para uma entrada degrau é necessário 


que Ga(z) = 1. 
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É. PROJETO DE CONTROLE DIRETO 


* Desvantagens: 
e O maior problema desse método é a escolha adequada de uma 


função de transferência em malha fechada. 


e A escolha dos polos e zeros que atendam às especificações para 


sistemas de ordem mais elevada pode ser difícil de ser feita. 


e O desempenho deste tipo de controle depende muito da precisão do 


modelo utilizado para a planta. 
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É. PROJETO DE CONTROLE DIRETO 


* Para se obter a posição desejada para os polos em malha fechada, pode-se 
utilizar os polos de Bessel. 


e Levando em conta o tempo de acomodação desejado: 


Localização dos polos de Bessel para o tempo de acomodação igual a 1 s 


[D —— 
=i | И 
Dweisme| cmm | _ 
Lame a oras amens 
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É. PROJETO DE CONTROLE DIRETO 


Localização dos polos de Bessel para largura de banda de 1 rad/s 


—1,0000 
—0,8660 + j0,5000 


—0,5906 + j0,9072 | —0,8516 + j0,4427 —0,9246 
—0,5385 + 0,9617 | —0,7998 + j0,5622 | —0,9093 + j0,1856 
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—0,6573 + j0,8302 | —0,9047 + j0,2711 


E PROJETO DE CONTROLE DIRETO 
ge PROCEDIMENTO 


* Escolha o tempo de acomodação (ou largura de banda) desejado. 


*$ Escolha a localização dos n polos a partir da tabela e divida pelo tempo 


de acomodação (ou multiplique pela largura de banda). 


* Obtenha С. (2) convertendo a localização dos polos e zeros no plano s 


para o plano z, usando a transformação z = e*”. 


** Verifique se Ga(z) atende as condições de causalidade, estabilidade e 
erro em regime permanente. Caso contrário, modifique G.,(z) até que as 


o. ИИИ sejam atendidas. 
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E PROJETO DE CONTROLE DIRETO 
ge EXEMPLO 


« Usando o método de controle direto, projete um controlador discreto para 


G(s) que atenda às seguintes especificações: 


1 


CS) =G+D(G+10) 


e a) erro nulo em regime permanente para uma entrada degrau; 


e b) tempo de acomodação de 4 s. 


e OBS.: Usar período de amostragem de 0,02 5. 
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E PROJETO DE CONTROLE DIRETO 


EXEMPLO 


Solucáo: 

% Sistema analógico 
g = ti(1,poly([-1 -10])); 
zpk(g) 


Y% Taxa de Amostragem 
Ts = 0.02; 


% Sistema dicreto 
gd = c2d(g, Ts, zoh’); 
zpk(gd) 


88 
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E PROJETO DE CONTROLE DIRETO 
ge EXEMPLO 


*$ Sistema de 2? ordem sem polos ou zeros instáveis. 


** A partir das tabelas de Bessel, escolhe-se os polos levando em conta o 


tempo de acomodacáo desejado: 
(—4.0530 + у 2.3400)/4 


*$ Escolher o ganho para erro de regime nulo para uma entrada degrau: 


1,369 


G ETO = ----ЕЕЗЕРНИЕ 
als) = 2420275 + 1369 
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E PROJETO DE CONTROLE DIRETO 


EXEMPLO 


% Pólos desejados -> ver tabela de Bessel 
pdes1 = (-4.0530 + j*2.3400)/4; 

pdes2 = (-4.0530 - j*2.3400)/4; 

den = poly([pdes1 pdes2]); 

num = polyval(den,0); 

gcl = tfínum, den) 

% FT discreta com os polos desejados (tabela de Bessel) 
gclz = c2d(gcl, Ts); 

zpk(gclz) 

% Obtencáo do controlador direto 

с = 1/gd"(gclz/(1-gclz)); 

zpk(c) 
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E PROJETO DE CONTROLE DIRETO 
к= EXEMPLO 


* Resposta ao degrau: 


Step Response 


system: mft 
setiling Time (sec): 3.7 1 


ü 1 Es a 4 5 E 


Time (sec) 
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É. CONTROLE DEADBEAT 


% Sistemas de controle contínuo só atingem a saída desejada 
assintoticamente após um tempo infinito. 


* Sistemas de controle discreto são capazes de atingir o sinal de 
referência depois de um intervalo de tempo finito. 


* O controle deadbeat ou controlador de tempo finito é uma abordagem 
utilizada em controle digital para fazer com que o sistema atinja o sinal 
de saída após um tempo pré-determinado. 


* Baseia-se no controle direto, porém escolhe-se como função de 
transferência em malha fechada desejada G.,(z) = 27°. 
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É. CONTROLE DEADBEAT 


* O expoente k deve ser maior ou igual ao atraso intrínseco do processo 
discretizado, ou seja, deve ser maior que a diferença entre os graus do 


denominador e numerador da função de transferência discreta. 


*$* Desta forma, o controlador obtido vale: 


1 z^ 
G (z) Р =z 


C(z) = 
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É. CONTROLE DEADBEAT 


* Vantagem: 
e Projeto extremamente simples, pois a única variável de projeto é o período de 


amostragem T.. 


* Desvantagem: 
e Pode exibir um comportamento inadequado entre intervalos de amostragem 
(oscilações ocultas), porque o sinal de controle se mantém inalterado entre 


dois intervalos de amostragem consecutivos. 


e O sinal de controle pode facilmente assumir valores elevados, saturando o 


conversor DA ou excedendo os limites dos atuadores presentes no sistema. 
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m CONTROLE DEADBEAT 
es EXEMPLO 


* Para o sistema G(s) mostrado abaixo, obtenha um controlador deadbeat 


para: 
G(s) == 
° (+ DG + 10) 
ea)7T,=0,02s 
° b) T. = 0,1 s 
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m CONTROLE DEADBEAT 


EXEMPLO 


Solucáo para Ts = 0,02s: 
%Sistema analógico 

g = ti(1,poly([-1 -10])); 
zpk(g) 


% Taxa de Amostragem 
Ts = 0.02; 


% Sistema dicreto 
gd = c2d(g, Ts); 
zpk(gd) 
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m CONTROLE DEADBEAT 
ee EXEMPLO 


% Planta discreta: z^-k (К --> Número de polos - número de 
zeros) 

e gclz = tf(1,[1 0], Ts); 

e Zpk(gclz) 


e % Obtencáo do controlador direto 


e С = 1/gd*(gclz/(1-gclz)); 
° zpk(c) 
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m CONTROLE DEADBEAT 
we EXEMPLO 


* Resposta ao Degrau: 


Step Response 
I 


0.8 - 


Amplitude 


0.6 - 


0.4 - 


0.2 - 


0 
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 


Time (sec) 
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0.4 


m CONTROLE DEADBEAT 
we EXEMPLO 


4* Sinal de controle: 


x 10 Step Response 
T T 


Amplitude 
o 
zl 
d 
J 


| 
0 0.2 0.4 0.6 0.8 1 1.2 1.4 


Time (sec) 
Eng. de Controle e Automação – Rejane de Barros Araüjo/André Ferreira 


m CONTROLE DEADBEAT 
ilt EXEMPLO 


e %Solucáo para Ts = 0,15: 
• % Sistema analógico 

° g-tf(1;poly([-1 -10])); 

° zpk(g) 


e % Taxa de Amostragem 
e [5 = 0.1; 


e % Sistema dicreto 
e gd = c2d(g, Ts); 


° Zpk(gd) 
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m CONTROLE DEADBEAT 
ee EXEMPLO 


% Planta discreta: z^-k (k --> Nº de polos — Nº de zeros) 
gclz = tf(1,[1 0], Ts); 
e zpk(gclz) 


e % Obtenção do controlador direto 


c= 1/gd*(gclz/(1-gclz)); 
° zpk(c) 
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E CONTROLE DEADBEAT 
= EXEMPLO 


* Resposta ao Degrau: 


Step Response 
| | I 


0.8 - 


Amplitude 


0.4 - 


0.2 - 


Time (sec) 
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m CONTROLE DEADBEAT 
we EXEMPLO 


> Sinal de controle 


Step Response 
| 


300 r 


| 1 


200 | 


100 + 


Amplitude 


-100 - 


-200 - 


-300 | | | | 
0 0.5 1 1.5 2 


Time (sec) 
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É. ExERCÍCIOS 


** Projete o controlador deadbeat para o sistema anterior de tal forma que 


o sinal de controle esteja na faixa —15 < u(k) < 15. 


e Solução: 


e Usar T; > 1,25 s. 
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É. CONTROLADOR PID CONTÍNUO 


* A estrutura ideal (equação matemática clássica) do controlador PID tem 


pouca utilização em aplicações práticas. 


** As estruturas de implementação do PID podem ser classificadas: 
** pelo tipo de interação entre as ações P, I, D. 


** pela distribuição das ações P e D na malha de controle. 


** PID ideal 


U (s) 
E(s) 


1 


= С (5) = «[ m Ts 
Г. 
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E: CONFIGURAÇÃO E VARIAÇÃO DO PID 
т PID IDEAL 


proporcional integral 


e(t) u(t) 


derivativo 


u(t) — K, («co + - j e(t)dt + T, 2) 
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E: CONFIGURAÇÃO E VARIAÇÃO DO PID 
e PID SERIE 


= 
derivativo 


e(t) e 


a u(t) 


proporcional 
integral 
U(s x 1 ^ 
) 6.) = К.| 1+5 (1 +7,5) 
к) Tis 
i 
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E: CONFIGURAÇÃO E VARIAÇÃO DO PID 
т PID PARALELO 


proporcional 


integral 


u(t) 


derivativo 
U(s) 
E(s) 


Ki 
= Ge(S) = K< 1 Kas 
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E: CONFIGURAÇÃO E VARIAÇÃO DO PID 
т PI+D 


r(t) e(t) 
E integral proporcional 


u(t) j y(t) 
Эс ИСЭ 


derivativo 
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E: CONFIGURAÇÃO E VARIAÇÃO DO PID 
т I-PD 


integral proporcional 


0 mna ш 
Š 


r(t) m e(t) 


| 
ө 


derivativo 
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E: CONFIGURAÇÃO E VARIAÇÃO DO PID 
т OBSERVAÇÕES 


« Os parâmetros de projeto dos controladores são diferentes e as funções 


de transferência não são sempre equivalentes. 


* No PID série, os zeros do controlador somente podem ser reais, 
enquanto que no PID paralelo e ideal, podem ser reais ou complexos 


conjugados. 


* A caracteristica de polos complexos conjugados pode ser útil quando o 


sistema a ser controlado tem polos complexos conjugados dominantes. 
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E: CONFIGURAÇÃO E VARIAÇÃO DO PID 
ae OBSERVAÇÕES 


% О PID série é denominado de interativo, pois o ajuste de parcela 


derivativa afeta a parcela integral. 


* A estrutura paralela é a mais real, porém seus parâmetros não têm 


interpretação física. 


* Quando se trabalha com conceito de polos e zeros, a estrutura série tem 
a vantagem de que cada zero do PID está relacionado com um 


parâmetro do controlador (T; ou T4). 
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E: CONFIGURAÇÃO E VARIAÇÃO DO PID 
ae OBSERVAÇÕES 


** Na estrutura I+PD, o “kick” da parcela proporcional é subtraindo. 


* Exemplo : PID discreto 


* Exemplo : simulação 
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E: CONFIGURAÇÃO E VARIAÇÃO DO PID 
= PID INDUSTRIAL 


PID série — 


U(s) |. 1 || Tas +1 
E(s) Ke E a E + 1 


Бю аз — “кые г. 
нек Es) “ Tis aT4s-c1 
PID ideal U(s) 


com filtro ” ES). 
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Agenda 


“ Representação de sistemas discretos no Espaço de Estados 


< Observabilidade e Controlabilidade 


* Alocação de polos utilizando realimentação de estados vt, 


м. 


< Observador de estados Y 
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@ 
E SISTEMAS DIGITAIS NO ESPAÇO DE ESTADOS 


4 Equações no espaço de estados são compostas de três tipos de 


variáveis utilizadas para modelar a dinâmica de um sistema: 
Y Entrada; 
“ Saída; 


Y Estado. 


“ Representação não é única. 


€ Número de estados é o mesmo em qualquer representação utilizada. 
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E SISTEMAS DIGITAIS NO ESPAÇO DE ESTADOS 


* Para sistemas digitais variantes no tempo (linear ои não linear) 


Y Equação de estados 


x(k +1) = flx(k),u(k), К] 


Y Equação de saída 


y(k) = g|x(k), u(k), К] 
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SISTEMAS DIGITAIS NO ESPACO DE ESTADOS 


* Para sistemas lineares variantes no tempo, as equações de estado e de 


saída podem ser simplificadas para: 


x(k + 1) = G(k)x(k) + H (k)u(k) 


y(k) = C(k)x(k) + D(k)u(k) 


x(k) — vetor de estados (ordem n) 
y (k) – vetor de saídas (ordem m) 


u(k) — vetor de entradas (ordem r) 


v 

v 

v 

Y G(k)- matriz de estados (n x n) 

Y H(k) — matriz de entradas (n x r) 

Y C(k) — matriz de saídas (m x n) 

Y D(k)- matriz de transmissáo direta (m x r) 
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E SISTEMAS DIGITAIS NO ESPAÇO DE ESTADOS 


% Caso o termo k não apareça nas matrizes, supõe-se que as mesmas 


são constantes (sistema linear invariante no tempo): 


x(k + 1) = Gx(k) + Hu(k) 


y(k) = Cx(k) + Du(k) 
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- SISTEMAS DiGITAIS X SISTEMAS CONTÍNUOS NO 


ID» ESPACO DE ESTADOS 


x(t) = fix), u(6) t] x(k + 1) = fIx(k), u(k), К] 
y(t) = glx(t), u(t), t] y(k) = glx(k), u(k), К] 


х@) = А@)х@) + B(t)u(t) x(k + 1) = G(k)x(k) + H(k)u(k) 
y(t) = C(t)x(t) + D(t)u(t) y(k) = C(k)x(k) + D(k)u(k) 


x(t) = Ax(t) + Bu(t) x(k + 1) = Gx(k) + Hu(k) 
y(t) = Cx(t) + Du(t) y(k) = Cx(k) + Du(k) 


Observacáo: As matrizes C e D sáo as mesmas para os dois casos! 
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E REPRESENTAÇÃO EM DIAGRAMA DE BLOCOS 


(b) 
(a) Digital (b) Contínuo 
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E FORMAS CANÓNICAS PARA EQUACOES NO 


m» ESPACO DE ESTADOS DISCRETAS 


Considere um sistema discreto representado por: 


y(k) + a4y(k — 1) + azy(k — 2) + + any(k — n) 
= bgu(k) + byu(k-1)+--+bpu(k—n) 


* 


% u(k) e y(k) são os valores de entrada e de saída do sistema no k-ésimo 


instante de amostragem. 


% Colocando na forma de função de transferência pulsada: 


Y(z) _ bo + bz l +- + 27" Y(z) _ boz” + bz" 1 L... + bn 
U(z) 1+az™t+ + ал27" U(z) zn +azn-1 +... + QA 
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E FORMAS CANÓNICAS PARA EQUAÇÕES NO 


m ESPACO DE ESTADOS DISCRETAS 


% Das infinitas representações по espaço de estados, destacam-se: 


$* 1) Forma canónica de controlador 
* 2) Forma canónica de observador 
$* 3) Forma canónica diagonal 


« 4) Forma canónica de Jordan 
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- FORMA CANÓNICA DO CONTROLADOR 


x(k + 1) 0 1 0 е $) xı (Kk) 0 
x(k + 1) 0 0 1 + 0 || x2(k) 0 
: =) É .. +lilu(k) 
Xn-1(k + 1) 0 0 O >= 1 xaa(k)| |0 
Xn(k + 1) “ln Uns  —0gp-2 `7 Uil x,(k) 1 
xy (k) 
y(k) = [b; i bzi c ba] 20) + ue 
Xy (k) 
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- FORMA CANÓNICA DO CONTROLADOR 


+ Invertendo a ordem dos estados, isto é, fazendo: 


£ (k) 0.0 .. 0 1][%xƏi (k) 
ewj jo 0 = 1 о 


eco) li 0 = о ою 
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FORMA CANÓNICA DO CONTROLADOR 


#(К+1)] pa -a | ааа AA] p 
X5(k + 1) 1 0 -. 0 о || 22k) | |0 
X3(k+1)|=]| 0 1 B 0 O || Xa(k) |+ |O[u(k) 
„(К+ 1) 0 0 ut 1 04/2 ССК) 0 
£k) 
I „(К 
y() = [ii bzi oc ba] O| + bouo 
£, (k) 
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FORMA CANÓNICA DO OBSERVADOR 


x(k + 1) 00-00 на, || x(k) b, — a, bg 
x(k + 1) 1 0 didis 0 0 Un? x(k) b,.4 — Qn-1 bo 
: =j: : . : : + : u (Kk) 

Xn-ı(k T 1) 0 0-1 0 — o Xp-1(Kk) b, — d bo 
Xalk + 1) 0.0 --- 0 1 -a Xn(k) b, — а, б, 

x (K) 

i x(k) 

y(k)=[0 0 +++» 0 JJ : + byul(k) 
Xn-1(k) 
x,(k) 
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FORMA CANÓNICA DO OBSERVADOR 


4 Fazendo a mesma inversão de ordem dos estados mostrada 


anteriormente: 
X (k + 1) — a 1 0 0 0 X (k) b, — 41 bo 
Xo(k + 1) 022 0 1 --- 0 0 (К) b, = a bo 
O E o : Р |+ ll) 
Xn-1(k T 1) —4n-1 0.0 -. 0 1 Xn-1(Kk) b,4 Qn-1 bo 
£,(k + 1) sa 00 o 0] £,(k) b, — an bo 
X (Kk) 
(К) 
у(к)=[1 0-- 00] : |+bou(k) 
Xi - (К) 
Xn(k) 
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@ A 
- FORMA CANÓNICA DIAGONAL 


% Se os polos da função de transferência pulsada de um sistema forem 


todos distintos: 


х1(К + 1) p 0 +. O]pa(k) 1 
dp СУШ 
adn 0 о .. р, xp (k) 1 
х1(К) 
у(к) = [а o 7 6| 209 + ul) 
Xp (k) 


“ onde р; são os polos do sistema e c; são os resíduos associados a cada 


polo. 
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- FORMA CANÓNICA DE JORDAN 


% Se a função de transferência pulsada do sistema apresenta m polos 


repetidos, entáo pode-se colocar o sistema na forma 


xi(k + 1) p 1 0 --- 010 O | | x(k) 0 
xk + 1) 0 Pi 1 --- 0 0 0 xo(K) 0 
x4K*1]| [00 0--pi0 xml) | | 1 


tL jJu(k) 


Xm+i(k + 1) 00 0..- 0 Pm+1 cC 0 Xm+1[k) 1 
x(k + 1) 000-010 -piix | li 
x (k) 
у(к)=[с o D+ bulto) 
x, (k) 
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E 
E EXEMPLOS 


HEP — | y 


Y(z) |. z+1 
U(z) 22 + 1.32 + 0.4 


< Forma Canônica do Controlador: 
xi(k+1)] po 1 ][x(k)] , [0 
bs "d Lia mis P + [ije 


›®=[1 up 


4 Forma Canónica de Observador: 


a + | ú Ë el 2102) | 


xk + DE la Ml Wol ¡Juas 


y()=1[0 1] Е 
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E 
E EXEMPLOS 


NEP o 


Y(z) |. z+1 
U(z) 22 + 1.32 + 0.4 


€ Expandindo em frações parciais: 


Y(z) _ 5/3 —2/3 
U(z) RE aa 


“ Forma Canônica Diagonal: 


Berol e^ -ollel * ho 


xo-E Ia 


Eng. de Controle e Automacáo — André Ferreira/Rejane de Barros Araüjo 


TRANSFORMACOES LINEARES 


x(k + 1) = Gx(k) + Hu(k) 
y(k) = Cx(k) + Du(k) 


* 


Ф 


* Utilizando uma matriz Р não singular: 
x(k) = PX(k) — PX(k + 1) = GPX(k) + Hu(k) 


* Pré-multiplicando cada lado da equação por Рт: 


£(k + 1) = P l1GP£(k)+ P ! Hu(k) 


* Definindo: P -1GP = 6, P^OH-H 
~ Obtém-se: 


&(k + 1) = GÉ(k) + Hu(k) 
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- TRANSFORMACOES LINEARES 


« Utilizando a mesma transformação de estados na equação de saída: 


y(k) = CPX(k) + Du(k) 


4. Definindo: 


O) 


CP=Ê D= 


4$ Chega-se a seguinte relação: 
y(k) = ERC) + Du(k) 
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š TRANSFORMAÇÕES LINEARES 


“ Desta forma: 

x(k + 1) = Gx(k) + Hu(k) x(k) = РК) — “(k+1)= GRC) + Hu(k) 
y(k) = Cx(k) + Du(k) — y (lO) = 2k) + Dulk) 
Caso se deseje diagonalizar a matriz de estados, pode-se escolher uma 

matriz P, de tal forma que: 
Р-Р = matriz diagonal 
Nos casos onde isso nào seja possível, pode-se usar a forma canónica 
de Jordan: 
P^ !GP = J = matriz de Jordan 


% Ver Apéndice A (Аб) - Ogata (Discrete Time Control Systems) 
“ Мега função "jordan" no MATLAB 
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n" TRANSFORMAÇÕES LINEARES 


ВЕСУ EXEMPLO 


* Dado o sistema representado no espaço de estados, transforme о 
sistema em um novo conjunto de variáveis de estado, £(k + 1), onde as 
novas variáveis de estado se relacionam com as variáveis de estado 


originais, x(k + 1), através da seguinte relacáo: 


x4 (k + 1) 0 1 0 ] [xi CK) 0 
x (k +1)] =| 0 0 1 | x2(k)|+|0|u(k) 
x3(k + 1) —2 —5 —71|x3(k) 1 
£4A(k + 1) = 2x4 (k) 
x1 (k) 
у(к) = [1 0 0] [x2(k) Z2(k + 1) = 3x4 (k) + 2x2(k) 
x3(k) 


X3(k + 1) = x,(k) + 4x5(k) + 5x3(k) 


Eng. de Controle e Automação — André Ferreira/Rejane de Barros Araújo 


n" TRANSFORMAÇÕES LINEARES 
EXEMPLO - SOLUÇÃO 


2 0 0 
X(k + 1) = з 2 ] x(k) = P^!x(k) 
1 4 5 


2 0 0J[O 1 0 0.5 0 0 =15 1 0 
ғо = | 2 JE 0 1 ||-075 0.5 0 |=|-1.25 07 04 


—2 —5 =L 05 —0.4 0.2 —25 04 —62 


atu 


0.5 0 0 
CP=[1 0 ЛЕ 0.5 0 |= [0.5 0 0] 
05 —04 02 


P ` 
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- TRANSFORMACOES LINEARES 
EE» EXEMPLO - SOLUCAO 
< O sistema transformado é: 


&(k + 1) = GÉ(k) + Hu(k) 
y(k) = C£(k) + Du(k) 


£ (k + 1) —15 1 0 7[%1(*%) 0 
Zo (k + D| = E 0.7 o4 | &2(к) | + |O| u(k) 
£, (k + 1) —25 04 -—6.21|£3(k) 5 


£ (k) 
y(k)= [0.5 0 0] |X>(k) 
Хз(К) 
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E SOLUÇÃO DA EQUAÇÃO DE ESTADOS PARA SISTEMAS 
ED DISCRETOS LINEARES INVARIANTES NO TEMPO 


x(k + 1) 2 Gx(k) + Hu(k) 
y(k) = Cx(k) + Du(k) 


4 Pode ser resolvida através de métodos iterativos. 


4 Para К inteiro 


x(1) = Gx(0) + Hu(0) 
x(2) = Gx(1) + Hu(1) = G*x(0) + GHu(0) + Hu(1) 
x(3) = Gx(2) + Hu(2) = G?x(0) + G^Hu(0) + GHu(1) + Hu(2) 
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E: _ SOLUÇÃO DA EQUAÇÃO DE ESTADOS PARA SISTEMAS 
me» DISCRETOS LINEARES INVARIANTES NO TEMPO 


4 Assim: 


x(k) = G*x(0) + X5-d G*-^1Hu(j), К = 1,2,3,.... 


% Substituindo x(k) na equação de saída, tem-se: 


y(k) = CG*x(0) + Су G*-1Hu(J) + Du(k) 
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E SOLUÇÃO DA EQUAÇÃO DE ESTADOS PARA SISTEMAS 


кз DISCRETOS LINEARES INVARIANTES NO ТЕМРО 


% É possível escrever a solução para a equação de estados homogênea 


x(k + 1) = Gx(k) como sendo: 


x(k) = V(k)x(0) 


% onde Y(k) é uma matriz nxn única que satisfaz a relação: 


W(k +1) = GU(k) W(0) = 1 
W(k) = GF 
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E SOLUÇÃO DA EQUAÇÃO DE ESTADOS PARA SISTEMAS 


IES DISCRETOS LINEARES INVARIANTES NO TEMPO 


4 A matriz (k) é chamada de matriz de transição de estados ou matriz 


fundamental. 


% Esta matriz contém todas as informações sobre o movimento livre do 


sistema definido pela equação x(k + 1) = Gx(k). 


% Desta forma e usando a matriz 'P(k), pode-se reescrever a equação 


abaixo: 
x(k) = G*x(0) + Djzo GHU), k = 1,2,3,... 


k—1 | k—1 


x(k) = W(k)x(0) + > V(k — j — DHu(j) = Y(k)x(0) + > W()Hulk — j — 1) 
j=0 j=0 
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E SOLUÇÃO DA EQUAÇÃO DE ESTADOS PARA SISTEMAS 


EE DISCRETOS LINEARES INVARIANTES NO TEMPO 


€ Fazendo a mesma substituição na equação 


y(k) = CGEx(0) + C X5-d G 7 Hu(j) + Du(k) 


| 


y(k) = СФ(К)х(0) + СУ) Y(k – j – 1)HuG) + Du(k) 


y(k) = СФ(К)х(0) + С Уо Wü)Hu(k — j — 1) + Du(k) 
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E SOLUÇÃO DA EQUAÇÃO DE ESTADOS PARA SISTEMAS DISCRETOS 


EER LINEARES INVARIANTES NO TEMPO USANDO TRANSFORMADA Z 


x(k + 1) = Gx(k) + Hu(k) 
% Aplicando a transformada Zem ambos os lados: 


ZX(z) — zx(0) = GX(z) + HU(z) 


* com X(z) = z[x(k)] e U(z) = z[u(k)] 
* Então: (zl — G)X(z) = zx(0) + HU(z) 
% Desta forma, tem-se: X(z) = (zI — G) !zx(0) + (zl — G) !HU(z) 


* Aplicando-se a transformada inversa em ambos os lados, obtém-se: 


x(k)-2z [GI —G) !z]x(0) +z ! [GI — G) 1 HU(z)] 
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E SOLUÇÃO DA EQUAÇÃO DE ESTADOS PARA SISTEMAS DISCRETOS 


ES LINEARES INVARIANTES NO TEMPO USANDO TRANSFORMADA Z 


% Comparando as equações abaixo, conclui-se que: 
k—1 


x(k) = G*x(0) + > GF Hu(j) 


J=0 


x(k) = z ![GI — G) 1z|x(0) + z ! [GI — G) 1 HU(z)] 


* 1) Gk = z ![GI — G) 12] 
+ 2) Xio Gi Hu) = z ![GI — G) 1 HU(z)] 
% OBS.: A utilização desse método requer a inversão de algumas matrizes, 


o que pode ser feito analiticamente ou com auxílio de um computador. 
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E 
E EXEMPLOS 


Do —  ... y 


% Obtenha a matriz de transição de estados para o sistema abaixo: 


x(k + 1) = Gx(k) + Hu(k) 


y(k) — Cx(k) 
d Los 4 m H c-h o] 


* A seguir, obtenha o estado x(k) e a saída y(k), para uma entrada 
u(k) = 1, parak = 0,1,2,... 


Considere que: 


0 
AS Ba E E m 
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E EXEMPLOS - SOLUÇÃO 


% Solução: 
W(k) = Gk = z-1[(zI — G) 12] 


*$* Assim: 
" -1 — Z —1 
(zI — G) ie 23:1 
z+1 1 
_ | (z + 0.2)(z + 0.8) (z + 0.2)(2 + 0.8) 
—0.16 E 


(z + 0.2)z + 0.8 (2 + 0.2)(z + 0.8) 


4 =} 5 5 
3 de a — NE um 3 
z+0.2 z+0.8 2302 + + 0.8 
mi — 0.8 0.8 "S 4 
3 „өтөгү ENS IP Бас — s SN 
z + 0.2 + 0.8 $0105" z+ 0.8 
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EXEMPLOS - SOLUÇÃO 


% Solução: 


* A matriz de transição de estados V (k) vale: 


W(k) = С^ = z-1[(zI — G) `1z] 


4 z _l z 322 _ 5 Z 
= 9-1 32+02 3z+0.8 3z2+0.2 32 + 0.8 
0.8 2 0.8 2 1 2 4 2 


32102 32+08 3z*02 3z408 


4(—0.2)* — 1(—0.8)“ s(-0.2)' — $(-0.8)* 


-98(—0.2)* + 28(—0,8)* —1(—0.2)* + 2(—0.8)* 


Eng. de Controle e Automação — André Ferreira/Rejane de Barros Araújo 


Ё. EXEMPLOS - SOLUÇÃO 


“ Usando o MATLAB: 
% >> G = [0 1;-.16 -1] 
>> Н = [1 1] 
% >> С = [1 0] 
% >> syms z К 


* 


“ >> Iztrans(inv(z*eye(2)-G)*z) 


Eng. de Conirole e Automação — André Ferreira/Rejane de Barros Araújo 


g РЕ 
- EXEMPLOS - SOLUCAO 


% Cálculo de x(k): 
z|x(k)] = X(z) = GI — G) !zx(0) + (zl — G) ! HU(z) 
= (zI — G) ![zx(0) + HU(z)] 


“ Рага uma entrada degrau unitário, tem-se: 


1 
U(z) = 1— = 


% O segundo termo da multiplicação 


—z2 + 27 


zx(0) + HU(z) = | ME EE zd 
— 1 
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E. 


4 Desta forma: 


EXEMPLOS - SOLUÇÃO 


X(z) = (zI — G) ![zx(0) + HU(z)] 


(22 + 2)z 
[162+ 0.2) (z + 0.8) (х — 1) 
(=z? + 1.84z)z 
(z + 0.2)(z + 0.8)(z — 1) 


€ Expandindo em frações parciais: 

17 

z + 0.2 

3.4 
6^ 

z + 0.2 


X(z) = 


9^ | 18 
zZ+08 z-1 
176, 7 

02 18* 


z+ 0.8 7—1 
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g РЕ 
- EXEMPLOS - SOLUCAO 


« O vetor de estados x(k) vale desta forma: 


17 „ 22 „25 
Y e + == 
5 1 ы 17.6 7 
> ок — 2E (0.8) + + 
— (0.2 — (—0.8)* +— 
* Easaída 
25 


-U Cozy 4I Cog 42 
у(Е) 2Cx(k) =[1 01|, ° 9 18 


= Lc 02* ro 
67 9^ ` 18 
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: _ : 
E MATRIZ FUNÇÃO DE TRANSFERÊNCIA PULSADA 


% Sistema discreto linear invariante no tempo de ordem n, com r entradas 


е т saídas: 
х(К + 1) = Gx(k) + Hu(k) 


y(k) = Cx(k) + Du(k) 


% Aplicando-se a transformada Z: 
zX(z) — zx(0) = GX(z) + HU(z) 
Y(z) = CX(z) + DU(z) 
€ Da definição de função de transferência, sabe-se que as condições 


iniciais são nulas. Então: 
X(z) = (zI — G) ! HU(z) 


Y(z) = [C(zl — G) 1H + D]JU(z) = F(z)U(z) 
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: _ : 
E MATRIZ FUNÇÃO DE TRANSFERÊNCIA PULSADA 


4 Assim: 


F(z) = C(zI =G) 1H + D 


* onde F(z) é a matriz função de transferência pulsada, de ordem mxr. 


adj(zI-G) 


À | _осү-1_ 
* Lembrando que: (zI — G) = EET 


C adj(zI-G)H 


4 Assim: F(z) = |21—С | 


+ р 


e Desta forma, ё fácil notar que os polos de F(z) são as raízes de 


4 


|21— G| = 0. De forma alternativa, pode-se usar os elementos a, da 
matriz б: zt Eq Zl аз eta a, = Ü 
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E OBTENCAO DAS MATRIZES DE ESTADO USANDO 


| 45 O MATLAB 


x = Ax + Ви —— x(k + 1) = Gx(k) + Hu(k) 


* Para obter as matrizes de estado discretas a partir do caso contínuo, 
utiliza-se a função “c2d”, onde T, é o período de amostragem em 


segundos: 


[G, H] = c2d (A, B, T.) 
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E 
E EXEMPLOS 


Do — 


% Obtenha as matrizes de estado discretas para o sistema abaixo, 


supondo um período de amostragem de 50 ms. 


| E pr he [s] a H н 


% Solução: 


* A=[0 1;-25 -4]; 

_ [0.9709 0.0448 
*% В = [0;1]; — 1—1.1212 0.7915 
«e [G,H] = c2d(A,B,0.05) 0.0012 


= К 
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E. 
a CONTROLABILIDADE 


% Um sistema é dito ser controlável se para qualquer condição 
inicial x(kg) existe uma sequência de controle u(k), k = ko, 
ko +1,...,kf — 1, tal que um estado final arbitrário x(k,) 


pode ser alcançado em um finito k,. 
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Ё. MATRIZ DE CONTROLABILIDADE 


x(k + 1) = Gx(k) + Hu(k) 


y(k) = Cx(k) + Du(k) 


€ Um sistema linear de ordem n invariante no tempo é completamente 
controlável se sua matriz de controlabilidade С tem rank completo (igual 
an). 


% C=[H GH .. g«-Dg] 
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- MATRIZ DE CONTROLABILIDADE 


46, EXEMPLO 


4 Determine a controlabilidade do sistema abaixo: 


х1(К +1) -2 —2 0 1[х:(К) 1 0 
x (k + 1)| = | 0 0 1 | x5(k) | + |0 iuo 
x3(k + 1) O —0.4 —0.51|x3(k) 1 1 
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- MATRIZ DE CONTROLABILIDADE 


К, EXEMPLO 


“ Solução: 
% >> G = [-2 -2 0;0 0 1;0 -0.4 -0.5] 
% >> Н = [1 0;0 1;1 1] 
# >> aux = [H G*H G^2*H] 


** >> rank(aux) 


*• No MATLAB, pode-se também utilizar a função “ctrl”: 
* >> aux = (ctrb(G,H)) 
** >> rank(aux) 


O rank da matriz é completo (3) Ə sistema é completamente controlável 
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- 
а OBSERVABILIDADE 


« Um sistema é dito ser observável se qualquer estado inicial 
x(ko) pode ser estimado a partir da sequência de controle 
u(k), k = ko, ko + 1,..., К, — 1 e das medidas y(k), k = ko, 
ko lot 
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Ё. MATRIZ DE OBSERVABILIDADE 


x(k + 1) = Gx(k) + Hu(k) 


y(k) = Cx(k) + Du(k) 


€ Um sistema linear de ordem n invariante no tempo é completamente 


observável se sua matriz de observabilidade O tem rank completo. 
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- MATRIZ DE OBSERVABILIDADE 


BED EXEMPLO 


< Dadas as matrizes de estado e de saída, determine a observabilidade do 


sistema abaixo: 
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- MATRIZ DE OBSERVABILIDADE 


шэ EXEMPLO 


% Solução: 
4$ >> G = [0 1 0;0 0 1;0 -3 4] 
* >> С = [0 0 1] 
>> aux = [C;C'G;C*G^2] 


** >> rank(aux) 


* No MATLAB, pode-se também utilizar a função “obsv”: 
>> aux = (obsv(G,C)) 
% >> rank(aux) 


O rank da matriz não é completo (3) Ə sistema apresenta um modo não observável 
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ü ~ 
F. REALIMENTACAO DE ESTADOS 


Do —.— A  h 


% Consiste em utilizar o vetor de estados para calcular a ação de controle 


do sistema. 


vík)  u(K) x(+1) x(k) v(k) 


* Para um sistema linear e invariante no tempo de ordem n е com m 
entradas, a matriz de realimentação de estados K é constante e possui 


ordem m x n. 
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F. REALIMENTACAO DE ESTADOS 


NEP — ll L| Ó 


v(k) uw% x(k+1) x(k) У(А) 
+ 


** A partir do diagrama de blocos acima, pode-se escrever as seguintes 


equações: 


x(k + 1) = Gx(k) + Hulk) 


y(k) = Cx(k) 
u(k) = —Kx(k) + v(k) 
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ü ~ 
F. REALIMENTACAO DE ESTADOS 


шу oo 


% Manipulando essas equações, pode-se definir a matriz de estados em 
malha fechada: 
x(k + 1) = Gx(k) + H|-Kx(k) + v(k)] 


= [G — HK |x(k) + Hv(k) 


Ga = G — HK 


* Assim, o sistema em malha fechada pode ser representado por: 


x(k + 1) = G¿x(k) + Hv(k) 


y(k) = Cx(k) 
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F. REALIMENTACAO DAS SAÍDAS 


< Para muitos sistemas físicos, medir todas as variáveis de estados é 


economicamente inviável ou até mesmo impossível de ser 


implementada. 


% Uma alternativa é se fazer a realimentação dos sinais de saída, como 


mostrado a seguir: 


v(k) u(k) x(k*1) x(K) v(k) 
— G 
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F. REALIMENTACAO DAS SAÍDAS 


NEP — 3101 h3Jho— 


v(k)  u(k) x(k+1) x(k) v(K) 
D + 


% O controle usando a realimentacáo das saídas vale: 
u(k) = —Kyy(k) + v(k) = —K,Cx(k) + v(k) 
« Assim: 


x(k + 1) = Gx(k) + H|-K,Cx(k) + v(k)] = [G — HK,C|x(k) + Hv(k) 


€ Com a matriz de estados neste caso sendo dada por: 


Су = G — HKyC 
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- REALIMENTAÇÃO DE ESTADOS X 


EEE» REALIMENTAÇÃO DAS SAÍDAS 


* O controle usando realimentação das saídas é inferior ao controle 
através da realimentação de estados: 
4$ A matriz C restringe a escolha da dinâmica desejada em malha 


fechada. 


* Realimentacáo de estados é um caso particular da realimentação das 


saídas, quando a matriz С é uma matriz identidade. 
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E _ CONTROLABILIDADE X ESTABILIZABILIDADE 


ED OBSERVABILIDADE X DETECTABILIDADE 


% 1) Sistema Controlável > matriz de controlabilidade apresenta rank 


completo. 


4$ 2)Gistema Estabilizáve)> todos os seus modos nào controláveis são 


assintoticamente estáveis. 


4$ 3)CSistema ObserváveD> matriz de observabilidade apresenta rank 


completo. 


“ 4) Sistema Detectável > todos os seus modos não observáveis decaem 


assintoticamente a zero. 
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E _ CONTROLABILIDADE X ESTABILIZABILIDADE 


IES OBSERVABILIDADE X DETECTABILIDADE 


4 CONCLUSÃO: 
% a) Se um sistema tem modos não controláveis, então esses modos 
não sofrem influência do sinal de controle aplicado e podem (ou não), 


decair a zero assintoticamente. 


% b) Se um sistema é não observável, deseja-se que todos os modos 


não observáveis sejam estáveis. 


Eng. de Conirole e Automação — André Ferreira/Rejane de Barros Araújo 


E ALOCAÇÃO (IMPOSIÇÃO) DE POLOS 


“ Utilizando a realimentacáo de estados ou realimentação das saídas, os 
polos do sistema podem ser escolhidos para que o sistema em malha 


fechada apresente as especificações de desempenho desejadas. 


“* Objetivo: 
* Escolher a matriz de ganhos K (ou К,) que faz os polos do sistema 


serem dados pelo conjunto arbitrário: 


(Ai = 1, aon} 
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E ALOCACAO DE POLOS USANDO 


BED REALIMENTACAO DE ESTADOS 


x(k + 1) 2 Gx(k) + Hu(k) 
y(k) = Cx(k) 
“ Definição: 
4 Se o par (GH) é controlável, então existe uma matriz de 


realimentação de ganhos K que é capaz de impor arbitrariamente 


os polos do sistema em qualquer conjunto (A, i = 1,...,n). 


4 Se o par (G,H) é apenas estabilizável, então todos os modos 


controláveis podem ser alocados arbitrariamente. 
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E ALOCACAO DE POLOS USANDO 


Э REALIMENTACAO DE ESTADOS 


* Procedimento geral: 


e, 


* 1) Obtenha o polinômio característico desejado: 


д.0) = | [a -20 
і=1 


% 


* 2) Calcule o polinômio característico em malha fechada usando а 


expressão: 


det(AL, — (G — HK)) 


e 


+ 3) Iguale os coeficientes das duas equações anteriores e obtenha n 


equações cujas soluções são os elementos que formam a matriz K. 
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E EXEMPLO 


Do — 53 1 hh. 


* Рага o par (G, H) mostrado abaixo, determine a matriz К que faz com 


que os polos em malha fechada estejam localizados em 0,3 + j0,2: 
_ (0 1 _ [0 
m P 4 im | | 


% Solução: 
* a) Sistema é controlável? 
% >> aux=ctrb(G,H) 
% >> rank(aux) 


4 ans=? €————— Sistema controlável 
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E EXEMPLO 


шо hA. 3 


% Solução: 


« b) Polinómio característico desejado: 


A," (A) = (A — 0.3 — j0.2)(A — 0.3 + j0.2) = 42 — 0.62 + 0.13 


** c) Matriz de estados em malha fechada: 
= „yr .[0 1] [0 _ | 0 1 | 
кеа zi 4 MIL k2] = 3—-k 4—К, 
« d) Polinómio característico em malha fechada: 


1 


det(Al, — (G — НК} = det a kı) A— (4 — k2) 


|- 4 @ — k1- G - ka) 
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E EXEMPLO 


BED — 


“ Solução: 


« e) Igualando-se os coeficientes das equações anteriores: 


1.4—k; 206 > k=34 


———» КТ = [3, | 
2.-3+k,=0.13 > k = 3.13 uS ud 
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- ALOCACAO DE POLOS USANDO O MATLAB 


“ Para sistemas de ordem elevada, o método anterior pode ser muito 


trabalhoso. 


% O MATLAB possui duas funções que fazem o cálculo da matriz K, tanto 
para sistemas contínuos quanto para sistemas discretos: 


+ Acker; Place 


% A função “acker” baseia-se na fórmula de Ackermann e é satisfatória para 


sistemas SISO de ordem « que 10, podendo lidar com raízes repetidas. 


% A função "place" é melhor para sistemas de ordem mais elevada e pode 


lidar com sistemas MIMO, embora náo funcione para raízes repetidas. 
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Ei. ALOCAÇÃO DE POLOS USANDO O MATLAB 


ВЕСУ EXEMPLO 


% Determine a matriz de ganhos K para o sistema abaixo, sabendo que 


deseja-se que os polos em malha fechada estejam localizados em (0.1, 


0.4 + j0.4). 
01 0 01 0.01 

G=|0 0.5 од Н = 0 | 

02 0 04 0.005 


“ Solução: 
# >> G = [0.1 0 0.1;0 0.5 0.2;0.2 0 0.4] 
% >> Н = [0.01 0 0.005] 
% >> Pd = [0.1 0.4+j*0.4 0.4-j*0.4] 
% >> K = acker(G,H,Pd) 
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E ESCOLHA DOS AUTOVALORES DE 


BED MALHA FECHADA 


% A escolha da localização dos polos (autovalores) do sistema em malha 
fechada está diretamente ligada à resposta transitória desejada para o 


sistema. 


“ Deve-se levar em conta que polos associados com modos mais rápidos 
apresentaráo valores elevados na matriz de realimentacáo de estados: 
% Gera um grande esforço de controle; 
* Degradação no desempenho devido a elementos não-lineares. 


*Ех.: saturação nos atuadores e no conversor AD. 


4 Se todos os polos desejados forem alocados em zero, o controle será do 


tipo deadbeat. 
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- ESCOLHA DOS AUTOVALORES DE MALHA FECHADA 


EDD EXEMPLO 


* Para um sistema representado pelas matrizes no espaço de estados 


discreto G e H abaixo, com T — 0,01 s, 


1.0 0.1 0.0 1.622x1078 
G = É 0.9995 0.0095 Н = |4.821x10 * 
00 —0.0947 0.8954 9.468х107? 


« Determine a matriz de estados рага alocar os polos em: 


1.£0.1,0.4 + j0.4) 


3. (0, 0,0) (deadbeat control) 
2.{0.4,0.6 + j0.33] 


% Obtenha a resposta dos З casos para um estado inicial: 


1 
au) = B 
1 
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- ESCOLHA DOS AUTOVALORES DE MALHA FECHADA 


NE» EXEMPLO 


“ Solução: 
# >> G = [1 0.1 0;0 0.9995 0.0095;0 -0.0947 0.8954] 
* >> Н = [1.622е-6 4.821e-4 9.468e-2] 
% >> P1 = [0.1 0.4+j*0.4 0.4-j*0.4] 
# >> P2 = [0.4 0.64j*0.33 0.6-j*0.33] 
% >> РЗ = [0 0 0] 


* Matriz de ganho de realimentacáo рага cada caso: 
#>> К1 = acker(G,H,P1) 
“>> K2 = acker(G,H,P2) 
«>> КЗ = acker(G,H,P3) 
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- ESCOLHA DOS AUTOVALORES DE MALHA FECHADA 


| omg EXEMPLO 


“ Solução: 
« Sistema em malha fechada para o caso 1: 
> [s = 0.01; 
% >> 5у51 = ss(G-H'K1,H,-K1,0,Ts) 


* Resposta para x(0) do caso 1: 
%>>xX0 = [1 1 1] 
“>> [y1,t1,x1] = initial(sys1,x0,0.2); 
>> subplot(221),stairs(t1,y1),xlabel("Tempo (s)),ylabel('Caso 1 - u(k)') 
( (:,1)),xlabel("Tempo (s)),ylabel('Caso 1 - x1(k)') 
,Stairs(t1,x1(:,2)),xlabel('Tempo (s)),ylabel('Caso 1 - x2(k)') 
( (:,3)),xlabel("Tempo (s)),ylabel('Caso 1 - x3(k)') 
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>> subplot(222),stairs(t1,x1 


) 
(222) 
«>> subplot(223) 
(224) 


>> subplot(224),stairs(t1,x1 


E ESCOLHA DOS AUTOVALORES DE MALHA FECHADA 


ЕЕЗ EXEMPLO 


% Solução: 


* Resposta para x(0) do caso 1: 


4 


x 10 
1 1.5 
С 0.5 = 1 
5 x 
+ 0 T “os 
2 9 
% d 
1 -0.5 
0 0.05 0.1 0.15 0.2 0.05 0.1 0.15 0.2 
Tempo (s) Tempo (s) 
2 500 
— 0 E = 
X X 
Y e 0 
x x 
— -2 ú 
9 9 -500 
O -4 Ó 
-6 | | | -1000 : | | 
0 0.05 0.1 0.15 0.2 0 0.05 0.1 0.15 0.2 
Тетро (5) Tempo (5) 


Eng. de Controle e Automação — André Ferreira/Rejane de Barros Araújo 


- ESCOLHA DOS AUTOVALORES DE MALHA FECHADA 


NE EXEMPLO 


“ Solução: 
“ Sistema em malha fechada para o caso 2: 
$»» [s = 0.01; 
% >> sys2 = ss(G-H'K2,H,-K2,0,Ts) 


* Resposta para x(0) do caso 2: 
%>>xX0 = [1 1 1] 
“>> [y2,12,x2] = initial(sys2,x0,0.2); 
>> subplot(221),stairs(t2,y2),xlabel("Tempo (s)'),ylabel('Caso 2 - u(k)') 
( (:,1)),xlabel("Tempo (s)),ylabel('Caso 2 - x1(k)') 
,stairs(t2,x2(:,2)),xlabel("Tempo (s)),ylabel('Caso 2 - x2(k)') 
( (:,3)),xlabel("Tempo (s)),ylabel('Caso 2 - x3(k)') 
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>> subplot(222),stairs(t2,x2 


) 
(222) 
«>> subplot(223) 
(224) 


«>> subplot(224),stairs(t2,x2 


E ESCOLHA DOS AUTOVALORES DE MALHA FECHADA 


EED EXEMPLO 


% Solução: 
* Resposta para x(0) do caso 2: 


2000 1.5 


1000 


Caso 2 - u(k) 
3 
e 
o 


Caso 2 - x1(k) 
o 
q 


2000 9 a _— 
3000 0.5 
0 0.05 0.1 015 0.2 0.05 0.1 0.15 0.2 
Tempo (s) Tempo (s) 
1 100 == 
= 0 [ — g 0 A 
Y e 
x x 
ej + e “100 
О О 
102) [^] 
3 2 & -200 
-3 : | -300 . | 
0 0.05 0.1 0.15 0.2 0 0.05 0.1 0.15 0.2 
Тетро (5) Тетро (5) 
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- ESCOLHA DOS AUTOVALORES DE MALHA FECHADA 


NES EXEMPLO 


“ Solução: 
“ Sistema em malha fechada para o caso 2: 
> [s = 0.01; 
% >> sys3 = ss(G-H*K3,H,-K3,0, Ts) 


* Resposta para x(0) do caso 3: 
%>>xX0 = [1 1 1] 
>> [y3,13,x3] = initial(sys3,x0,0.2); 
>> subplot(221),stairs(t3,y3),xlabel("Tempo (s)),ylabel('Caso З - u(k)') 
( (:,1)),xlabel("Tempo (s)),ylabel('Caso З - x1(k)') 
,stairs(t3,x3(:,2)),xlabel("Tempo (s)),ylabel('Caso З - x2(k)') 
( (:,3)),xlabel("Tempo (s)),ylabel('Caso З - x3(k)') 
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«>> subplot(222),stairs(t3,x3 


) 
(222) 
«>> subplot(223) 
(224) 


«>> subplot(224),stairs(t3,x3 


E ESCOLHA DOS AUTOVALORES DE MALHA FECHADA 


ЕЕЗ EXEMPLO 


% Solução: 


* Resposta para x(0) do caso 3: 


4 


x 10 
3 15 
2 - 
= = 'r 
3 1 x 
со e 0.5 
o 0 o 
o o) 
Š " a o 
-2 | : | -0.5 : | : 
0 005 0.1 015 0.2 0 0.05 0.1 015 0.2 
Tempo (s) Tempo (s) 
2 2000 
ZO « 1000 
бү e 
x x 
e 2 РА 0 
о о 
3 3 
3 4 © -1000||| 
-6 | | | -2000 : | 
0 0.05 0.1 0.15 0.2 0 0.05 0.1 0.15 0.2 
Тетро (5) Тетро (5) 
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- ESCOLHA DOS AUTOVALORES DE MALHA FECHADA 


= EXEMPLO 


% Solução: 
* Conclusão: 
«Quando escolhe-se autovalores associados com modos rápidos: 


$* 1) Ganhos da matriz de realimentacáo serão elevados; 


* 2) Resposta transitória terá grandes oscilações. 
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Е. _ ALOCAÇÃO DE POLOS USANDO 


EE REALIMENTACAO DAS SAÍDAS 


“ А realimentação das saídas apresenta uma capacidade limitada de se 
atribuir os autovalores do sistema, comparado com o que é possível se 


fazer utilizando realimentacáo de estados. 


“ Geralmente não é possível atribuir arbitrariamente a localização dos 
polos do sistema, mesmo se o sistema é completamente controlável e 


totalmente observável. 


% É possível atribuir arbitrariamente os polos para controladores 
dinâmicos usando realimentação das saídas, caso o sistema seja pelo 


menos estabilizável e detectável. 
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E ALOCACAO DE POLOS USANDO 


EE» REALIMENTACAO DAS SAÍDAS 


4€ Existem várias abordagens disponíveis para o projeto de controladores 
dinâmicos: 
+ Estima-se o valor do estado usando os sinais de entrada e saída 


do sistema; 


* Utiliza-se esse estado estimado para fazer a realimentação de 


estados. 
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E RASTREAMENTO (TRACKING) 


« А estrutura de alocação de polos mostrada permite que o sistema vá para 


um estado zero, partindo de um estado inicial diferente de zero. 


% Um problema mais geral é quando se deseja que o sistema siga uma 


determinada referência r(k). 


% Este problema pode ser resolvido com um esquema de controle com 
DOF, pois é necessário determinar duas matrizes: matriz de ganho de 


realimentacáo (K) e matriz de ganho de referência (F). 


r(k) v(k)  u(k) x(k+1) x(k) y(k) 
r| oa e 


/Rejane de Barros Araújo 


E RASTREAMENTO (TRACKING) 


% A entrada de referência é dada por v(k) = Fr(k), e a lei de controle é 


escolhida como: 
u(k) = —Kx(k) + Fr(k) 
% onde r(k) é a referência a ser rastreada. 
« As equações do sistema em malha fechada são: 
x(k + 1) = Gax(k) + HEr(k) 
y(k) = Cx(k) 


Ф, 


* com a matriz de estados em malha fechada sendo igual a 


Ga = G — HK 
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E RASTREAMENTO (TRACKING) 


4 Atransformada Z da saída é dada por: 


Y(z) = Clzl, — Ga] ! HFR(z) 


% О erro de rastreamento para uma entrada degrau unitário é dado por: 


lim(z — 1)(Y(z) — R2) = lim(C[zI — Ga] ^ HF — П 


= CU, — Ga] HF — I 


Ф, 


% Para um erro de regime nulo, é necessária a condição: 
Clh — Ga] ^ HF = I, 


4$ Se o sistema for quadrado (nº de entradas = nº de saídas) e G, for estável, 
acha-se a matriz de ganho de referência como: 


F = [C T, — Gae) ill 
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- RASTREAMENTO (TRACKING) 


WE EXEMPLO 


4 Considere o modelo de um motor DC, discretizado com um período de 


amostragem de 20ms, que é representado no espaco de estados por: 


ms = [1799 Es d x + [201563] q) 


x (k+1)] L1 0 x2(k) 


x4 (k) 


y(k) = [0.01191 0.01107] ns 


4$ Projete um controlador usando alocação de polos, tal que o sistema 


tenha as seguintes especificacóes de desempenho: 
** a) Erro de regime permanente nulo para uma entrada degrau; 
* b) Coeficiente de amortecimento de 0,7 e tempo de acomodação em torno 
de 1s. 
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- RASTREAMENTO (TRACKING) 


E» EXEMPLO 


“ Solução: 
4$ % Especificações 
** csi = 0.7; 
« ts = 0.85; Ytempo de acomodação 
* Mp = exp(-csi"pi/sqrt(1-csi^2)) 
« Wn = 4/(csi*ts) 


% Ts = 0.02; % período de amostragem 


4$ % Pólos discretos desejados 
* pdes = exp(-csi Wn*Ts + j"Wn'*sqrt(1-csi^2)* Ts); 
* P = [pdes conj(pdes)] 
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- RASTREAMENTO (TRACKING) 


BEP EXEMPLO 


“ Solução: 
+ % Sistema dicreto no espaço de estados 
4$ G = [1.799 -0.8025;1 0]; 
< H = [0.01563 0]; 
* С = [0.01191 0.01107]; 
% D = 0; 
% sysd = ss(G,H,C,D,Ts) 


4$ % Matriz de realimentação de estados 
% K = acker(G,H,P) 
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E. RASTREAMENTO (TRACKING) 


BED EXEMPLO 


% Solução: 
4 96 Matriz em malha fechada 
% QGcl = G-H*K 


* 96 Matriz de ganho de referência 
% F=inv(C*inv(eye(max(size(G)))-Gcl)*H) 


* 96 Sistema em malha fechada 
% sysdmf = ss(Gcl,H*F,C,D,Ts) 
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E. RASTREAMENTO ( TRACKING) 


EF» EXEMPLO 


% Solução: Resposta ao degrau do sistema em malha fechada 


Step Response 


System: sysdmf 
Peak amplitude: 1.05 
Overshoot (95): 4.6 System: sysdmf 


At time (sec): 0.66 Setting Time (sec): 0.889 


Amplitude 


Time (sec) 
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- CONTROLE INTEGRAL 


88 


« O método de rastreamento mostrado anteriormente é uma ação de 


controle feedforward. 
% Este tipo de controle não usa qualquer tipo de realimentacáo: 
4$ Estratégia esteja sujeita a erros de modelagem (que sempre 
ocorrem na prática); 


% O erro em regime permanente não será nulo. 


* Para evitar esse problema, pode-se utilizar um controle integral. 


Eng. de Conirole e Automação — André Ferreira/Rejane de Barros Araújo 


- CONTROLE INTEGRAL 


* Incluindo o controle integral, o esquema de controle pode ser 


representado como mostrado abaixo: 


(К)  x(k+1) 
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E CONTROLE INTEGRAL 


« As equações do sistema passam a ser: 
x(k + 1) = Gx(k) + Hu(k) 
x(k + 1) = x(k) + r(k) — y(k) 
y(k) = Cx(k) 
u(k) = —Kx(k) — Kx(k) 


4 Рага um sistema com l saídas, o estado x tem dimensão l x 1. 


* A equação de estados pode ser reescritas em função do vetor de 


estados aumentado: 


&(k) = [x(k) x(k)]' 
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E CONTROLE INTEGRAL 


4 Desta forma: 


оз а= Le | - lol. mimo) + res 
уб =1с als 


« Arrumando as equações, tem-se: 
_ n 0 
#(К + 1) = (6 — HK)X(k) + | [ræ 
t 
y(k) = [C 0]Х(Е) 


* onde 
6-|^ | й = [5] К=[к К] 
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- CONTROLE INTEGRAL 


EDS EXEMPLO - SOLUCAO 


% Resolva o exemplo anterior usando controle integral. 


4* As matrizes de estado são: 


p ws K 


H = s C = [0.01191 0.01107] 


% Incluindo o controle integral, obtém-se: 


1.799 —08025 0 0.01563 
б=|©. = 1 0 0 Bee 0 | 
—0.01191 —0.01107 1 0 
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- CONTROLE INTEGRAL 


EFD EXEMPLO - SOLUÇÃO 


4$ % Especificações 
* csi= 0.7; 
* ts = 0.85; Ytempo de acomodação 
* Mp = exp(-csi"pi/sqrt(1-csi^2)) 
« Wn = 4/(csi*ts) 


$* Ts = 0.02; % período de amostragem 


% % Sistema discreto no espaço de estados 
% G=[1.799 -0.8025;1 0]; 
* H = [0.01563 0]; 
* C = [0.01191 0.01107]; 
% D = 0; 
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- CONTROLE INTEGRAL 


EE» EXEMPLO - SOLUÇÃO 


% 9% Inclusão do controle integral 
4 n=max(size(G)); 
$ ОШ = [G zeros(n,1);-C eye(1)]; 
Ні = [H;zeros(1,1)]; 


( 
** Ctil = [C zeros(1)]; 
> Dtil = D; 


Eng. de Controle e Automação — André Ferreira/Rejane de Barros Araújo 


- CONTROLE INTEGRAL 


ES» EXEMPLO - SOLUCAO 


4$ 9e Polos desejados 
* % 1) A partir da especificações 
* pdes = exp(-csi Wn*Ts + j"Wn'*sqrt(1-csi^2)* Ts); 
« 96 2) Devido ao integrador (colocado em 0.2 para ter um efeito 
desprezível na resposta do sistema) 
% pint = 0.2; 
* Р = [pdes conj(pdes) pint] 


4$ % Matriz de realimentação de estados 
4$ КШ = acker(Gtil, Htil,P) 
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- CONTROLE INTEGRAL 


ES EXEMPLO - SOLUCAO 


4$ 96 Matriz em malha fechada 
+ QGcl = Gtil - Htil*Ktil 


> Hcl = [zeros(n,1);eye(1)] 


+ 96 Sistema em malha fechada 
% sysdmf = ss(Gcl,Hcl,Ctil,Dtil, Ts) 
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= CONTROLE INTEGRAL 


>» EXEMPLO - SOLUCÁO 


% Resposta ao degrau do sistema em malha fechada 


step Response 


System: sysdmf 
Peak amplitude: 1.05 
Overshoot (95): 4.55 
: : System: sysdmft 
AE CES Setting Time (sec): 0.915 


Amplitude 


Time (sec 
(sec) g de Barros Araújo 


E ESTIMAÇÃO DE ESTADOS 


« A medição do vetor de estados completo é muitas das vezes impossível 


ou inviável economicamente. 


“ Para implementar um controle usando realimentacáo de estados, uma 


estimativa do vetor de estados X(k) pode ser utilizada. 


« O vetor de estados pode ser estimado a partir de sinais de entrada e 
saída do sistema, utilizando-se um estimador de estados, também 


conhecido como observador de estados. 
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- OBSERVADOR DE ORDEM COMPLETA 


* Para estimar todos os estados de um sistema, uma possibilidade seria 


usar as mesmas equacóes de estado da planta a ser observada. 
% Usando o sistema em malha aberta, tem-se: 
X(k + 1) = GX(k) + Hu(k) 
4€ Este estimador de malha aberta presume um perfeito conhecimento da 


dinámica do sistema e não possui realimentacáo para corrigir os erros 


inevitáveis que ocorrem em qualquer implementação. 
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- OBSERVADOR DE ORDEM COMPLETA 


« As limitações desse tipo de observador ficam claras quando se calcula o 


у= 


£(k + 1) = GX(k) 


erro de estimação: 


% 1) O erro de estimação depende da matriz de estados do sistema e não 


pode ser escolhido arbitrariamente; 


“ 2) Para um sistema instável, o observador também será instável, sendo 


incapaz de rastrear os estados do sistema. 
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- OBSERVADOR DE ORDEM COMPLETA 


% Uma alternativa prática é realimentar a diferença entre a saída medida e 


a saída estimada, como mostrado abaixo (observador de Luenberger): 


£(k + 1) = G&(k) + Hu(k) + L[y(k) — C2(k)] 


% Desta forma, a dinámica do erro de estimação é dada por: 


Z(k + 1) = (G — LC)&(k) 


« Assim, a dinámica do erro de estimação depende dos autovalores da 


matriz do observador, que é dada por: 


Go = G — LC 
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E OBSERVADOR DE ORDEM COMPLETA 


% Diagrama de blocos do observador de ordem completa. 


x(k + 1) = Gx(k) + Hu(k) 
y(k) = Cx(k) 
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E OBSERVADOR DE ORDEM COMPLETA 


* Pode-se mostrar que os autovalores da matriz transposta 
T 
Gy = G! —C'L 


% possui os mesmos autovalores da matriz do observador. 
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E ESTIMAÇÃO DE ESTADOS 


X(k + 1) = GX(k) + Hu(k) + L]y(k) — CX(k)] 
4$ Definição: 
% Se o par (G,C) é observável, então existe uma matriz de ganho de 
realimentação L que posiciona arbitrariamente os polos do 


observador em qualquer conjunto (li, i = 1,2, ...,n]. 


% Se o par (G,C) for apenas detectável, então todos os modos 


observáveis podem ser alocados arbitrariamente. 
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E ESTIMAÇÃO DE ESTADOS 


4 Observação: 
% 1) O par (G,C) é observável (detectável) se, e somente se, o par 


(GT, СТ) for controlável (estabilizável). 


* 2) Pode-se utilizar o MATLAB para resolver o problema da alocação 


de polos do observador usando: 


L = place(G', C',poles)' 
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E. ESTIMACÁO DE ESTADOS 


106 4 EXEMPLO 


€ Determinado sistema discreto (Т = 0,01 s) é representado pelas 


seguintes matrizes de estado: 


1.0 0.1 0.0 1.622x10 ° 
G = É 0.9995 00095 Н = |4.821x10 * C-[1 0 0] 
0.0 —0.0947 0.8954 9.468х1072 


% Calcule a matriz de ganhos L, рага que os polos do observador estejam 


localizados em (0.1, 0.2 + j0.2]. 
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E. ESTIMACÁO DE ESTADOS 


ERTA EXEMPLO - SOLUCAO 


% >> G = [1 0.1 0;0 0.9995 0.0095;0 -0.0947 0.8954] 
4 >> Н = [1.622е-6 4.821e-4 9.468e-2]' 

% >>C=[100] 

4$ >> рде = [0.1 0.24j*0.2 0.2-j*0.2] 

% >> L = place(G',C',pdes)' 


* >> eig(G-L*C) 
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Е. OBSERVADOR DE ORDEM COMPLETA 


O observador que utiliza a expressáo 
£(k + 1) = G£(k) + Hu(k) + L[y(k) — C&(k)] 
% é chamado de observador preditor (prediction observen. 

% Vetor de estados estimados (e ações de controle associadas) em 
determinado instante de amostragem não dependem do valor de medida 
atual da saída do sistema. 

% Um outro tipo de observador, chamado de filtering observer, estima о 


vetor de estados utilizando o valor da saída atual do sistema (supondo 


que o tempo de processamento é desprezível). 
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Е. OBSERVADOR DE ORDEM COMPLETA 


4$ Desta forma: 
£(k + 1) = GX(k) + Hu(k) + Ll[y(k + 1) — C(GX(k) + Hu(k))] 
% A dinâmica do erro de estimação neste caso é dada por: 


Z(k + 1) = (G — LCG)X() 


* Comparando os dois tipos de observadores, pode-se concluir que para 
o filtering observer a matriz C passa a ser substituída pelo produto CA. 


< Nova matriz de observabilidade: CA 
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- OBSERVADOR DE ORDEM COMPLETA 


4$ Se o par (G,C) é observável, о par (С, Сс) também será observável, 


exceto quando a matriz A possuir autovalores nulos. 


“ Neste caso, o раг (С, CG) será detectável, porque os autovalores nulos 


estáo associados com modos estáveis. 


% Além disto, autovalores nulos estão associados com os modos mais 
rápidos: 
% Projeto do observador pode ser feito escolhendo-se uma matriz L que 


forneca valores adequados para os autovalores restantes de (G — LCG). 
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- OBSERVADOR DE ORDEM COMPLETA 


EE EXEMPLO 


* Рага o exemplo anterior, determine o valor da matriz de ganhos do 


filtering observer. 


1.0 0.1 0.0 1.622х10-6 
a= [oo 0.9995 00095 Н = |4.821х107* С= [1 0 0] 
0.0 —0.0947 0.8954 9.468х1072 


% Solução: 
* Usando a mesma localização dos polos definida anteriormente: 
% >> L = place(G',(C*G)';pdes) 
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Е. OBSERVADOR DE ORDEM REDUZIDA 


Observador de ordem completa: 
% O vetor de estados é estimado a partir das entradas e saídas do 


sistema. 


€ Para um sistema de ordem n e com l saídas (que são funções lineares 
dos estados), porque estimar n variáveis? 
> É possível estimar n—1 variáveis e a partir dessas estimativas 


reconstruir o vetor de estados completo. 


“ Desta forma, para sistemas SISO o observador de ordem reduzida terá 


ordem n — 1. 
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- OBSERVADOR DE ORDEM REDUZIDA 


4 Suposicoes: 


« Matriz de entrada B e matriz de saída C possuem rank completo. 


% Desta forma, os elementos do vetor de saída y são linearmente 
independentes, formando um vetor de comprimento l, deixando n — l 


variáveis a serem determinadas. 


=} 


* onde M é uma matriz (n — Dxn de rank completo, cujas linhas são 


4 Fazendo: 


linearmente independentes de C, e zé o estado parcial desconhecido. 
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E OBSERVADOR DE ORDEM REDUZIDA 


% As matrizes em espaço de estados para as variáveis de estado 


transformadas são: 


B, i 4d 
В, = Q; B= ES 
2 in-i 
A = 07140, = паті 
AO 0 A; | As tn-l 
—— 
| n-i 
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- OBSERVADOR DE ORDEM REDUZIDA 


« A equação de estados para o estado parcial desconhecido е: 
4 Define-se também a variável de saída: 
Ya (k) = y(k + 1) — Gy(k) — H,u(k) = +G2z(k) 
* A dinâmica do observador, já incluindo o erro no cálculo de y, é 
assumida como sendo linear e invariante no tempo, sendo dada por: 
2(k + 1) = Gsy(k) + G42(k) + H;u(k) + L[y, (k) — G22(k)] 
= (G, — 1б,)#(®) + бзу(®) + L(y(k + 1)- Gy (k) — Нуш) 1 + Н;щ(®) 
% onde ¿Z representa o valor estimado para o vetor de estados parcial z. 
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E OBSERVADOR DE ORDEM REDUZIDA 


* Para evitar ter que usar o termo y(k + 1), estima-se a seguinte variável: 
x(k) = Z(k) — Ly(k) 


“ Desta maneira, obtém-se o observador de ordem reduzida: 
X(k + 1) = (G, — LG;)x(k) + (G4L + G4 — LG, — LG;L)y(k) + (H, — LH Ju(k) 
= GoX(k) + Gyy(k) + Hou(k) 
% onde: 
Go = G4 — LG; 
Gy = GoL + 63 — LG, 


Ho — H, — LH, 
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ый OBSERVADOR ОЕ ORDEM REDUZIDA 


* Diagrama de blocos do observador de ordem reduzida: 


x(k +1) x(k) 


x(k + 1) = Gx(k) + Hu(k) 
y(k) = Cx(k) 
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Е. OBSERVADOR DE ORDEM REDUZIDA 


“ A dinámica do observador de ordem reduzida depende da matriz A,, 
cujos autovalores devem: 
« 1) Estar dentro do círculo; 
* 2) Ser suficientemente rápidos, para rastrear os estados do sistema 


observado. 
O projeto do observador de ordem reduzida depende basicamente da 


matriz L, já que a partir desta matriz todas as outras matrizes 


necessárias podem ser obtidas. 
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E OBSERVADOR DE ORDEM REDUZIDA 


% Assim, o vetor de estados pode ser calculado através de: 


R 
х0) 0 У | 
Z(R) 


= š k De Мер T bes 
“rh dx) “|х(®) 


% lembrando que: 
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= OBSERVADOR DE ORDEM REDUZIDA 


NiTTUTOFID 
ми 


ED OBSERVAÇÃO 


4 1) Fazendo б = б, = Г, pode-se mostrar que os polos de 4,7 


podem ser alocados arbitrariamente desde que o par (6,, G> ) seja 


controlável ou, na forma dual, que o par (G,, G2) seja observável. 


4 2) Pode-se demonstrar que para o par (G,,G,) ser observável, o par 


(G, C) também deve ser observável. 
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= OBSERVADOR DE ORDEM REDUZIDA 


NiTTUTOFID 
ми 


ШЕ EXEMPLO 


* Рага о sistema descrito nos exemplos anteriores, projetar um observador 


de ordem reduzida. 


1.0 0.1 0.0 1.622x10 ° 
G = É 0.9995 00095 Н = |4.821x107* C=[1 0 0] 
0.0 —0.0947 0.8954 9.468х1072 


% Solução: 


« Ver arquivo “10ргој observador ordem reduzida.m" 
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- PROJETO DO REGULADOR: COMBINACAO DA REALIMENTACAO 


ED DE ESTADOS COM O ESTIMADOR DE ESTADOS 


r(k) (К) х(К+1) x(k) 


OBSERVADOR 


REALIMENTAÇÃO DE ESTADOS 
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- REALIMENTACAO DE ESTADOS DO 


IED OBSERVADOR 


* Caso o vetor de estados não esteja disponível para ser realimentado, 


Y 


pode-se utilizar os estados estimados gerando a seguinte lei de controle: 


u(k) = —KX(k) + v(k) 


* Substituindo em: 
x(k + 1) = Gx(k) + Hu(k) 


y(k) = Cx(k) 
* Obtém-se: 
x(k + 1) = Gx(k) — HKX(k) + Hv(k) 
* Reescrevendo na forma do erro de estimação X = x — £: 


x(k + 1) = (G — HK)x(k) + HKX(k) + Hv(k) 
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- REALIMENTACAO DE ESTADOS DO 


ERP» OBSERVADOR 


% Lembrando que para um observador de ordem completa é válida a 
relação: 


X(k + 1) = (G — LC)X(k) 


« Assim, tem-se que: 


k k 
жут) 200] * [o] veo 


Polinômio característico: 
Aa(A) = det|AI — (G — HK)]|det|AI — (G — LC)] = А„(Л4)А— (A) 
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- PRINCÍPIO DA SEPARACAO 


% Os autovalores do sistema em malha fechada do sistema podem ser 
escolhidos de forma completamente independente dos autovalores do 
observador. 

4$ Controlador e observador são projetados separadamente е 


utilizados em conjunto! 
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E ESCOLHA DOS AUTOVALORES DO 


126 OBSERVADOR 


% A escolha dos pólos do observador não se baseia nas restrições para a 
acáo de controle que existe na escolha dos pólos do controlador. 


% Polos do Controlador ^ devem atender às especificações de 
desempenho desejadas. 


4$ Polos do Observador > como uma regra empírica, devem ser de З а 10 
vezes mais rápidos que os pólos do controlador. 

% Na prática, deve-se tomar cuidado para não escolher polos 
extremamente rápidos, pois assim o observador irá estimar os ruídos 
de medição, ao invés dos estados do sistema. 

* Evitar projetar um observador deadbeat! 
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EL. ESCOLHA DOS AUTOVALORES DO OBSERVADOR 


ERES EXEMPLO 


* Para o sistema descrito nos exemplos anteriores, fazer uma 
realimentação de estados usando um observador de ordem compleia, 
de tal forma que: 

*• Polos do sistema em MF: (0.6, 0.4 + j;0.33) 
% Estado inicial do sistema: x(0) = [11 1]? 


« Estado inicial do observador: X(0) = [0 0 0]”. 


1.0 0.1 0.0 1.622x1078 
G=|0.0 0.9995 00095| H = |4.821x10 ^ C=[1 0 0] 
0.0 —0.0947 0.8954 9.468х1072 
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Hi. ESCOLHA DOS AUTOVALORES DO OBSERVADOR 


ED EXEMPLO - SOLUCAO 


+ % Sistema discreto no espaço de estados 
* G = [1 0.1 0;0 0.9995 0.0095;0 -0.0947 0.8954]; 
Н = [1.622е-6 4.821e-4 9.468e-2]'; 
% C = [100]; 
% D = 0; 
« Ts = 0.01; % período de amostragem 


4$ n = max(size(G)); % ordem do sistema 


“+ 9o Pólos desejados рага o sistema em MF 
4 pdes c = [0.6 0.44j*0.33 0.4-j*0.33]; 
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К. __ ESCOLHA DOS AUTOVALORES DO OBSERVADOR 


ED EXEMPLO - SOLUCAO 


% % Polos desejados para o observador 
% pdes_o = [0.1 0.1+j*0.1 0.1-j*0.1]; 


+ % Matriz de ganhos para realimentação de estados 


* K = place(G,H,pdes с) 


% % Matriz de ganhos do observador 


% L=place(G',C',pdes o) 


% % Sistema Controlador + Observador 
* GG = [G-H*K H*K;zeros(n) G-L*C]; 
* HH = [H;zeros(n,1)]; 
* CC = [C zeros(1,n)]; DD = 0; 
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Hi. ESCOLHA DOS AUTOVALORES DO OBSERVADOR 


ED EXEMPLO - SOLUCAO 


% sys = ss(GG,HH,CC,DD,Ts) 


* 9oEstado inicial do sistema 
$ x0=[111]; 


* 96 Estado inicial do observador 
$* x0o = [0 0 0]; 


* % Estado inicial total 
% xx0 = [x0;x0-x00] % Estado e erro de estimação 
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E... ESCOLHA DOS AUTOVALORES DO OBSERVADOR 


ED EXEMPLO - SOLUCAO 


« Resposta livre ao estado inicial: 
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Hi. ESCOLHA DOS AUTOVALORES DO OBSERVADOR 


ED EXEMPLO - SOLUCAO 


“ %Rastreamento da referência: 
4$ % Matriz de ganho de referência 
4$ F = inv(C'inv(eye(max(size(GA)))-(G-H*K))*H) 


4 HHr = [H*F;zeros(n,1)] 


+ 96 Sistema em malha fechada 
4» sysmf = ss(GG,HHr,CC,DD,Ts) 
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E... ESCOLHA DOS AUTOVALORES DO OBSERVADOR 


ED EXEMPLO - SOLUCAO 


« Resposta ao degrau: 


Step Response 
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Time (sec) 
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ASSUNTOS COMPLEMENTARES 


% 1) Invariáncia dos zeros do sistema: 

“* Uma limitação da realimentação de estados é a incapacidade desse 
tipo de controle em modificar os zeros do sistema em malha aberta, o 
que pode afetar significativamente a resposta do sistema em malha 
fechada (Cf. Fadali, capítulo 9, seção 9.4, p.353). 


* 2) Alocação de polos usando função de transferência: 
* Envolve a resolução de uma equação diofantina; 
« Método bastante utilizado para projetar controlador a partir de 
modelos paramétricos de sistemas reais. 
(Cf. Fadali, capítulo 9, seção 9.7, p.370). 
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